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I.  INTRODUCTION 


This  report  covers  two  separate  efforts  conducted  by  the  author 
under  AFOSR  Contract  F44620-76-C-0023  from  October  1,  1975  to  September 
30,  1977.  The  first  and  larger  of  these  efforts  was  an  analytical  and 
experimental  study  of  the  erosive  burning  of  composite  propellants,  while 
Lhe  second  was  an  analytical  modeling  study  of  the  combustion  of  aluminum 
particles  in  nitrogen-carbon  dioxide  mixtures. 

Erosive  burning,  the  augmentation  of  solid  propellant  burning 
rate  by  the  flow  of  products  across  a burning  surface,  is  becoming  in- 
creasingly important  with  use  of  lower  port-to-throat  area  ratio  motors 
and  nozzleless  motors  which  result  in  high  velocity  crossflows.  The 
response  of  various  propellants  to  such  crossflows  must  be  known  by  the 
motor  designer  in  order  for  him  to  perform  adequate  motor  design.  In 
addition,  it  is  important  ^?at  the  propellant  formulator  understand  the 
effect  of  various  formulation  parameters  on  the  sensitivity  of  a propellant 
to  crossflows  ;o  that  he  may  tailor  his  propellants  to  the  desired  charac- 
teristics. (For  example,  in  a nozzleless  rocket  motor,  the  decrease  in 
pressure  from  the  head  end  to  the  aft  end  of  the  grain  tends  to  result  in 
slower  burning  at  the  aft  end  in  the  absence  of  erosive  effects.  Depending 
upon  the  sensitivity  of  the  formulation  to  crossflow,  the  increasing  Mach 
Number  along  the  grain  port  may  lead  to  undercompensation,  exact  cancella- 
tion, or  over compensation  of  the  pressure  effect.)  In  this  program,  a 
test  apparatus  was  constructed  for  the  study  of  the  effects  of  crossflows 
up  to  Mach  1 on  propellant  ballistics,  and  seven  propellants  (with  system- 
atically varied  formulation  parameters)  were  characterized  in  this  apparatus. 
In  addition,  a first  generation  analytical  model  for  erosive  burning  of 
composite  propellants  based  on  bending  of  columnar  diffusion  flames  by 
crossflow  was  developed  and  considerable  progress  toward  development  of 
a more  fundamental  second  generation  model  (still  based  on  flame-bending) 
was  made.  Results  of  the  experimental  and  analytical  efforts  to  date 
are  described  in  Section  II.  This  work  is  being  continued  under  AFOSR 
Contract  F49620-78-C-0016  in  fiscal  year  1978. 
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Existing  published  models  of  the  combustion  of  aluminum  particles 
employ  a flame-sheet  approximation  wherein  the  reaction  of  aluminum  vapor 
and  oxidizer  is  assumed  to  occur  instantaneously  in  an  infinitesimally 
thin  reaction  zone  located  some  distance  from  the  particle  surface,  with 
the  combustion  thus  being  controlled  by  the  diffusion  of  oxidizer  and 
fuel  species  to  this  zone.  Preliminary  calculations  (neglecting  condensa- 
tion effect)  of  aluminum  particle  combustion  with  consideration  of  finite- 
rate  kinetics  using  kinetic  data  recently  measured  by  Fontijn  for  the 
reaction  of  aluminum  vapor  with  carbon  dioxide  indicated  that  the  flame- 
sheet  approximation  is  quite  poor  for  aluminum  particles  smaller  than 
approximately  30  to  50  microns  in  diameter,  a size  regime  of  considerable 
interest  regarding  aluminum  combustion  in  solid  propellant  rocket  motors. 

In  fact,  this  analysis  indicated  that  aluminum  particle  combustion  more 

115  2 

closely  followed  a d - to  d ' -burning  rate  law  than  the  d -burning  rate 

law  which  follows  from  the  assumption  of  infinite  kinetics.  Since  most 
laboratory  burning  rate  data  available  for  aluminum  particle  r.^bustion 
are  for  particles  of  50  microns  and  larger  diameter  while  aluminum  cast  in 
motors  is  typically  in  the  5 to  20  micron  diameter  range,  calculation  of 
degree  of  combustion  of  aluminum  in  motors  is  quite  sensitive  to  the  power 
law  used  to  extrapolate  the  laboratory  data  to  motor  conditions.  Accord- 
ingly, a model  allowing  for  finite-rate  kinetics  for  the  initial  aluminum 
vapor-oxidizer  reaction  and  also  incorporating  a postulated  description 
of  aluminum  oxide  condensation  was  developed  for  further  examination  of 
the  sensitivity  of  predicted  burning  rate  (and  burning  rate-particle 
diameter  dependency)  to  the  assumption  of  infinite  kinetics.  Model  devel- 
opment results  are  presented  in  Section  III. 
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II.  EROSIVE  BURNING  OF  COMPOSITE  PROPELLANTS 
A.  Introduction  and  Background 

Erosive  burning  refers  to  the  augmentation  of  "normal"  burning 
rate  of  a solid  propellant  by  action  of  a flow  of  gas  parallel  to  the  burning 
surface.  Several  empirical  relations  appear  in  the  literature  to  describe 
the  total  burning  rate,  which  is  the  sum  of  the  normal  and  erosive  rates. 


Data  are  often  correlated  by  one  of  the  following  expressions: 

e - r/r  - 1 + K.  (V  - V )m,  m < 1 (1) 

o 1 t — 

e = r/r  - 1 + K,(M  - Mjm,  m < 1 (2) 

o l t — 

e - r/r  «=  1 + K,(G  - Gjm,  m < 1 (3) 

o J t — 

e = r/r  = 1.  + (aG°‘8L  °’1 2/r  ) exp  (-8r  p /G)  (4) 

o o ° P 


e “ erosive  burning  rate/burning  rate  at  same  pressure 
in  the  absence  of  crossflow 

r ■ propellant  linear  burning  rate,  including  erosive 
effects 

r * propellant  linear  burning  rate  in  the  absence  of 
° crossflow 

V - mainstream  crossflow  velocity 
M - mainstream  crossflow  Mach  Number 
G ■ crossflow  mass  flux 

L * length  parameter  defining  crossflow  Reynolds'  Number 
a,  6,  K^,  1^,  K3,  m ■ empirical  constants 

the  fourth  of  these  being  a form  of  the  Lenoir-Robillard  expression.  The 
subscript  "t"  refers  to  threshold  crossflow  conditions  below  which  erosion 
does  not  occur.  (Some  propellants  have  been  correlated  with  non-zero 
threshold  values,  while  others  have  been  correlated  with  threshold  values 
set  equal  to  zero.) 

General  observations  of  importance  from  the  past  experimental 
studies  (1-10)  include: 

1.  Plots  of  burning  rate  versus  gas  velocity  or  mass  flux  at  constant 

pressure  are  usually  not  fitted  best  by  a straight  line. 
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2.  Threshold  velocities  and  "negative"  erosion  rates  are  often 
observed. 

3.  Slower  burning  propellants  are  more  strongly  affected  by  cross- 
flows  than  higher  burning-rate  formulations. 

4.  At  high  pressure,  the  burning  rate  under  erosive  conditions  tends 
to  approach  the  same  value  for  all  propellants  (at  the  same  flow 
velocity)  regardless  of  the  burning  rate  of  the  propellants  at 
zero  crossflow. 

5.  Erosive  burning  rates  do  not  depend  upon  gas  temperature  of  the 
crossflow  (determined  from  tests  in  which  various  "driver  propell- 
ant's" products  are  flowed  across  a given  test  propellant). 

There  is,  however,  very  little  data  available  for  high  crossflow  velocities 
(greater  than  M ^ 0.3).  In  addition,  there  has  been  no  study  in  which 
various  propellant  parameters  have  been  systematically  varied  one  at  a 
time.  Such  a study  is  necessary  for  determination  of  erosive  burning 
mechanisms  and  proper  modeling  of  the  erosive  burning  phenomena.  Much  of 
the  past  work  has  not  resulted  in  Instantaneous  (as  opposed  to  averaged 
over  a range  of  pressure  and  crossflow  velocity)  measurements  of  erosive 
burning  rates  under  well-characterized  local  flow  conditions. 

Erosive  burning  can  have  a strong  influence  on  rocket  motor 
performance  (ballistics)  as  reflected,  for  example,  in  pressure-time  curves 
for  low  port-to-throat  area  ratio  motor  designs.  Such  effects  are  parti- 
cularly important  in  the  case  of  nozzleless  rocket  motors,  where  the  gas 
flow  is  choked  near  or  at  the  aft  end  of  the  grain  through  most  or  all  of 
the  motor  operation.  Since  this  point  is  thus  the  effective  throat,  and 
the  throat  area  is  thus  a function  of  regression  rate  of  a propellant 
surface  being  subjected  to  a very  high  crossflow  velocity,  the  result 
is  a chamber  pressure-time  history  which  is  very  dependent  upon  the  erosive 
burning  characteristics  of  the  propellant.  The  effects  of  erosive  burning 
on  solid  propellant  rocket  interior  ballistics  for  low  port-to-throat 
area  ratio  motors  and  nozzleless  motors  are  discussed  in  some  detail  by 
this  author  in  Reference  11.  As  an  example,  results  of  an  analysis  of 
the  ballistics  of  a nozzleless  motor  with  initially  uniform  port  area, 
shortly  after  ignition,  with  an  assumed  no-crossflow  burning  rate  pressure 
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relationship  of  rQ  - bp°  and  an  erosivity  relationship  of  e * 1 + ^M, 

are  presented  in  Table  II-I.  Values  of  (r  , / r„  )._,.,  are  presented 

as  a function  of  the  erosivity  constant  (Kj)  and  the  burning  rate  exponent 

(n) . As  may  be  seen,  for  the  case  of  no  erosion  (Kj  * 0),  the  aft  end  will 

recede  more  slowly  than  the  fore  end,  due  to  lower  pressure  at  the  aft  end. 

As  K„  increases,  the  r ,./rr  ratio  also  increases,  going  through  unity 
i art  fore 

(generally  desirable)  at  a value  of  Kj  which  depends  on  the  burning  rate 
exponent.  The  results  of  Table  II-I  give  some  indication  of  the  sensitivity 
of  nozzleless  motor  design  to  the  erosive  burning  characteristics  of  the 
propellant,  and  thus  point  out  the  importance  of  information  regarding  the 
propellant's  erosive  burning  characteristics  to  the  designer  and  the 
propellant  formulator. 

Since  there  is  such  a strong  interaction  between  the  local  flow 
environment  and  the  propellant  burning  rate,  it  is  necessary  to  be  able 
to  predict  this  interaction  in  order  to  design  and  calculate  the  performance 
of  a low  port/ throat  area  ratio  rocket  (particularly  a nozzleless  rocket 
with  a port/throat  area  ratio  of  unity).  A review  of  the  literature  has 
indicated  that  there  is  no  unifying  model  or  theory  which  can  be  used  to 
reliably  predict  propellant  burning  rates  in  an  erosive  situation,  nor  is 
there  a supply  of  systematic  experimental  data  characterizing  the  erosive 
burning  behavior  of  propellants  as  a funciton  of  compositional  variables. 
Thus,  development  of  an  analytical  model  of  erosive  burning  properly  des- 
cribing the  physical  effects  which  result  in  crossflow  velocities  augmenting 
solid  propellant  burning  rate,  coupled  with  an  experimental  program  to 
define  the  effects  of  various  parameters  on  erosive  burning,  is  important 
to  the  design  and  development  of  advanced  solid  rocket  systems. 

B.  Analytical  Modeling 

1.  Review  of  Past  Work 

Over  the  years,  a large  number  of  models  of  erosive  burning  of 
composite  (heterogeneous)  and  double-base  (homogeneous)  propellants  have 
been  developed.  A list  of  models  by  author's  name,  divided  into  four 
categories  of  models,  is  presented  as  Table  II-II.  These  models  have  been 
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TABLE  II-I.  Zero-Time  Ballistic  Analysis  of  Nozzleless 
Motor  with  Initially  Uniform  Port  Area 


e = r/r  ■ 1 + K..M 
o Z 

FORE  END:  M **  0 


rFORE  * ro,FORE  " bP  STATIC  FORE 


AFT  END:  M = 1 


rAFT  “ (K2  + 1)ro,  AFT  " <K2  + 1)bP  STATIC  AFT 


PSTATIC  AFT/PSTATIC  FORE  * 1/(y  + 
- (K2  + 1)b  pn 

' • rAFT  , . .,n  STATIC  FORE 

(Y  + 1) 


THEREFORE  rAFT/rF0RE  ^ <K2  + D/(Y  + D 


FOR  Y * 1*25 


n 

0.4 


0.6 


0.8 


K r/r 

2 AFT  FORE 


0 

0.72 

0.5 

1.08 

1.0 

1.45 

1.5 

1.80 

0 

0.61 

0.5 

0.92 

1.0 

1.23 

1.5 

1.54 

0 

0.52 

0.5 

0.78 

1.0 

1.05 

1.5 

1.31 
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TABLE  II-II.  Erosive  Burning  Models  by  Category 


1.  Models  Based  on  Heat  Transfer  from  a "Core  Gas"  in  the 
Presence  of  Crossflow: 

Lenior  & Robillard  (15) 

Burick  and  Osborn  (16) 

Zucrow,  Osborn  and  Murphy  (17) 

Saderholm  (3) 

Marklund  (7) 

Jojic  & Blagojevic  (18) 

2.  Models  based  on  Alteration  of  Transport  Properties  in  Region 
from  Surface  to  Flame  Zone  by  Crossflow,  Generally  Due  to 
Turbulence  Effects.  Includes  Effects  on  Conductivity  from 
Flame  Zone  Back  to  Propellant  and  Effects  on  Time  for 
Consumption  of  Fuel  Pockets  Leaving  Surface: 

Saderholm,  Biddle,  Caveny,  et  al.  (19) 

Lengelld  (20) 

Corner  (Double-Base)  (21) 

Vandenkerckhove  (Double-Base)  (22) 

Zeldovich  (Double-Base)  (23) 

Vilyunov  (Double-Base)  (10) 

Geckler  (24) 

3.  Models  Based  on  Chemical  Reacting  Boundary  Layer  Theory 
(Homogeneous  Systems  Only) : 

Tsuji  (25) 

Beddini,  et  al.  (26) 

Kuo,  et  al.  (27) 

4.  Other 

Klimov  (18) 

Molnar  (29) 

Miller  (30) 

King  (12,13) 


discussed  by  this  author  In  References  11  - 14,  and  will  not  be  discussed 
In  detail  here.  A general  weakness  of  the  models  in  the  first  category 
(models  based  upon  added  heat  transfer  from  a "core  gas"  In  the  presence  of 
crossflow)  is  that  they  predict  substantial  dependence  of  the  erosive  burning 
contribution  on  the  temperature  of  the  core  gas:  such  dependence  was  found 
by  Narklund  and  Lake  (7)  to  be  absent,  as  discussed  in  References  11  - 13. 

Of  the  models  other  than  the  one  developed  by  this  author  and  discussed 
below,  those  of  Lengelld  (20),  Beddini,  et  al.  (26),  and  Kuo,  et  al.  (27), 
appear  to  be  the  most  advanced  of  the  models  listed.  Common  to  all  three 
of  these  models  is  the  assumption  that  the  Increase  in  propellant  burning 
rate  associated  with  crossflow  results  from  turbulence  associated  with 
this  crossflow  penetrating  between  the  propellant  gas  flame  zone(s)  and 
the  surface,  causing  increases  in  mass  and  energy  transport  rates.  However, 
for  a typical  propellant  containing  oxidizer  particles  with  diameters  of 
from  10  to  50  microns,  diffusion  flame  distances  may  be  calculated  to  be 
typically  of  the  order  of  one-quarter  to  one-half  of  the  particle  diameter, 
or  2.5  to  25  microns.  On  the  other  hand,  for  a crossflow  velocity  of  200 
meters/second  (650  ft/second),  the  universal  u f y correlation  (trans- 
piration effects  neglected)  indicates  a laminar  sublayer  thickness  of 
approximately  10  microns  and  a buffer  zone  thickness  of  about  50  microns, 
full  turbulence  not  being  achieved  closer  than  60  microns  from  the  propellant 
surface.  Moreover,  transpiration  of  the  binder  and  oxidizer  decomposition 
gases  from  the  propellant  surface  will  tend  to  increase  the  thickness  of 
these  zones.  Thus,  it  is  not  at  all  certain  that  crossflow- induced  turbu- 
lence does  penetrate  into  the  zone  between  the  propellant  surface  and  the 
gas-phase  flame  zone(s) . In  addition,  even  if  the  turbulent  region  does 
extend  into  this  zone,  in  order  for  the  eddies  to  have  significant  effect 
on  mixing  and  thus  on  heat  and  mass  transfer,  they  must  be  considerably 
smaller  than  the  flame  offset  distance;  that  is,  they  must  be  on  the  order 
of  one  micron  in  diameter  or  less.  It  is  not  clear  to  this  author  that 
a significant  amount  of  turbulence  of  this  scale  will  be  induced  in  the  zone 
between  the  propellant  surface  and  the  gas-phase  flame  zone(s)  by  crossflows 
up  to  Mach  1,  more  than  an  order  of  magnitude  above  typical  erosive  burning 
threshold  velocities.  Accordingly,  an  alternate  possible  mechanism  for 
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erosive  burning  of  composite  propellants  based  upon  bending  of  columnar 
diffusion  flames  by  a crossflow  has  been  postulated  and  model  developemnt 
based  upon  this  picture  has  been  carried  out  by  this  author. 

In  the  development  of  a proper  model  of  erosive  burning  of  a 
given  class  of  propellants,  it  is  necessary  that  a physical-chemical  mech- 
anism for  the  "normal"  (no  crossflow)  burning  of  such  propellants  be 
specified,  that  the  boundary  layer  flow  be  properly  described  (theoretically 
or  empirically)  and  that  the  description  of  these  processes  be  properly 
coupled. 

Considering  first  the  flow  field,  it  is  informative  to  estimate 

flow  profiles  and  angles  near  the  surface  of  a composite  propellant  for 

a typical  erosive  burning  situation.  As  an  example,  let  us  examine  a case 

6 2 

where  the  operating  pressure  is  6.89*10  N/ni  (1000  psi) , the  propellant 
flame  temperature  is  3000°K,  the  crossflow  mainstream  velocity  is  200  m / 
sec  (650  ft/sec),  the  characteristic  length  dimension  for  determining 
Reynold's  Number  is  15  cm  (0.5  feet)  and  the  propellant  burning  rate  is 
1.25  cm/sec  (0.5  in/sec).  In  this  case,  the  gas  velocity  away  from  the 
surface  calculated  at  the  flame  temperature  Is  approximately  4 m/sec  (13 
ft/sec).  Using  Mickley  and  Davis's  (31)  flow  profile  data  for  boundary 
layer  velocity  profiles  in  the  presence  of  transpiration,  we  estimate 
that  the  crossflow  velocity  10  pm  from  the  propellant  surface  is  about 
10  m/sec  (30  ft/sec).  A simplified  energy  balance  equating  the  heat 
feedback  flux  from  a flame  sheet  above  a propellant  surface  to  the  value 
required  for  preheating  and  vaporizing  the  solid  ingredients  at  a regression 
rate  of  1.25  cm/sec  (0.5  in. sec)  indicates  that  the  gas-phase  flame  must 
be  on  the  order  of  10  pm  from  the  surface.  Thus,  t the  position  of  the 
gas-phase  flame  front,  the  velocity  component  away  from  the  propellant  is 
about  4 m/sec,  while  the  velocity  component  parallel  to  the  surface  is 
10  m/sec,  and  the  resultant  flow  vector  makes  an  angle  with  the  propellant 
surface  of  only  22  degrees.  While  this  vector  will  vary  with  distance 
from  the  surface,  since  the  velocity  components  normal  to  and  parallel 
to  the  surface  do  not  scale  with  distance  from  the  surface  in  exactly 
the  same  way,  the  variation  will  not  be  great.  Thus,  fuel  and  oxidizer 
gas  columns  leaving  the  surface  will  not  flow  perpendicular  to  the  surface 
(as  they  would  in  the  absence  of  crossflow,  but  at  an  angle  of  approxi- 
mately 20  to  25  degrees  from  parallel  with  the  surface  for  this  typical  case. 
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The  important  feature  of  this  picture  is  that  any  diffusion 
flame  at  the  AP-binder  boundaries  is  bent  over  toward  the  propellant 
surface  by  the  crossflow  velocity.  Since  the  deflection  of  this  mixing 
column  or  cone  can  be  shown  to  cause  the  distance  from  the  base  to  the  tip, 
measured  perpendicular  to  the  propellant  surface,  to  decrease,  the  height 
above  the  propellant  at  which  any  given  fraction  of  the  mixing  of  AP  pro- 
ducts and  fuel  decomposition  products  is  complete  should  therefore  be  de- 
creased and  the  distance  from  the  propellant  surface  to  the  "average" 
location  of  the  diffusion  flame  should  also  be  decreased.  This  in  turn 
will  increase  heat  feedback  and  thus  increase  the  burning  rate.  A first 
generation  model  based  upon  this  picture,  which  permits  prediction  of  burning 
rate-pressure-crossflow  velocity  relationships,  given  only  no-crossflow 
burning  rate  versus  pressure  data  over  a wide  pressure  range,  has  been 
developed.  This  model,  described  briefly  below,  has  been  presented  in 
detail  in  Reference  13.  In  addition,  considerable  progress  has  been  made 
on  development  of  a second  generation  model  of  a more  fundamental  nature, 
requiring  as  input  only  propellant  composition  and  particle  size  data: 
formulation  of  this  model  is  also  discussed  below. 

2.  First  Generation  Model 

A schematic  depicting  the  first  generation  composite  propellant 
erosive  burning  model  is  presented  in  Figure  II-l.  In  the  first  part  of 
the  figure,  we  picture  the  flame  processes  occurring  in  the  absence  of 
crossflow.  There  are  two  flames  considered:  an  ammonium  perchlorate 
deflagration  monopropellant  flame  close  to  the  surface;  and  a columnar 
diffusion  flame  resulting  from  mixing  and  combustion  of  the  AP  deflagra- 
tion products  and  fuel  binder  pyrolysis  products  at  an  average  distance 
somewhat  further  from  the  surface.  Three  Important  distance  parameters 
considered  are  the  distance  from  the  propellant  surface  to  the  "average" 
loation  of  the  kinetically  controlled  AP  monopropellant  heat  release  (L^) , 
the  distance  associated  with  mixing  of  the  oxidizer  and  fuel  for  the  diffu- 
sion flame  an<3  the  distance  associated  with  the  fuel-oxidizer 

reaction  time  subsequent  to  mixing  (L^n).  A heat  balance  between  heat 


10 


= <3  “ ¥ 

* §233 

n . 


<r 

Ui 

Q 

Z 

00 

ac 

UJ 

> 

o 

t/l 

UJ 

I/I 

o 

_J 

(_> 

UJ 

2 

u. 

z 

o 

V) 

3 

ll 

u. 

5 

UJ 

tr 

O 

LL 

uj 

(C 

UJ 

z 


o 

UJ 

ft 

3 

r 

u 

CL 

UJ 

3 

o 

z 

a 

- 1 

J 

UJ 

3 


UJ 

!| 

> N 
ff  o 

< z 

o X 
O 2 

a * 


UJ 

ui 

UM 

N 

2 

►- 

a. 

o 

IL 

5 

x 

UJ 

2 

if 

1a 

E 

z 

UJ 

0 

j 

5 

H 

o 

in 

u. 

♦ 

H 

u 

< 

UJ 

a. 

< 

UJ 

J 

UJ 

a; 

= i_ 

UJ 

3 

u. 

u. 

0 

K 

Z 

s 

-1 

l.l 

K 

IN 

Jf 

1 

-J  UJ 

z 


Iz  g 

UJ  - o 

H X ’ 
< 


Q 

Z 

UJ 

u. 

O 


3 Z 

8S  5 

z (7  uj 
< uj  > , 

H £ < < 
i/i  ® 


— — < 
O i/i 


uj  ui 

o o 
cc 

UJ 
> 


a. 

O 

cr 

a. 

O 

z 

o 

2 

a. 

< 


UJ 

> 

3 


(A 

»- 

4 

J- 


H 

■ cr 
a 


+ + * 


u 

3 

UJ 

> 


u! 

A 

o 

oc. 

O 

O 


0 


li 


0 


Figure  11-1.  Schematic  of  Geometrical  Model  of  Erosive  Burning  (Two-Flame  Model). 


feedback  from  these  two  flames  and  the  energy  requirements  for  heating 
the  propellant  from  Its  Initial  temperature  to  the  burning  surface  tempera- 
ture and  decomposing  It  yields  (assuming  that  the  heat  feedback  required 
per  unit  mass  of  propellant  consumed  Is  Independent  of  burning  rate) : 


kl(TAP  ‘ V k2(Tf  ' V 

Vqfeedback  “ Lj  * ^iff  + LKln 


(5) 


T^p  ■ ammonium  perchlorate  flame  temperature 
■ propellant  flame  temperature 


Tg  ■ surface  temperature 


^feedback 


heat  feedback  flux  from  gas  flames  to  propellant  surface 


k ■ ratio  of  heat  feedback  flux  from  AP  flame  to  the 

average  temperature  gradient  from  that  flame  to  the 
surface 


k2  * ratio  of  heat  feedback  flux  from  columnar  diffusion 
flame  to  the  average  temperature  gradient  from  that 
flame  to  the  surface. 


The  situation  pictured  as  prevailing  with  a crossflow  is  shown 

in  the  second  part  of  Figure  II-l.  Since  L_  and  L„.  are  both  kinetically 

I Kin 

controlled  and  are  thus  simply  proportional  to  a characteristic  reaction 

time  (which  is  assumed  to  be  unaffected  by  the  crossflow)  multiplied  by 

the  propellant  gas  velocity  normal  to  the  surface  (which  for  a given 

formulation  is  fixed  by  burning  .rate  and  pressure  alone),  these  distances 

are  fixed  for  a given  formulation  at  a given  burning  rate  and  pressure 

independent  of  the  crossflow  velocity.  Of  course,  since  c.-ossflow  velocity 

affects  burning  rate  at  a given  pressure  through  its  influence  on  the 

diffusion  process  as  discussed  below,  lj  and  L^n  are  influenced  through 

the  change  in  burning  rate,  but  this  is  simply  coupled  into  a model  by 

expressing  L and  L as  explicit  functions  of  burning  rate  and  pressure 
I Kin 

in  that  model.  The  important  point  is  that  they  can  be  expressed  as  functions 


of  these  two  parameters  alone  for  a given  propellant.  However,  the 
distance  of  the  mixing  zone  from  the  propellant  surface  is  directly  affected 
by  the  crossflow.  It  may  be  shown  through  geometrical  arguments  coupled 
with  the  columnar  diffusion  flame  height  analysis  presented  by  Schultz, 
Penner  and  Green  (A),  that  L^^  measured  along  a vector  coincident  with 
the  resultant  crossflow  and  transpiration  velocities  should  be  approximately 
the  same  as  L^^  in  the  absence  of  a crossflow  at  the  same  burning  rate 
and  pressure  (except  at  very  high  ratios  of  local  crossflow  velocity  to 
transpiration  velocity).  That  is,  the  magnitude  of  1^  ^ is  essentially 
independent  of  the  crossflow  velocity  although  its  orientation  is  not. 

Thus,  the  distance  from  the  surface  to  the  "average"  mixed  region  is  de- 
creased to  sin  0 where  0 represents  the  angle  of  the  average  flow 

vector  in  the  mixing  region.  (See  Figure  II— 1. ) The  heat  balance  at 
the  propellant  surface  now  yields: 

ki(TAP-  v . k2»f  - y ... 

q feedback  LI  lDiff  ,ln  9 + ‘‘Kin 

This  picture  has  been  used  as  the  basis  of  development  of  a 
first  generation  flame  bending  model  for  prediction  of  burning  rate  versus 
pressure  curves  at  various  crossflow  velocities,  given  only  a curve  of 
burning  rate  versus  pressure  in  the  absence  of  crossflow.  This  model 
employs  no  empirical  constants  other  than  those  obtained  from  regression 
analysis  of  the  no-crossflow  burning  rate  data.  Thus,  although  it  is  not 
as  powerful  as  a model  which  would  permit  prediction  of  erosive  burning 
phenomena  with  no  burning  rate  data  at  all  (but  only  propellant  composition 
and  ingredient  size  data),  it  is  still  a very  useful  tool,  in  that  it 
permits  prediction  of  erosive  burning  characteristics  given  only  relatively 
ealsly  obtained  strand-bomb  burning  rate  data.  (By  comparison,  the 
Lenoir  and  Robillard  model  employs  two  free  constants  which  are  adjusted 
to  provide  a best  fit  of  erosive  burning  data  for  a given  propellant  and 
since  these  constants  vary  from  propellant  to  propellant,  the  Lenoir  and 
Robillard  model  does  not  permit  a priori  erosive  burning  predictions  for 
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new  propellants  without  some  erosive  burning  data,  whereas  the  first 
generation  model  developed  on  this  program  does  not  require  such  data. 

The  general  approach  followed  in  development  of  this  model 

was : 


1.  Derive  expressions  for  L,  and  L , as  functions  of 

1 Diff*  . Kin 

burning  rate  (or  burning  mass  flux,  m^  ),  pressure,  and 
propellant  properties  and  substitute  these  into  a propellant 
surface  heat  balance. 


2.  Work  the  resulting  equation  into  the  form: 

-|l/2 


r * A P 
o 3 


1 + 


2 2 

1 + A,d  P 
5 P 


tq  ■ burning  rate  in  absence  of  crossflow 


P - pressure 

d * particle  size  of  oxidizer  (diameter) 
P 

A3>  A^ , A,.  - empirical  constants 


(7) 


for  burning  in  the  absence  of  crossflow  and  perform  a regression 
analysis  using  no-crossflow  burning  rale  data  to  obtain  best 
fit  values  for  Aj,  A^,  and  A^.  (d  is  the  average  ammonium 
perchlorate  particle  size.  For  a given  propellant,  the  burning 
rate  data  may  be  just  as  effectively  regressed  on  A^,  A^  and 
A5  dp2,  eliminating  the  necessity  of  actually  defining  an  effect- 
ive average  particle  size.) 


3. 


4. 


From  these  results,  obtain  expressions  for  L^.,  and 

as  functions  of  burning  rate  farm,  ) and  pressure. 


Combine  these  expressions  with  an  analysis  of  the  boundary  layer 

flow  which  gives  the  crossflow  velocity  as  a function  of  distance 

from  the  propellant  surface,  mainstream  velocity,  and  propellant 

burning  rate,  to  permit  calculation  of  the  angle  d,  Lj,  LDiff, 

L,.  , and  m,  for  a given  pressure  and  crossflow  velocity. 

Kin  burn 
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Detailed  equation  development  for  this  model  is  presented  in  Reference 
13.  Initial  testing  of  the  model  was  carried  out  using  a systematic  ero- 
sive burning  data  set  taken  by  Saderholm  (2) . (This  was  the  only  systematic 

erosive  burning  data  set  found  in  the  literature  with  sufficient  zero 

2 

crossflow  data  to  permit  evaluation  of  A^,  A^,  and  A^d^  .)  Results  are 
presented  in  Figure  II-2.  As  may  be  seen,  agreement  between  experiment 
and  theory  is  excellent  over  a wide  range  of  pressure  and  crossflow  velo- 
cities. Further  testing  of  this  model  against  data  obtained  in  the  exper- 
imental part  of  this  program  (described  below)  is  presented  later  in  this 
report . 


3. 


Second  Generation  Model 


As  indicated  above,  the  first  generation  model  does  require 
as  input  burning  rate  versus  pressure  data  with  no  crossflow  (l.e. , strand 
data)  for  each  propellant  for  which  erosive  burning  predictions  are  to  be 
made.  A more  fundamental  model  (with  explicit  calculation  of  the  distances 
of  flames  from  the  surface  rather  than  Inference  of  these  distances  from 
zero  crossflow  data)  of  the  propellant  combustion  process  which  would 
permit  prediction  of  burning  rate  versus  pressure  with  or  without  cross- 
flow,  given  only  propellant  composition  and  ingredient  size  data.  Is  highly 
desirable.  During  the  second  year  of  this  program,  this  author  has  been 
working  on  development  of  such  a model  for  the  rather  limited  case  of  pro- 
pellants containing  unimodal  ammonium  perchlorate  oxidizer  and  no  metal 
additives  (with  plans  to  extend  this  model  later  to  treat  multimodal  oxidizer 
and  metallized  propellant  cases).  The  first  thought  was  to  simply  modify 
the  Beckstead-Derr-Price  (BDP)  model  (32)  for  prediction  of  burning  rates 
of  composite  propellants  as  a function  of  pressure  (in  the  absence  of 
crossflow)  to  allow  for  bending  of  the  diffusion  flame (s)  considered  in 
that  model.  However,  upon  careful  review  of  the  BDP  model,  this  author 
found  sufficient  problems  and  areas  of  disagreement  with  that  model  that 
it  was  decided  to  develop  an  entirely  new  composite  propellant  combustion 
computer  code  (embodying  many  of  the  BDP  concepts,  while  modifying  or 
replacing  others)  with  the  flame-bending  mechanism  described  in  the  pre- 
ceding section  embedded  in  the  mathematical  analysis.  Major  modifications 
to  the  BDP  model  included  are: 
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PRESSURE  (atm) 

Figure  11-2.  Eros iv«  Burning  Modal  Predictions  and  Comparisons  with 
Sadarholm  Data.  Transpiration  Effacts  includad. 


1) 


Variation  in  the  ratio  of  local  oxidizer-surface  intersectional 
area  to  the  binder  surface  area  as  the  propellant  surface  regresses 
past  an  oxidizer  particle  is  considered.  (In  the  BDP  model, 
a geometrical  average  ratio  is  used;  this  involves  an  assumption 
that  a lot  of  very  nonlinear  processes  can  be  linearly  averaged.) 

2)  The  kinetics  of  subsurf  ace/ surf  ace  exothermic  reactions  are 
considered,  with  use  of  rate  expressions  based  upon  the  work 
of  Waesche  and  Wenograd  (33).  (In  the  BDP  model,  subsurface/ 
surface  heat  release  is  Included  with  the  endothermic  ingredient 
vaporization  heats,  with  the  resultant  implicit  assumption  that 
the  amount  of  heat  release  in  these  reactions  per  unit  mass  of 
propellant  is  independent  of  such  parameters,  as  burning  rate.) 

3)  A correction  of  an  inconsistency  in  definition  of  areas  in  the 
BDP  model  is  made. 

4)  The  calculation  of  the  dimensionless  stoichiometric  group  needed 
for  calculation  of  the  diffusion  flame  height  via  the  Burke- 
Schumann  (34)  analysis  is  modified.  (The  group  uBed  in  the 

BDP  model  is  inconsistent  with  that  defined  in  the  original 
work  of  Burke-Schumann. ) 

5)  A two-flame  (fuel-gas/oxldizer-gas  columnar  diffusion  flame  and 
ammonium  perchlorate  monopropellant  flame),  rather  than  a 
three-flame  model,  is  used.  (With  correction  of  the  calculation 
of  the  stoichiometry  dimensionless  group  for  the  Burke-Schumann 
analysis,  it  no  longer  appears  necessary  to  differentiate  between 
the  parts  of  the  diffusion  flame  inside  and  outside  of  an  ammonium 
perchlorate  monopropellant  flame.) 

6)  The  procedure  for  calculation  of  heat  feedback  from  the  diffusion 
flame  and  the  AP  monopropellant  flame  is  modified.  (In  the  BDP 
model,  all  flames  are  considered  to  occur  in  flame  sheets  at 
discrete  distances  from  the  surface:  in  the  current  model,  the 
AP  monopropellant  heat  release  is  treated  as  a flame-sheet 

type  heat  release,  but  the  diffusion  flame  heat  release  is  considered 
to  occur  over  a finite  range  of  distances  from  the  propellant 
surface.) 
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7)  The  distance  associated  with  oxidizer-binder  gas  interdiffusion 
in  the  presence  of  crossflow  is  assumed  to  be  reduced  by  s 
factor,  sin  6,  where  6 is  the  angle  of  the  resultant  of  the 
crossflow  and  transpiration  velocities  relative  to  the  surface, 
as  in  the  first  generation  model. 

A major  assumption  made  in  the  BDP  model  (and  variants  thereof) 
is  that  one  may  work  in  terms  of  an  average  oxidizer-fuel  ratio  for  a given 
size  oxidizer  particle.  In  reality,  however,  an  oxidizer  particle  and  the 
fuel  surrounding  it  (and  associated  with  it)  will  be  receiving  heat  feedback 
from  a diffusion  flame  of  strongly  varying  oxldizer/fuel  ratio  during  its 
burning.  As  the  oxidizer  particle  first  becomes  exposed  to  the  surface, 
with  only  its  tip  showing,  the  local  oxidizer-fuel  ratio  will  be  quite  low. 

As  the  burning  surface  passes  the  equator  of  the  particle,  however,  the 
oxidizer-fuel  ratio  will  be  comparatively  high,  and  as  the  particle  burns 
out,  the  ratio  will  again  be  low.  Implicit  in  the  BDP  use  of  an  "average" 
oxidizer-surface  planar  intersectional  area  is  the  assumption  that  all  of 
the  highly  non-linear  dependencies  of  burning  rate,  flame  temperature,  and 
consequently  heat  feedback  from  the  diffusion  flame  can  be  linearly 
averaged  over  the  range  of  the  variations  during  regression  of  the  propellant 
surface  through  the  oxidizer.  Things  may  work  out  this  way,  but  this 
appears  to  this  author  to  be  a somewhat  risky  a priori  assumption. 
Accordingly,  in  this  model  (limited  thus  far  to  unimodal  oxidizer)  an 
attempt  is  made  to  allow  for  the  variation  in  local  oxidizer/fuel  ratio 
associated  with  the  burning  of  an  individual  oxidizer  particle  due  to  the 
variation  in  relative  oxidizer-fuel  surface  intersectional  areas  as  the 
surface  regresses  through  the  particle. 

In  deciding  how  to  treat  this  variation  (or,  indeed,  whether  to 
treat  it)  one  must  first  address  the  question  of  propellant  surface  and 
subsurface  response  to  variation  in  heat  feedback  flux  from  the  varying 
oxldizer/fuel  gas-phase  diffusion  flame.  If  the  burning  rate  response  is 
very  slow,  such  variations  in  feedback  flux  are  damped  out  and  the  aver- 
aging procedure  of  BDP  is  probably  adequate.  If,  on  the  other  hand, 
response  of  burning  rate  variations  to  heat  feedback  flux  variations  is 
sufficiently  fast,  one  may  use  quasi-steady  state  calculations  of  the  burning 
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rate  at  each  fuel/oxidizer  area  ratio  during  the  regression  of  the  burning 

surface  through  the  particle  and  then  properly  average  these  to  arrive  at 

an  average  burning  rate.  In  between  these  extremes  lies  great  difficulty. 

A transient  heat  conduction  program  allowing  for  surface  ablation  was 

employed  to  examine  the  response  of  ablation  rate  to  variation  in  heat 

6 2 2 

flux  to  the  surface.  Variations  in  heat  flux  up  to  10  cal/cm  sec 
(corresponding  to  approximate  doubling  of  heat  feedback  flux  from  a 
typical  steady-state  value  in  0.50  msec,  the  time  required  for  a propellant 
burning  at  2 cm/sec  to  regress  10  microns)  were  examined.  In  all  cases, 
the  burning  rate  response  was  found  to  track  the  feedback  flux  variation 
within  10  percent.  Accordingly,  it  was  concluded  that  use  of  a quasi- 
steady-state  approach  to  calculation  of  propellant  burning  rate  at  various 
oxidizer/fuel  ratios  associated  with  different  intersections  of  the  pro- 
pellant burning  surface  with  a given  oxidizer  particle  would  not  be  seriously 
in  error. 


As  mentioned  earlier,  this  second  generation  model  is  presently 
limited  to  unimodal  oxidizer  particle  size.  Having  concluded  that  one  can 
use  a quasi-steady-state  approach  to  calculating  burning  rate  as  a function 
of  the  ratio  of  planar  areas  of  oxidizer  and  associated  fuel  intersected 
by  the  regressing  surface,  one  is  next  faced  with  the  question  of  how 
to  calculate  the  distribution  of  these  areas.  Since  composite  propellants 
are  normally  quite  highly  loaded  with  solid  oxidizer  in  the  rubber  fuel 
binder,  and  since  with  unimodal  oxidizer  propellants  the  desire  for  these 
high  loadings  tends  to  lead  to  loadings  approaching  maximum  theoretical 
loading,  it  was  decided  that  as  a reasonable  approximation,  one  might 
assume  a regular  packing  of  oxidizer  crystals  in  the  binder  corresponding 
to  the  arrangement  of  a cubic  closest  packing  array,  though  with  the  spacing 
larger  than  that  for  a true  cubic  closest  packing,  corresponding  to  less 
than  100  percent  of  theoretical  maximum  loading.  Simple  geometrical 
considerations  then  permit  one  to  calculate  the  characteristic  lattice 
dimension  (where  lattice  spacings  in  three  mutually  orthogonal  planes 
are  given  by  D^,  0.866  D^,  and  0.82  D^)  as: 


/0. 737 \1/3  n „ / 
°1  " \ VL0  / °o  \ 


(0.737[WFO/pqx  + (1  - WF0)/pf]p 


ox  I 


,1/3 


WFO 


(8) 
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VLO  = volumetric  fraction  oxidizer  in  propellant 

D = oxidizer  diameter 
o 

WFO  - weight  fraction  oxidizer  in  propellant 

3 

p ■ oxidizer  density  (1.95  gm/cm  for  ammonium  perchlorate) 

OX  M 

p^.  = fuel  (binder)  density  (0.92  gm/cin  for  HTPB  or  CTPB 
binders) 

It  is  arbitrarily  assumed  that  the  propellant  burns  in  the  direction 
in  which  the  planes  of  oxidizer  are  separated  by  0.82  D^.  This  distance 
is  broken  up  into  equally  spaced  increments  and  straightforward  geometrical 
relations  are  then  used  to  calculate  the  planar  intersection  area  of  the 
burning  surface  with  the  oxidizer  (APOX)  and  its  associated  fuel  planar 
area  (AFU)  at  each  of  the  intersection  planes,  with  the  assumption  that 
wherever  two  layers  of  oxidizer  overlap  the  fuel  is  apportioned  between 
them  in  the  ratio  of  their  planar  surface  intersection  areas.  The  result 
of  these  calculations  is  a table  of  planar  oxidizer-surface  intersectional 
area  (APOX)  and  associated  fuel  surface  area  (AFU)  versus  distance  of  the 
intersection  plane  from  the  top  of  the  particle  (XDTOP) . Results  of  a 
typical  calculation  are  presented  in  Table  II-III.  Burning  rates  for  each 
of  these  conditions  (starting  at  the  top  of  the  particle  since  one  must 
allow  for  different  regression  rates  of  fuel  and  oxidizer)  are  then  calcu- 
lated as  described  below  and  an  averaging  procedure,  also  described  below, 
is  then  used  to  calculate  the  propellant  burning  rate. 

Next  let  us  address  the  question  of  calculation  of  propellant 
burning  rate  at  each  of  the  conditions  defined  by  the  various  distances 
of  the  burning  surface  Intersection  plane  from  the  top  of  the  oxidizer 
particle,  as  listed  for  the  example  in  Table  II-III.  First,  since  as  men- 
tioned above,  different  oxidizer  and  fuel  regression  rates  are  to  be  allowed, 
one  must  address  rather  carefully  the  questions  of  surface  geometry  and 
mass  conservation  at  the  surface.  In  this  model,  as  in  the  BDP  model,  the 
fuel  is  assumed  to  regress  in  a planar  manner,  and  the  oxidizer-fuel  surface 
is  forced  to  be  continuous  at  their  intersection.  These  restrictions, 
coupled  with  the  fact  that  the  linear  regression  rates  of  fuel  and  oxidizer 
parallel  to  their  directions  of  regression  are  allowed  to  differ,  force 
the  oxidizer  surface  to  assume  a curved  shape  as  it  regresses.  Oxidizer 
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TABLE  II-III.  Variation  of  Oxidizer  Planar  Surface  Intersection 

Area  and  Associated  Fuel  Surface  Area  with  Distance 
of  Intersection  Plane  from  Top  of  the  Particle  - 
Typical  Case 

Particle  Diameter  * 20  microns 

Weight  Fraction  Oxidizer  ■ 0.73 

3 

Oxidizer  Density  ■ 1.95  gm/crc 

3 

Binder  Density  ■ 0.92  gm/cm 

Volumetric  Oxidizer  Loading  ■ 56.05  percent 

D^  ■ 21.91  microns 


XDTOP  APOX  AFU 


Increment 

Number 

Distance  from  Top 
of  Particle  (microns) 

Planar  Oxidizer  2 

Intersectional  Area  (microns) L 

Planar  Fuel  ^ 
Area  (microns) 

1 

0.119 

7.4 

19.1 

2 

1.017 

60.6 

147.2 

3 

1.915 

108.8 

280.4 

4 

2.813 

151.9 

263.8 

5 

8.712 

18.19 

225.8 

6 

4.610 

222.9 

192.8 

7 

5.508 

250.8 

165.0 

8 

6.407 

273.6 

142.1 

9 

7.305 

291.3 

124.4 

10 

8.203 

304.0 

111.7 

11 

9.101 

311.6 

104.1 

12 

10.0 

314.2 

101.5 

13 

10.899 

311.6 

104.1 

14 

11.797 

304.0 

111.7 

15 

12.695 

291.3 

124.4 

16 

13.593 

273.6 

142.1 

17 

14.492 

250.8 

165.0 

18 

15.390 

222.9 

192.8 

19 

16.288 

189.9 

225.8 

20 

17.187 

151.9 

263.8 

21 

18.085 

108.8 

280.4 

22 

18.983 

60.6 

147.2 

23 

19.881 

7.4 

19.1 
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mass  fluxes  may  be  expressed  relative  to  either  the  actual  curved  surface 
area  or  the  planar  projection  of  this  area,  the  two  values  being  related  by: 

r 

ra  (APOX)-m  (ASOX)  - — (ASOX)  (9,9a) 

OX , P ' OX, B Pox 


tn 

ox,p 

• 

m 

ox,s 


APOX 

ASOX 


oxidizer  mass  flux,  based  on  planar  surface  projection 

oxidizer  mass  flux,  based  on  actual  total  curved  surface 
area 

planar  projection  of  oxidizer  surface  area 
total  curved  oxidizer  surface  area 

linear  regression  rate  of  oxidizer,  normal  to  its  surface 
oxidizer  density 


The  average  mass  flux  of  fuel  and  oxidizer  normal  to  the  mean  regression 
plane  is  given  by: 


m 


m (AFU)  + m (APOX) 
fuel  ox,pv 

AFU  + APOX 


m.  . (AFU)  + m (ASOX) 
fuel  ox, s 

AFU  + APOX 


(10) 


It  is  important  to  know  the  value  of  ASOX  at  each  plane  since  the  Arrenhius 

expression  relating  oxidizer  mass  flux  to  surface  temperature  must  be 

written  in  terms  of  m to  be  meaningful: 

ox,s 


m - B ex  (-E  /RT  ) 

OX, 8 OX  OX  S 


(ID 


Bqx  - pre-exponential  rate  factor 

E 11  activation  energy  of  oxidizer  surface  ablation  reaction 
ox 

Ts  ■ surface  temperature 


A similar  expression  for  the  fuel  pyrolysis  rate: 
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(12) 


m,  , ■ B-  , exp(-E,  ,/RT  ) 
fuel  fuel  r fuel  s' 


enables  one  to  calculate  the  ratio  of  oxidizer  and  fuel  regression  rates 
as  a function  of  surface  temperature: 


ox 


' fuel 


m P,  , p,  , B exp(-E  /RT  ) 
ox.s  fuel  fuel  ox  ox  s 

.1 — — — m ■ ■ 


m 


fuel  Pox 


p B,  .exp(-E,  ,/RT 
Kox  fuel  fuel  s 


(13) 


There  is  considerable  uncertainty  as  to  best  values  to  be  used  for  B , 

ox 

Bfud»  Eqx*  and  Efuei:  thus,  parametric  study  of  the  effects  of  these 
values  is  required.  Note  that  it  has  been  assumed  here  that  the  oxidizer 
and  fuel  surface  temperatures  are  equal.  This  is  probably  not  a particularly 
good  assumption,  but  relaxing  it  requires  a rather  complex  three-dimensional 
heat  transfer  analysis. 

Now,  how  does  one  go  about  calculating  ASOX  for  succeeding  regression 
intervals  through  the  oxidizer  particle?  First,  it  is  assumed  (approximated) 
that  at  the  first  increment  after  the  tip  of  the  particle  becomes  exposed 
(in  Table  II-III,  when  the  distance  from  the  top  of  the  particle  is  0.119 
microns)  the  oxidizer  surface  is  planar.  The  procedure  outlined  below  for 
calculation  of  burning  rate,  given  the  local  oxidizer/fuel  area  ratio,  is 
then  used  to  calculate  the  oxidizer  and  fuel  linear  regression  rates  under 
the  conditions  given  for  this  first  increment.  The  fuel  regression  rate 
is  then  used  to  calculate  the  time  for  the  regressing  fuel  to  reach  the 
second  Increment  (distance  from  the  initial  particle  top  of  1.017  microns 
in  Table  II-III)  as: 


TAU2-(XDTOP2  - XDT0P1)/rf  (14) 

The  distance  which  the  center  of  the  AP  particle  peak  regresses  in  that 
time  is  then  calculated  as: 


A(DELOX)  - rQX(TAU2) 


(15) 
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Similar  procedures  are  followed  for  each  succeeding  Increment,  yielding 
for  each  XDTOP  (distance  of  fuel  surface  from  the  initial  top  of  the 
oxidizer  particle)  a value  of  DELOX  (distance  of  the  center  of  the  oxidizer 
surface  from  the  initial  top).  The  geometrical  method  depicted  in  Figure 
11—3  (for  a case  where  the  oxidizer  regresses  more  slowly  than  the  fuel) 
is  then  used  to  calculate  ASOX  at  each  calculational  increment.  Applying 
the  Pythagorean  Theorem  to  the  larger  tiangle  and  using  |h|  ■ | XDTOP  - DELOX |, 

2RX  | XDTOP  - DELOX | - A2  + (XDTOP  - DELOX)2  (16) 

while  similar  analysis  of  the  smaller  triangle  yields: 

A2  - D (XDTOP)  - XDTOP2  (17) 

o 

Elimination  of  A from  these  equations  and  use  of  ASOX  ■ 2ttR^|X  TOP  - DELOX| 
then  gives: 

ASOX  - tt[(XDTOP)D  - 2 (XDTOP)  (DELOX)  + DELOX2]  (18) 

o 

For  the  calculation  of  the  columnar  diffusion  flame  between  fuel 
gases  from  the  binder  pyrolysis  and  oxidizer  gases  from  the  oxidizer  decom- 
position (discussed  later)  several  other  parameters  associated  with  the 
surface  configuration  of  the  oxidizer  particle-associated  fuel  combination 
at  each  Increment  must  be  calculated.  First,  the  combined  radius  of  the 
oxidizer  and  binder  gas  streams  (in  this  model,  a modified  Burkp-Schumann 
analysis  with  a fuel  annulus  surrounding  an  oxidizer  gas  core  is  employed 
for  the  columnar  diffusion  flame  calculation)  is  calculated  as: 

Rgg  - /(AFU  + APOX)  h (19) 

RgS  ■ outer  radius  of  fuel  annulus 

In  line  with  the  requirement  in  the  Burke-Schumann  analysis  that  the  linear 
velocity  of  the  fuel  and  oxidizer  streams  have  the  same  Initial 
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Figure  11*3.  Schematic  Demonitrating  Calculation  of  Oxidizer  Surface  Area  at  Some 
Time,  TAU,  After  First  Exposure  of  the  Top  of  the  Oxidizer  Particle. 
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value,  the  oxidizer  and  fuel  gas  streams  leaving  the  surface  are  assumed 
to  adjust  their  areas  quickly  from  the  planar  areas  of  the  solids  to  meet 
this  requirement.  Since  the  temperature  at  the  surface  Is  assumed  to  be 
the  same  for  the  fuel  and  oxidizer,  and  pressures  are  equal,  this  leads  to 
an  expression  for  the  radius  of  the  Inner  oxidizer  gas  jet  of: 


m,  , (AFU) (MW) 
fuel  ox 

m (ASOX) (MW) , . 

ox,s  fuel 


(20) 


Lnc  ■ oxidizer  gas  column  radius 

Ob 

MW  - molecular  weight  of  oxidizer  gases  (estimated  to  be 
ox  approximately  35  for  ammonium  perchlorate) 

MW  = molecular  weight  of  fuel  pyrolysis  gases  (estimated  to 
ue  be  approximately  19  for  CTPB  or  HTPB  binder) 


The  linear  gas  velocity  away  from  the  surface,  also  required  in 
the  modified  Burke-Schumann  analysis,  as  well  as  for  calculation  of  charac- 
teristic reaction  distances  (products  of  reaction  times  and  this  velocity) 
is  calculated  as: 


V 


gas, surf 


"o,.B(AS0X)RI.urf 

<MWox’(LBS)2p 


(21) 


V f ■ gas  velocity  away  from  propellant  our  face 

y sur t 

P » pressure 


Finally,  the  ratio  of  the  molar  fuel/oxidizer  ratio  to  stoichio- 
metric molar  fuel/oxidizer  ratio  for  the  combined  fuel  and  oxidizer  streams 
(<f>),  also  required  in  the  modified  Burke-Schumann  analysis,  is  calculated  as: 


♦ 


"bs2  V 
Ss2 


(22) 


Cj  ■ initial  concentration  of  fuel  in  the  binder  pyrolysis 
product  gases 

C ■ initial  net  (excess  of  oxidizer  over  fuel  molecules)  con- 
centration of  oxidizer  in  the  oxidizer  deompostion  product 
gases 

1 « stoichiometric  moles  of  fuel  per  mole  of  oxidizer. 
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Since  the  initial  pressure  and  temperature  of  the  fuel  and  gas  streams  are 
the  same,  the  concentration  ratio  Cj/C^  may  be  replaced  by  a ratio  of  mole 
fractions  ^/Y^.  The  mole  fraction  of  fuel  in  the  binder  pyrolysis  products 
for  HTPB  and  CTPB  binder  systems  in  the  absence  of  subsurface  reactions 
depleting  some  of  this  fuel  (as  discussed  later)  is  approximately  0.96, 
while  the  initial  net  mole  fraction  of  oxidizer  molecules  in  the  decompos- 
ition products  of  ammonium  perchlorate  is  approximately  0.28.  For  an 
HTPB/AP  or  CTPB/AP  system,  the  stoichiometric  moles  of  fuel  per  mole  of 
oxidizer  are  about  0.6,  calculated  on  the  basis  of  1^0,  CO^  stoichiometry. 
(These  calculations,  as  well  as  the  calculations  of  molecular  weights  for 
ammonium  perchlorate  and  HTPB  decomposition  gas  streams  mentioned  earlier 
are  based  on  HTPB  binder  yielding  approximately  0.16  moles  CO,  0.03  moles 
N-,  and  5.0  moles  C H per  100  grams;  and  AP  yielding  approximately  5 

^ 1*4  4 • 1 

moles  H2O,  2 moles  NH^,  2 moles  C10^,  1 mole  HC10^,  1 mole  HC1,  1 mole  ^0, 
1.5  moles  0^,  0.5  moles  Cl^,  1 mole  NO,  and  1 mole  per  585  grams.) 

With  subsurface  reactions,  the  fuel  mole  fraction  is  reduced  by  a factor 
(1  - 6)  while  the  oxidizer  mole  fraction  is  reduced  by  a factor  (1  - a) 
where  B is  the  fraction  of  fuel  consumed  in  subsurface  reactions  and 
a is  the  fraction  of  oxidizer  consumed  in  these  reactions.  (See  later 
discussion.)  With  these  substitutions,  Equation  22  becomes: 


4> 


V2.o(1  - 8)(RBS2  - V> 
1 Yl,o(1  - “>  LBS2 


(23) 


mole  fraction  of  fuel  in  binder-pyrolysis  products  in 
absence  of  subsurface  reactions 

mole  fraction  of  oxidizer  in  oxidizer  decomposition 
products  in  absence  of  subsurface  reactions. 


The  burning  rate  of  the  propellant  at  any  given  set  of  oxidizer/ 
fuel  conditions  (any  regression  increment)  is  controlled  by  heat  releases 
(exothermic  reactions)  at  various  locations.  In  this  model,  we  consider 
three  principal  heat  release  zones:  (1)  heat  release  in  a thin  subsurface 
zone  quite  near  (and  including)  the  propellant  surface;  (2)  heat  release 
in  the  gas-phase  above  the  propellant  from  ammonium  perchlorate  decomposition 


27 


products  burning  as  a monopropellant;  and  (3)  heat  release  from  a diffusion 
flame  between  AP  decomposition  (and  monopropellant  flame)  products  and  fuel 
vapor  released  by  binder  pyrolysis. 

The  subsurface/surface  heat  release  is  calculated  by  an  iterative 
process,  coupled  with  the  remainder  of  the  model,  In  which  an  estimate  of 
the  subsurface  temperature  profile  is  made  and  substituted  into  an  Arrenhlus 
rate  expression  representing  subsurface  heat  release  rate  data  measured  by 
Waesche  and  Wenograd  (33),  which  is  then  integrated  from  the  surface  into 
the  propellant  to  obtain  the  total  subsurface  heat  release  per  unit  mass  of 
propellant.  This  procedure  differs  quite  a bit  from  that  of  the  BDP  model, 
in  which  the  amount  of  subsurface  heat  release  per  unit  mass  of  propellant 
is  assumed  to  be  a constant,  independent  of  such  parameters  as  burning  rate, 
and  is  included  with  the  binder  heat  of  vaporization.  Since  the  subsurface 
temperature  profile  steepens  rapidly  with  increasing  burning  rate,  our 
procedure  results  in  the  subsurface  heat  release  per  unit  mass  of  propellant 
decreasing  with  Increasing  rate.  As  will  be  discussed  later,  the  surface 
energy  balance  in  this  model  is  written  with  the  surface  area  of  the  oxidizer 
and  associated  fuel  as  the  basis:  thus  all  terms  appear  in  the  units  of 
energy/time.  For  bookkeeping  convenience,  the  surface/subsurface  heat 
release  term  is  written  as: 


SMRBO  °*»* 


m (ASOX) 


(25) 


”fuel(AFU) 


SMRBO  ■ stoichiometric  ratio  (mass)  of  binder  to  oxidizer 
(0.111  for  HTPB/AP  or  OTPB/AP  systems) 


mc  , - mass  flux  of  fuel 
fuel 


AFU 


fuel  surface  area 


Based  upon  Waesche  and  Wenogiad  subsurface  reaction  rate  data,  the  fraction 
of  oxidizer  reacted  per  unit  time  is  given  for  AP/CTPB  systems  as  a function 
of  temperature  by: 


sub 


sub 


R ■ B . exp(-E  . / RT) 
o Sub  K sub 

5.1(10)12/D  for  D > 20p 
o o 

2.55(10U)  for  D < 20y 

o 


(26) 


Dq  ■ oxidizer  particle  diameter  in  microns 


Bsub  " Pre“exPonential  for  subsurface  reaction  rate  equation 
Esub  " activation  energy  for  subsurface  reaction  rate  equation 
T ■ temperature 


The  unperturbed  (uncoupled)  subsurface  temperature  profile  is  given  by: 


T - (T  - T )exp(r  p C x/X  ) + T 
s o ox  ox  pox  ox  o 


(27) 


surface  temperature 


Tq  * propellant  bulk  temperature 

r - linear  regression  rate  of  oxidizer 
ox 

■ oxidizer  density  (1.95  gm/cm  for  ammonium  perchlorate) 

C “ oxidizer  heat  capacity  (approximately  0.4  cal/gm  °K  for  AP) 
pox 


ox 


oxidizer  thermal  conductivity  (approximately  0.001  cal/cm  sec 
*K  for  AP) 

x ■ distance  from  surface  (sign  convention  such  that  it  is  nega- 
tive below  the  surface). 
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Substitution  of  Equation  27  into  Equation  26  and  Integration  of: 


/ 


X* 

Rdx 

a 


r 

ox 


(28) 


x'  s distance  below  surface  (sign  convention  such  that  it  is 
negative)  at  which  reaction  rate  drops  to  ^1  percent 
of  its  surface  value 


yields,  with  use  of  the  approximation  e U % 1 - u,  in  Equation  27  and 

x’  *-0.25  X /r  p C : 

ox  ox  ox  pox 


B X e DUU  "RT  ‘[1 
^ sub  ox  8 

a 'v  2 


/RT  2 

sub'  sOT  2.  -0.25E  . (T  - T )/RT  . 

RT  [1  - e sub  s os] 


r p C E . (T  - T ) 
ox  ox  pox  sub  s o 


(29) 


Thus  Equation  29  relates  the  mass  fraction  of  oxidizer  reacting  exothermically 
at  or  below  the  surface  to  the  surface  temperature  and  linear  regression 
rate  of  the  oxidizer  and  may  be  used  in  an  iterative  procedure  with  the 
surface  energy  balance  given  later  and  Equation  11  to  solve  for  all  three 
quantities. 

As  regards  gas-phase  heat  release  zones,  a two-flame  approach 
was  chosen  for  this  model,  the  two  flames  being  an  AP  monopropellant  flame 
and  a columnar  diffusion  (Burke-Schumann)  flame.  The  reasons  that  a two- 
flame  rather  than  a three-flame  model  (as  the  BDP)  was  chosen  were: 


1)  Mathematical  simplication. 

2)  Lack  of  apparent  difference  in  a diffusion  flame  between  AP 
decomposition  products  and  fuel  and  a flame  between  AP  monopro- 
pellant flame  products  and  fuel.  In  both  cases,  the  overall 
stoichiometry  is  the  same  since,  while  AP  decomposition  products 
bring  more  oxidizer  into  a binder  fuel  stream  then  do  AP  monopro- 
pellant flame  products,  they  also  bring  more  fuel,  with  the  result 
that  the  overall  mixture  ratio  at  a given  point  is  nearly  the 
same. 


30 


3)  Provisions  were  made  in  calculation  of  the  AP  monopropellant 

heat  release  to  allow  for  consumption  of  reactants  for  that  flame 
in  that  part  of  the  columnar  diffusion  flame  which  occurred 
inside  the  AP  flame. 


Three  distance  parameters  are  important  in  calculating  heat  feed- 
back from  these  gas  flames  to  the  propellant  surface.  These  are  pictured 
in  Figure  II-4.  These  distances  are  FH90  sin  0,  L^p  and  L . FH90  refers 
to  the  distance  associated  with  completion  of  90  percent  of  the  mixing  of 
fuel  and  oxidizer  gas  products.  (If  there  were  no  reaction  delay,  this 
would  be  equivalent  to  the  90  percent  heat  release  point.)  refers  to 

the  reaction  distance  associated  with  the  binder  gas-oxidizer  gas  flame, 
and  L^p  refers  to  the  reaction  distance  associated  with  the  monopropellant 
AP  product  flame  (both  being  characteristic  times  divided  by  the  gas  velo- 
city away  from  the  surface  as  in  the  first  generation  model.)  As  before, 
flame  bending  is  assumed  to  reduce  the  distance  td  the  end  of  the  columnar 
diffusion  heat  release  (90  percent  point)  by  reducing  FH90  to  FH90  sin  0, 
measured  perpendicular  to  the  surface.  FH90  is  calculated  as  a function  of 
various  parameters  using  a modified  Burke-Schumann  analysis  as  described 
below.  A series  of  calculations  of  FH90  as  a function  of  these  parameters 
were  generated  externally  and  correlations  of  the  results  were  used  in  the 
final  program.  In  this  model,  it  is  assumed  that  the  fraction  of  planar 
projection  of  surface,  AP0X/(AP0X  + AFU)  receives  flux  from  both  the  AP 
and  columnar  diffusion  flames  (the  latter  at  the  adiabatic  flame  tempera- 
ture, T*,  which  is  a function  of  the  oxidizer/fuel  ratio,  m [ASOX]/ 

1 OX  f 8 

m.fueifAFU])  while  the  remaining  fraction  of  the  surface  receives  flux  only 
from  this  diffusion  flame;  however,  these  fluxes  are  assumed  to  smear  out 
uniformly  in  the  propellant.  Thus,  the  total  heat  flaw  from  the  gas-phase 
heat  release  zones  is  given  as: 


APOX 


gas 


^series  flames^  + 


AFU 


^diffusion  flame^ 


& 

i 


4 
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Figure  li-4.  Schematic  Showing  Key  Dimensions  Relating  to  Gas-Phase  Heat  Release. 


Heat  release  from  the  AP  monopropellant  flame  la  assumed  to  be  a planar  one, 
resulting  In  a discontinuity  in  the  temperature  derivative  at  its  point  of 
release,  while  the  columnar  diffusion  flame  is  assumed  to  release  its 
heat  uniformly  between  x ■ L and  x ■ Lov  + FH90  sin  0.  Actually,  it  would 
be  better  (more  physically  realistic)  to  assume  different  uniform  heat 
releases  in  the  L_v<  x < L._  and  the  L < x < L„  + FH90  regions  with  the 
difference  related  to  the  heat  release  at  the  AP  reaction  plane.  Sample 
calculations,  however,  indicate  that  the  effect  on  predicted  burning  rates 
is  small,  while  the  added  complexity  of  this  approach  is  considerable. 

With  use  of  multiple  algebraic  manipulations,  we  arrive  at  the  following 
expressions  for  the  heat  fluxes  at  the  surface  (allowing  for  reduction 
of  reactants  available  for  the  AP  monopropellant  reaction  by  occurrence  of 
some  diffusion  flame  reaction  closer  to  the  surface): 

6 I1  ft  CPFH90  sin9  (e  L)f 

Z1  1 z2 

6 “ A C FH90  sin  $(e  * 15 

P 


^diffusion  flame 


i C (T  - T 

— ! k1 — 1 

(e  1 - 1) 


G ■ q 

^series  flames  diffusion  flame  + 


mQAp(l  - a) 


i (lap  ~ lrx^  m 
_ " FH90  sin  8J  L~ 


mC  FH90  sinB  z. 

2 „ 3 


- e 


m 9 

m 


C FH90  sin  9 z.  z „ 

p 12 

e - e 


+ 1 


(32) 


z,  - mC  (Lov  + FH90  sin  0)/X 

1 p KA 

z0  • mC  FH90  sin  0/X 

2 P 

z3  - mCp(FH90  sin  0 + - LAp  )/X 

m ■ average  surface  mass  flux  (based  on  planar  area) 

C ■ gas  heat  capacity  (function  of  oxldizer/fuel  ratio  as 
^ discussed  later) 

- flame  temperature 

T - surface  temperature 
8 

X ■ thermal  conductivity  of  gas  (discussed  later) 

Q * Heat  release  per  unit  mass  associated  with  HC10,(g)  + 
NHj(g)  -*•  Equilibrium  Products 

a - Fraction  of  AP  reacted  below  surface 
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Next,  let  us  consider  the  calculation  of  the  distances  L , 

Ka 

L^p,  and  FH90  (and  FH90  sin  9).  The  distance  FH90,  which  is  calculated  from 
a modified  Burke -Schumann  columnar  diffusion  flame  analysis  (modified  to 
allow  for  axial  diffusion)  is  defined  as  the  distance  from  the  starting  plane 
at  which  90  percent  of  the  fuel  (for  oxidizer-rich  cases)  or  90  percent  of 
the  oxidizer  (for  fuel-rich  cases)  will  be  consumed,  assuming  infinite 
reaction  kinetics.  This  definition  of  the  characteristic  diffusion  distance 
differs  from  that  of  the  BDP  model  where  the  characteristic  distance  is 
defined  as  the  distance  from  the  starting  plane  to  the  point  of  closure 
of  the  flame  over  the  oxidizer  (fuel-rich  cases)  or  the  fuel  (oxidizer- 
rich  cases).  One  serious  problem  with  use  of  the  flame  closure  point  to 
define  the  characteristic  distance  is  that  is  has  a singularity  for  stoichio- 
metric situations:  that  is,  for  stoichiometric  oxidizer-fuel  ratio,  the 
flame  does  not  close  and  this  characteristic  distance  goes  to  infinity. 

Since  most  of  the  heat  is  still  released  fairly  close  to  the  surface,  this 
latter  definition  of  a characteristic  diffusion  distance  leads  to  seriously 
misleading  results  as  regards  heat  feedback  to  the  propellant  surface  at 
near  stoichiometric  oxidizer-fuel  ratios:  at  oxidizer-fuel  ratios  far  from 
stoichiometric,  FH90  and  the  distances  associated  with  flame  closure  differ 
only  slightly. 

While  the  Burke-Schumann  analysis  (and  its  modification  allowing 
for  axial  diffusion)  deals  with  a central  fuel  circular  jet  surrounded  by 
an  annular  jet  of  oxidizer,  we  have  here  reversed  the  situation.  This  does 
not  basically  change  the  analysis,  but  could  lead  to  confusion  regarding 
nomenclature  for  anyone  comparing  equations.  The  modified  Burke-Schumann 
analysis  results  in  a series  expression  for  the  concentration  of  oxidizer 
as  a function  of  axial  position  (distance  from  the  propellant  surface) 
and  radial  position  (distance  form  the  center  of  the  oxidizer  particle 
surface: 

*LBS  \ LBS  ^BS 

2 21  2 2 2+ 

hs  -hs^BS  hs  -Ss 


f 
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where : 


■ (actual  molar  fuel/oxldizer  ratio) /(stoichiometric  molar 
fuel/oxidizer  ratio)  (It  appears  that  in  the  BDP  analysis 
the  stoichiometric  fuel/oxidizer  ratio  is  abritrarily  set 
equal  to  1-0  . ) 

L * oxidizer  jet  radius 

Rgg  * outer  radius  of  fuel  annular  jet 

C *>  initial  concentration  of  oxidizer  in  the  central  oxidizing 
gas  stream 

C •»  oxidizer  concentration  at  any  radial  position  r and  axial 
position  x (Negative  values  exist.  Physically  these  repre- 
sent fuel  regions,  with  the  magnitude  of  C representing  the 
concentration  of  oxidizer  needed  to  react  stolchiometrically 
with  the  fuel  at  that  location.) 

r * radial  distance  from  center  of  oxidizer  surface 

x * axial  distance  from  propellant  surface 

V * axial  gas  velocity 
gas 

D ■ gas  diffusivity  (Discussed  later.) 

■ nth  positive  root  of  J^p^R^,,)  * ® 

J ■ zero  order  Bessel  Function 
o 

* first  order  Bessel  Function 

This  equation  was  programmed  for  solution  for  oxidizer  concentration  as 

a function  of  r and  x for  any  given  set  of  four  independent  parameters: 

L , D/V  , <f>,  and  R /L__i  which  may  be  seen  to  specify  the  problem.  In 
BS  gas  BS  BS 

this  program,  the  first  nine  terms  of  the  series  expression  (n  ■ 1...9) 
were  used.  The  program  also  included  integration  of  the  profiles  from 
r ■ 0 to  r * at  various  values  of  x for  evaluation  of  the  fraction 
of  the  deficient  ingredient  reacted  as  a function  of  x.  (Negative  values 
of  C were  set  equal  to  zero,  of  course,  in  this  integration.)  These  fractions 
were  then  plotted  against  x to  determine  FH9Q  as  a function  of  the  four 
independent  parameters: 

FH90  - D/Vgas,  ♦.  VV  (34) 
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An  extensive  set  of  calculations  covering  a wide  range  of  each  of  these 
variables  was  carried  out  and  tabulations  and  correlations  of  the  results 
were  built  into  the  final  combustion  model. 

The  reaction  distances,  L and  L , are  calculated  as  the  products 

KA  Ar 

of  reaction  times  and  gas  velcoity  away  from  the  surface.  Using  the  Zeldo- 
vich  approach  for  premixed  flame  analysis  along  with  several  minor  approxi- 
mations which  will  not  be  detailed  here,  we  arrive  at  for  the  oxidizer-fuel 

reaction  distance,  LDV: 

KX 


2„  2, 


K0F  Vgas,aurf(1  + ^ Tf  (exptEACT,0F/RTF] ) 


RX 


P* 


(35) 


E ■ activation  energy  for  the  fuel-oxidizer  reaction, 

’ probably  about  11,000  calories/mole  for  ammonium 
perchlorate  systems  based  on  the  data  of  Powling 

Kqf  ” constant  including  the  pre-exponential  for  the 

reaction  rate  term  along  with  several  other  propor- 
tionality constants  such  as  the  gas  law  constant 


For  the  ammonium  perchlorate  reaction  distance,  we  find,  neglecting 
variation  of  temperature  at  the  AP  heat  release  site  (probably  not  a very 
bad  approximation  due  to  the  low  activation  energy  associated  with  the 
ammonia-perchloric  acid  reaction): 

KAT1  V 

L , ftp.  &a8.i.8u.r.f.  . (26 

LAP  p 

K^p  * constant  including  the  pre-exponential  for  the  reaction 

rate  term  along  with  several  other  proportionality  constants 


(Note  that  both  reactions  have  been  assumed  to  be  second  order,  as  shown 
by  the  dependence  of  the  L's  on  pressure.  This  assumption  could  easily 
be  relaxed  to  treat  an  nth  order  reaction  by  replacing  P in  the  denominator 

by  p""1.) 

The  same  approach  to  calculation  of  sin  6 was  used  in  this  model 

(12  13) 

as  in  the  first  generation  model,  ’ the  resulting  equations  being: 


U - v — 

gas,x  - FH90  sin  6 gas, surf  T 


(37) 
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y+  . 1°-5(TO0  sin  e)(Ucro„flov)0-9T(°-V'9 

D0;1  ,((T  + T.l/2)1'8 

channel  s f 


(38) 


U+  - Y+  for  Y+  <5 

U+  - -3.05  + 5.00  lnY+  for  5 < Y+  <_  30 
U+  - 5.5  + 2.5  lnY+  for  Y+  >30 


(39) 


U 


crossflow, x ■ FH90  sin  0 


0.023(0  )0. 9 T 0.18  1 

crossflow  f 


d0:1  , p0-1 

channel 


exp(  60  Ugas>x  m FH9Q  sin  O^crossflow^ 


(40) 


sin  6 


gas,x  - FH90  sin  8 


J U2 

» gas. 


(41) 


+ u 

X ■ FH90  sin  8 crossflow, x ■ FH90  sin  0 


U Q_  - linear  gas  flow  rate  away  from  propellant 

® ’ at  distance  FH90  sin  0 from  the  surface 


U 


crossflow 


.x 


FH90  sin  0 


linear  gas  crossflow  rate  at  distance 
FH90  sin  8 from  the  surface 


0 * angle  between  resultant  velocity  vector 
and  planar  surface  (See  Figure  II-4.) 

Y+  • dimensionless  value  of  FH90  sin  0 

U+  ■ dimensionless  crossflow  velocity  at  x ■ 
FH90  sin  0 

^crossflow  ■ mainstream  crossflow  velocity 

D . , ■ flow  port  hydraulic  diameter 

channel 


As  mentioned  earlier,  the  final  flame  temperature,  T^,,  depends 
on  the  relative  flow  rates  of  fuel  and  oxidizer  gases  at  each  calculational 
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increment  during  progression  of  the  propellant  surface  through  the 
oxidizer.  Accordingly,  a table  of  flame  temperature  calculated  with  a 
thermochemical  equilibrium  program  versus  a parameter  representing  the 
relative  flow  rates  is  generated  for  the  propellant  system  of  interest 
(e.g.,  HTPB-AP)  and  included  as  a tabular  look-up  in  the  final  program 
in  the  form: 


In  addition,  the  product  gas  heat  capacity  is  somewhat  dependent 
upon  this  parameter  and  an  additional  tabular  look-up,  based  on  thermo- 
chemical calculations,  of  the  form: 


C - 

P 


(43) 


is  included. 

At  this  point,  we  have  26  equations  (9,  9a,  10  - 13,  19  - 21, 

23  - 25,  29  - 32,  34-43)  in  27  unknowns  m^,  rox,  m,  mfuel»  T8, 


rfuel*  RBS*  LBS’  Vgas,  surf*  Cp’  ^sub’  ^gas’  ^series  flames, 

x,  FH90,  6,  Tf,  Ug#gjX  „ pjjgQ  8in  e»  Y * u * 

).  For  closure  of  the  problem,  we  finally  write 


^diffusion  flame’  LAP*  LRX’  FH9°’  6*  Tf’  Ugas,x  - FH90  sin  0’  Y » U ’ 
^crossflow, x ■ FH90  sin  0 

an  energy  balance  at  the  propellant  surface  as: 


m.  (AFU)  [C  , .(T  - T ) + Q 1fc.  . 

fuel  p.fuel  8 o melt,  fuel 


fuel 


L UCi  VOp 


m (ASOX) [C 
ox , s P , 


(T  - T ) + 0 ] + m (ASOX)Q  . . (1 

ox  s o Tnelt.ox  ox,s  subl 


- a)  ■ 


(44) 


Cp,fuel 


C 

p,ox 


• e 

q , + q 
sub  ngas 


solid  fuel  heat  capacity  (0.3  cal/gm  °K  for  HTPB  or 
CTPB) 

solid  oxidiaer  heat  capacity  (0.4  cal/gm  °K  for  AP) 


3B 


Q . ■ heat  of  melting  of  binder  (0  for  HTPB  and  CTPB  which 

e ’ are  assumed  not  to  melt) 


Q 


subl 


heat  of  melting  of  oxidizer  (30  cal/gm  for  AP) 


heat  of  pyrolysis  of  fuel  (433  cal/gm,  endothermic, 
for  HTPB) 

heat  of  sublimation  of  oxidizer  (450  cal/gm,  endother- 
mic, for  AP) 


The  resulting  27  equations  (some  of  which,  as  mentioned,  are 

tabular  look-ups  or  correlations  of  one  parameter  as  a function  of  others) 

are  solved  simultaneously  in  a computer  program  for  each  given  set  of  value 

of  ASOX  and  AFU  associated  with  each  surface  regression  Increment.  Among 

the  outputs  of  each  solution  are  values  for  r and  r,  , which  are  used  in 

ox  fuel 

calculation  of  ASOX  from  the  known  APOX  for  the  succeeding  Increment  via 
Equations  14  - 18. 

In  the  solution  procedure,  the  thermal  conductivity  of  the  gas 
and  the  ratio  of  diffusivity  to  gas  velocity  (one  of  the  independent  para- 
meters in  Equation  34),  both  proportional  to  the  square  root  of  temperature, 
are  evaluated  at  the  average  of  the  flame  and  surface  temperature. 

As  mentioned  earlier,  definitive  values  for  B , E , B,  , and 

’ ox  ox  fuel 

Efuei  are  hard  to  come  by,  as  are  values  for  the  gas-phase  reaction  pre- 
exponentials, K^p  and  K^.  In  addition,  further  study  of  the  Waesche  and 
Uenograd  data  (more  of  which  has  been  recently  received  by  this  author) 
may  result  in  changes  in  the  pre-exponential  and/or  activlatlon  energy  used 
for  the  subsurface  reaction  term.  Further,  the  thermal  conductivity  and 
diffusivity  probably  cannot  be  estimated  to  within  better  than  a factor 
of  2 or  3.  Thus,  considerable  testing  of  the  model  against  data  to  deter- 
mine the  best  values  of  these  parameters  (within  reasonable  ranges)  is 
required.  Such  testing  has  not  yet  been  carried  out. 

As  the  program  is  stepped  through  the  succeeding  increments  of 
fuel  plane  distance  from  the  initial  top  of  the  oxidizer  particle  (see 
Table  II-III),  the  oxidizer  surface  will  either  assume  a protruding  bulge 
or  a depression  relative  to  the  planar  fuel  around  it,  depending  upon 
whether  the  oxidizer  linear  regression  rate  is  slower  or  faster  than  the 
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binder  regression  rate.  This  raises  interesting  questions  regarding  the 
"end-game"  if  the  particle  burns  out  before  the  surrounding  fuel.  (Geo- 
metrical considerations  show  that  the  inverse  problem  cannot  occur  as  long 
as  increment  sizes  are  kept  sufficiently  small.)  In  this  case,  there  is 
no  oxidizer  to  burn  with  the  surrounding  fuel  in  succeeding  increments  and 
the  burning  rate  is  set  equal  to  zero  for  these  remaining  increments. 

Three  different  approaches  have  been  taken  to  calculating  the  average 
propellant  burning  rate  from  the  information  obtained  during  the  proced- 
ure of  stepping  through  the  increments  of  regression  of  the  fuel  planar 
surface  past  the  oxidizer  particle  (XDTOP)  increments).  In  the  first  of 
these,  the  burning  rate  is  calculated  by  statistically  averaging  all  of 
the  oxidizer  mass  fluxes  and  fuel  mass  fluxes  over  the  increments  and 
dividing  by  the  propellant  density: 


r 

avg 


;<mo*,p,jApoii)  + "w,j"v 

“propellant  + "V 


(45) 


while  in  the  second  approach  it  is  calculated  by  statistically  averaging 
all  of  the  oxidizer  mass  fluxes  and  then  dividing  by  the  overall  oxidizer 
mass  fraction  and  the  propellant  density: 


I m . APOX, 

i ox,p,J  j 


avg  p ..a  ..  KAPOX,  + AFU.) 

propellant  overall  ^ j j 


(46) 


aoverall 


overall  oxidizer-mass  fraction 


The  fact  that  these  two  procedures  do  not  always  give  the  same  result  is 
tied  in  with  the  "end-game"  problem  mentioned  above.  If  the  oxidizer  burns 
out  before  the  fuel  plane  reaches  the  bottom  of  the  oxidizer,  mass  fluxes  for 
succeeding  increments  are  set  equal  to  zero  in  the  procedure  currently  used. 
Not  only  does  this  result  in  different  answers  by  the  two  above  procedures, 
but  it  also  pulls  the  average  burn  rates  down.  One's  first  temptation  is 
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to  simply  perform  Che  sunning  procedure  over  jusc  chese  IncremenCs  for  which 
a burning  race  is  calculated,  boCh  in  Che  numeraCor  and  denominaCor  of 
Equacions  45  and  46,  buC  lc  is  noC  clear  vhecher  or  noc  Chis  is  more 
physically  realiscic  Chan  Che  procedure  of  summing  over  all  IncremenCs, 
wich  burning  race  sec  equal  Co  zero  for  IncremenCs  in  which  chere  is  no 
oxidizer  for  Che  fuel.  A chird  procedure  of  calculacing  average  burning 
race  was  developed  which  basically  does  Cake  Chis  approach,  however.  Chough 
in  a slighcly  differenC  manner.  In  Chis  procedure,  Che  burning  race  is 
calculaCed  by  dividing  Che  oxidizer  parcicle  diameCer  by  the  sum  of  Che 
Cimes  required  for  each  increment  until  Che  bottom  of  the  oxidizer  particle 
is  reached: 


r - D H TAU  - D /![  (A(DELOX))  ./r  .]  (47) 

avg  o j j o j j ox,j 

TAU,  ■ time  for  regression  plane  to  move  from  one  increment 
J to  the  next;  summation  carried  out  only  to  increment 
at  which  oxidizer  burns  out 


DELOX  ■ distance  of  center  of  oxidizer  surface  from  initial 

J oxidizer  peak  at  jth  calculation  increment  (See  Figure 
III-3. ) 


r 


ox,J 


linear  oxidizer  regression  rate  in  jth  calculation 
Increment 


This  procedure  begs  the  question  of  what  happens  to  the  fuel  "left  over" 
when  the  oxidizer  particle  burns  out  before  the  fuel.  Physically  we  can 
perhaps  just  assume  that  it  somehow  flakes  off.  This  question  needs  to  be 
addressed  further.  This  third  approach  is  intuitively  more  appealing  to 
this  author  in  that  it  allows  for  the  fact  that  the  particle  will  spend 
more  time  at  regression  increments  where  the  burning  rate  is  lower,  while 
the  first  two  procedures  Involve  an  implicit  assumption  that  each  of  the 
regression  Increments  is  equally  likely. 

To  date,  only  preliminary  calculations  have  been  run  with  this 

model,  with  no  attempt  at  optimizing  the  values  of  such  parameters  as  K^, 

Kn_,  B , E , B,  , and  E.  , . These  calculations  have  been  run  only  for 
OF  ox  ox  fuel  fuel 

Formulation  4525,  the  73/27  AP/HTPB  propellant  with  20  micron  diameter 
unlmodal  oxidizer  used  as  the  baseline  formulation  in  the  experimental 
phase  of  this  program.  The  equations  used  for  thermal  conductivity  and  gas 
diffusivlty  were: 
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5.5*10  & cal/cmaec  °K, 


T ln°K 


(48) 


-5  3/2  2 

D ■ 1.0*10  T ' /P  cm  /sec. 


T In  PK,  P in  atm. 


(49) 


The  activiation  energies  for  the  fuel  vaporization  and  oxidizer  sublimation, 

Ef  . and  E .were  chosen  as  16,900  cal/mole  and  22,000  cal/mole,  respectively, 
tuei  ox 

the  former  value  based  on  data  of  Cohen,  Fleming,  and  Derrv  ' and  the 

(35) 

The  pre-exponential 


latter  based  on  data  of  Andersen  and  Chaiken. 
factor  in  the  oxidizer  sublimation  expression,  B 


ox 


was  set  at  2*l(r  gm/ 


cm  sec  based  on  the  BDP  analysis  of  pure  ammonium  perchlorate  deflagration, 

while  Bfue2  was  set  at  5500  gm/cm  sec,  this  value  being  chosen  to  force 

the  linear  regression  rates  of  oxidizer  and  fuel  to  be  equal  at  Tg  * 900°K 

(somewhat  arbitrary,  to  say  the  least).  Based  on  the  BDP  analysis  for  pure 

ammonium  perchlorate  deflagration,  K^,  (Eq'n.  36)  was  chosen  to  be  0.0001 

to  0.0002  atm  sec.  The  diffusion  flame  reaction  distance  constant  KQF 

(Eq'n.  35)  was  chosen  to  yield  a reaction  rate  approximately  ten  times 

higher  than  the  AP  monopropellant  reaction  rate  at  T - 2000  °K,  resulting 

-13  -2 

in  a range  of  values  of  (0.3  - 0.5)*10  atm  sec  #K  (again,  rather 
arbitrarily).  Predicted  burning  rate  versus  pressure  curves  at  0 and  1000 
ft/sec  crossfljw  velocities  for  various  combinations  of  these  parameters 
are  plotted  in  Figure  II-5  along  with  data  for  Formulation  4525.  As  may  be 
seen,  agreement  between  theory  and  experiment  is  surprisingly  good.  However, 
it  is  relized  that  this  may  be  fortuitous  and  that  application  of  the  model 
to  other  formulations  may  not  give  such  satisfactory  results. 


4.  Examination  of  Effect  of  Postulated  Flow  Profiles 

During  the  course  of  this  program,  the  author  became  aware 
of  complaints  that  data  on  erosive  burning  taken  in  test  devices  where 
driver  grain  product  gases  were  passed  over  small  specimens  (strips 
or  tablets)  of  the  test  propellant  did  not  extrapolate  well  to  motor  con- 
ditions. the  erosive  effects  being  considerably  less  in  actual  motors  than 
anticipated  from  the  laboratory  results.  One  possible  explanation  for  this 
is  that  the  boundary  layer  flow  profiles  are  considerably  different  in  the 
test  device  flow  channel  than  in  a motor.  In  most  test  devices,  including 
the  one  used  in  this  program,  the  ratio  of  blowing  velocity  (gas  velocity 
normal  to  the  propellant  surface,  generated  by  the  combustion)  to  crossflow 
velocity  is  quite  small,  less  than  0.05,  lying  in  a range  vhc  s the  data 
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BURN  RATE  (tn/Mc) 


Figure  1 1-5.  Reiulti  of  Preliminary  Calculations  with  Second  Generation  Model, 
Formulation  4525  (73/27  AP/HTPB,  20  micron  Diameter  AP). 
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(31) 

of  Mickley  and  Davis  used  in  both  of  the  models  described  above  are 

(37)  (38) 

applicable.  Recent  work  by  Yamada,  ejt  al . and  Dunlap,  et  al.,  ' 

however,  indicates  that  in  cylindrically  perforated  motors,  where  the  ratio 
of  blowing  velocity  to  crossflow  velocity  tends  to  be  much  higher  (except 
at  the  aft  end  of  very  long  grains) , the  flow  profiles  are  considerably 
different,  approximating  those  of  an  lnviscld  flow  with  a no-slip  wall 
boundary  condition.  In  tb^s  case,  the  axial  velocity  flow  profile  is  given 
by  (referring  back  to  the  nomenclature  of  the  first  generation  model  re- 
grading diffusion  distance): 

. «.  J.  /;d/2)  - W1"  8 

“crossflow, y - sin  6 2 C°*  I 2 \ (D/2) 

u ■ average  crossflow  velocity 
D ■ flow  channel  diameter 

The  first  generation  model  was  modified  to  use  this  profile  rather  than 
the  one  described  by  Equations  ’ to  40  in  the  development  of  the  second 
generation  model  (which  was  also  k ed  in  the  original  version  of  the  first 
generation  model).  A set  of  calculations  was  then  run  for  a motor  with 
a port  diameter  of  1.2  inches  using  both  types  of  profiles  for  comparison. 
Formulation  4525  (73/27  AP/HTPB,  20  micron  diameter  AP)  was  used  for  these 
predictions  since,  as  will  be  shown  later,  good  agreement  was  found  between 
the  Generation  1 Model  using  the  Mlckley-Davis  profiles  and  data  taken  in  our 
test  apparatus  with  this  propellant.  Result*  of  these  calculations  are 
shown  in  Figure  II-6.  As  may  be  seen,  replacement  of  the  Mlckley-Davis 
profiles  with  the  lnviscld  no-slip  profiles  results  in  a considerable  reduc- 
tlonln  the  predicted  degree  of  erosive  burning.  This  is  a particularly 
lmportanc  result,  pointing  out  the  necessity  of  correct  definition  of  flow 
profiles  in  a given  motor  configuration  for  accurate  prediction  of  erosive 

P 

burning.  Thus  it  appears  that  further  attention  need  be  paid  to  accurate 
definition  of  profiles,  not  only  in  cylindrically  perforated  motors,  but  in 
wagon-wheel  perforations,  star  configurations,  and  any  other  configurations 
where  it  is  felt  that  erosive  burning  may  be  Important.  Modification  of 
the  first  generation  model  (and  the  second  generation  model)  to  accomodate 
various  flow  profiles  Is  quite  simple,  the  modification  of  the  first  genera- 
tion model  to  use  lnviscld  no-slip  profiles  having  taken  less  than  one 
man-day. 
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Figure  11-6.  Erosive  Burning  Rite  Predictions  with  the  First  Generation  Model  Using 
Mickley-Oavis  Boundary  Layer  Profiles  and  Using  Inviscid  No-Slip  Wall 
Boundary  Layer  Profiles.  Formulation  4525  (73/27  AP/HTPB,  20  micron  Diameter  AP). 


c. 


EXPERIMENTAL 


1.  Equipment 

In  the  past,  two  basically  different  approaches,  each  with 
advantages  and  disadvantages,  have  been  pursued  for  measurement  of  steady- 
state  erosive  burning  rates.  The  first  Involves  determination  of  ablation 
rates  of  grains  in  rocket  motors  under  actual  firing  conditions,  while  the 
second  involves  measurement  of  burning  rates  of  propellant  specimens  (tablets 
or  strips)  located  in  a stream  of  hot  flowing  gas  external  to  a motor. 

The  first  method  has  the  advantage  that  the  erosive  burning  process  is  being 
studied  under  actual  motor  test  conditions,  but  the  disadvantage  that  these 
conditions  are  difficult  to  determine  accurately.  Pressure  and  velocity 
vary  through  the  chamber  and  also  change  with  time.  Such  burn  rate  measure- 
ment techniques  as  interrupted  burning  must  use  time-averaged  values,  and 
the  time  periods  must  be  relatively  long.  Continuous  measurement  of  burning 
rates  within  a rocket  by  x-ray  requires  special  elaborate  equipment.  The 
use  of  probes  to  measure  burning  rates  and  pressures  is  difficult  because 
many  probes  are  required,  and  also  runs  the  risk  of  interfering  with  the 
chamber  flow  and  disturbing  normal  burning. 

The  major  disadvantage  of  the  second  method  of  studying  erosive 
burning  lies  in  possible  differences  in  detailed  flow  characteristics  for 
flow  over  stripb  or  buttons  of  propellant,  as  compared  to  flow  adjacent  to 
a grain  surface  in  a motor,  even  at  the  same  freestream  velocity  and  pressure. 
(In  fact,  while  flow  of  a hot  gas  stream  from  a gas  generator  over  a pro- 
pellant strip  or  button  preceded  by  some  sort  of  flow  channel  will  tend  to 
develop  turbulent  boundary  layer  type  profiles,  the  recent  work  by  Yamada, 
et  al.,  and  Dunlap,  et  al. , ' indicates  that  the  flow  profiles  in  typical 

cylindrically  perforated  motors  approximate  those  of  an  inviscid  flow  with 
a no-slip  wall  boundary  condition,  as  discussed  in  the  preceding  section 
of  this  report.)  Advantages  are  that  such  parameters  as  freestream  crossflow 
velocity  and  pressure  can  be  easily  controlled  and  measured.  Most  important, 
continuous  measurement  of  burning  rate  can  be  accomplished  with  high  speed 
visual  photography  through  quartz  windows  in  the  apparatus.  Hence,  results 
obtained  using  this  approach  are  likely  to  be  more  accurate  and  reliable 
than  those  from  motor  firings. 
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Previous  studies  have,  In  general,  not  resulted  In  Instantaneous 
measurements  of  erosive  burning  rates  under  well-characterized  local  flow 
conditions.  This,  coupled  with  a lack  of  data  In  the  high  Mach  number 
region  (M  ^ 0.5)  and  a lack  of  studies  In  which  propellant  parameters  are 
systematically  varied  one  at  a time  under  Identical  hydrodynamic  conditions 
has  resulted  In  an  incorplete  understanding  of  erosive  burning  phenomena. 

The  experimental  test  apparatus  and  procedures  employed  In  this 
study  of  erosive  burning  are  described  in  some  detail  in  Reference  12. 

A schematic  of  the  basic  test  apparatus  is  presented  as  Figure  II-7.  A 
cylindrically  perforated  6C4  driver  grain  (15.2  cm  outside  diameter,  10.2 
cm  inside  diameter)  whose  length  is  chosen  to  give  the  desired  operating 
pressure  for  a given  test,  produces  a high  velocity  gas  flow  through  a 
transition  section  into  a rectangular  test  section  which  contains  the  test 
grain  (generally  the  same  formulation  as  the  driver  grain).  The  contoured 
transition  section  is  approximately  10  cm  (4  inches)  long.  The  test  grain 
extends  from  the  test  section  back  through  the  transition  section  to  butt 
against  the  driver  grain  in  order  to  eliminate  leading  edge  effects  which  would 
be  associated  with  a test  grain  standing  alone.  The  test  grain  is  approxi- 
mately 30  cm  (12  inches)  long  (plus  the  10  cm  extending  through  the  transi- 
tion section)  by  1.90  x 2.50  cm  (3/4  inch  and  1 inch)  web  and  burns  only 
on  the  1.90  cm  face.  The  flow  channel  of  the  test  section  is  initially 
1.90  cm  x 1.90  cm  (3/4  inch  x 3/4  inch),  opening  up  to  1.90  cm  x 4.45  cm 
(3/4  inch  x 1-3/4  inch)  as  the  test  propellant  burns  back  through  its 
2.54  cm  (1  inch)  web.  For  high  Mach  number  tests,  the  apparatus  was  operated 
in  a nozzlelesB  mode  with  the  gases  choking  at  or  near  the  end  of  the  test 
grain,  while  for  lower  Mach  Number  tests,  a 2-dimensional  nozzle  was  installed 
at  the  end  of  the  test  channel. 

During  each  test,  pressure  and  crossflow  velocity  varied  with  time 
and  location  along  the  test  grain.  (For  the  nozzleless  tests,  pressure 
varied  significantly  with  time  and  location,  while  crossflow  velocity 
varied  considerably  with  location  but  not  significantly  with  time.  For 
tests  using  a nozzle  with  an  initial  port  to  throat  area  ratio  of  1.5  or 
higher,  on  the  other  hand,  pressure  did  not  vary  strongly  with  location  but 
did  rise  with  time  due  to  the  progresslvity  of  the  driver  grain,  while 
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Figure  11-7.  Sketch  of  Test  Hardware. 


crossflow  velocity  varied  strongly  with  tine  and  slightly  with  location.) 
These  variations  permitted  design  of  teats  to  yield  considerable  burning 
rate-pressure-crossflow  velocity  data  in  relatively  few  tests,  provided  that 
these  parameters  could  be  measured  continuously  at  several  locations  along 
the  test  grain.  These  parameters  were  measured  in  the  following  manner. 

The  burning  rate  was  directly  measured  by  phot  ographing  the 
ablating  grain  with  a high-speed  motion  picture  camera  through  a series 
of  four  quartz  windows  located  along  the  length  of  the  test  section.  (See 

Figure  II-7  and  II-8.)  Frame  by  frame  analysis  of  the  films  allowed 
determination  of  instantaneous  burning  rate  as  a function  of  time  at  each 
of  the  four  window  locations.  The  windows  were  flush-mounted  on  the  inside 
of  one  side  wall  and  sealed  with  "0"-ring  seals.  A detailed  drawing  of  the 
test  section,  emphasizing  the  viewing  window  layout  is  shown  in  Figure  II-8. 
The  flow  gap  was  surrounded  by  propellant  on  one  side  and  asbestos  phenolic 
on  the  other  three  sides,  with  circular  cutouts  in  the  asbestos  phenolic 
through  which  the  inner  part  of  the  windows  butted  flush  against  the  pro- 
pellant. (Inner  and  outer  window  sections  were  used,  since  under  the  more 
severe  test  conditions,  the  inner  window  surface  suffers  damage  during  the 
test  or  during  post-test  cooldown.)  To  date,  tests  with  either  quartz 
or  RTV  potting  compound  inner  windows  have  been  satisfactory  as  long  as  the 
side  of  the  test  propellant  was  adequately  inhibited  to  prevent  any  side- 
burning,  and  quite  satisfactory  films  have  been  obtained. 

For  nozzled  cases,  the  measured  location  of  the  burning  propellant 
surface  at  each  window  as  a function  of  time,  together  with  the  known 
constant  throat  area,  permitted  straightforward  calculation  of  the  cross- 
flow  velocity  as  a function  of  time.  However,  the  very  sensitive  dependence 
of  Mach  Number  on  area  ratio  for  M > 0.5  made  calculation  of  crossflow 
velocity  from  area  ratio  measurement  quite  poor  for  nozzleless  cases. 
Accordingly,  for  these  tests,  stagnation  pressure  was  measured  at  the  aft 
end  of  the  test  section  and  used  in  combination  with  the  measured  driver 
chamber  pressure  for  calculation  of  the  stagnation  pressure  in  the  test 
section  as  a function  of  time  and  position.  (Static  pressure  wall  taps 
at  each  window  location  were  used  for  measurement  of  static  pressure  as  a 
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Figure  11-8.  Test  Section  Detail,  Emphasizing  Viewing  Windows. 


function  of  time  for  both  nozzled  and  nozzleless  cases.)  From  the  static 
and  stagnation  pressure  values  determined  as  a function  of  time  and  posi- 
tion down  the  test  section,  crossflow  Mach  Number  and  velocity  were  calcu- 
lated as  a function  of  time  at  each  window  location  in  the  test  section 
for  the  nozzleless  cases. 

2.  Test  Matrix 

The  rationale  of  the  experimental  part  of  this  program  was  to 
measure  the  erosive  burning  characteristics! over  a wide  range  of  pressure 
and  crossflow  velocity,  of  a series  of  propellants  in  which  various  formu- 
lation parameters  were  systematically  varied.  A series  of  14  formulations 
(Table  II-IV)  was  originally  selected  for  investigation,  with  10  tests  to 
be  conducted  with  the  first  formulation,  and  5 tests  each  with  the  succeeding 
formulations.  During  the  course  of  this  program,  seven  of  these  formula- 
tions (1-5,  7 and  8)  were  characterized.  Additional  study  of  the  matrix 
will  be  carried  out  in  a follow-on  program. 

The  first  six  formulations  listed  in  Table  II-IV  are  "scholastic" 
formulations.  (These  are  referred  to  as  "scholastic"  formulations  in  that 
they  are  formulations  specifically  chosen  to  permit  systematic  variation  of 
well-defined  composition  and  ingredient-size  parameters,  including  the  use 
of  unlmodal  ammonium  perchlorate  particle  size,  but  as  a consequence  are 
not  formulations  being  currently  considered  for  mission  applications.) 

It  was  considered  that  the  use  of  unlmodal  oxidizer  in  early  testing  was 
important,  since  any  model  permitting  prediction  of  burning  rate-pressure- 
crossflow  velocity  characteristics  from  first  principles  will  almost  cer- 
tainly be  first  derived  for  unlmodal  oxidizer.  (Methods  of  handling  multi- 
modal oxidizer  sizes  for  predictions  of  burning  rates  even  in  the  absence 
of  crossflows  are  still  the  subject  of  considerable  debate.)  Formulation  1 
(also  referred  to  elsewhere  as  Formulation  4525)  is  a baseline  1667°K 
HTPB  formulation  (73/27  AP/HTPB)  for  the  initial  test  series.  Formulations 
2,  3,  and  4 (5051,  4685  and  4869)  were  selected  for  investigation  of  the 
interrelated  effects  of  oxidizer  particle  size  and  base  (no  crossflow) 
burning  rate.  Formulations  1 and  4 are  essentially  identical  except  for 
use  of  burning  rate  catalyst  to  change  base  burning  rate.  Formulations 
2 and  3 differ  from  Formulation  1 in  oxidizer  particle  size  (200,  5 and 


51 


w a 

5!  - 

* 

u u 

u fl  9 

u 9 to 

v u m 

«M  « M 

. s a j 

u i»  9 


kJ  u 

0*  o 

* < v 

M *4 

to  k- 

O to) 

*4 

a.  at 

< N A 


I 1 


X X X X 

to  w to  to  to. 

4*  -to  <to  -4  < 


• AX  H* 

J 1*  « 

A 0 v i 

H • to-  to  to4 

!s  c»  “ 

£5  >5  a 

to  A • to 

A to  to  *-4  0 

u A 9 to 

A *Q  5 • 

AC  A Q ^ 

X ^ to  * 

^ « A X 

0 dto  » 44 

fcS  f°.“  > 

r Hsc  . 

~ I «n  to 

•A  3 X «M  A 

AX  P. 


■ ■ ■ |S  S£  ? « i 

H Hn^  I2h 


cto  m 

U ^ H 

o o 3 u.  o 

n m v n 


•H  0Q 

2 

54  s 


£ £ 5 fc 

E E 5,  « 

5!  23  n 25 

r-.r-.%0fA 
N N <M  N 

V X N *S 

<n  m n n 


k ^ ^ to 

A lA  «H  (A  0N  «# 

X <N  »A  dp  NO  CM 

fl  m o <o  oq  ia 

3 ia  -*  -e  »a 

Z w w w w 

H N w *•  m 


« « s 

9 • s. 

X£  £ 

u 


Hs 

Z 

- « s 

■ u *> 

32.  ■? 

i!  I 


3 * 

la  a 


ia  <n 
«s  M «to 

5|  5 

O CO 


A p A 

ill 


- 0 to  0 

a os  h « 

m am 


X 

X 

H X 

A 

• A 

< A 

to 

H to 

O to 

a 

Si 

~ % S5! 

<N  H 

^ H >»  H 

ftj  4 0]  4 

li  1? 


I i 

4i  to 

1 I 

f f 


52 


20  micron  AP,  respectively),  and  as  a consequence,  also  In  base  burning 
rate.  Comparison  of  results  from  tests  with  these  four  formulations  should 
permit  isolation  of  the  oxidizer  particle  size  and  base  burning  rate  effects 
on  sensitivity  of  propellant  burning  rate  to  crossflow. 

In  terms  of  independent  variables.  Formulation  5 (5542T)  differs 
from  Formulation  1 only  in  terms  of  oxidizer/fuel  ratio  (>7/23  versus  73/ 
27),  yielding  a higher  flame  temperature  (20656K  vs.  1667°K)  and  a different 
burning  rate,  oxidizer  size  being  held  constant.  Thus  comparison  of  the 
results  for  these  formulation*  permits  definition  of  the  effect  of  oxidizer/ 
fuel  ratio  change  at  constant  oxidizer  particle  size.  With  Formulation  7 
(5565T),  on  the  other  hand,  oxidizer/fuel  ratio  is  varied  from  that  of 
Formulation  1 (82/18  vs.  73/27),  but  oxidizer  sizes  are  changed  to  give 
approximately  the  same  zero-crossflow  burning  rate-pressure  curve  for  the 
two  formulations  - this  permits  examination  of  the  effect  of  varying  oxi- 
dizer/fuel (and  thus  flame  temperature)  at  constant  base  burning  rate. 
Formulation  8 (5555T)  is  identical  to  Formulation  7 except  for  use  of  much 
finer  oxidizer  to  yield  higher  base  burning  rate. 

A total  of  45  tests  were  carried  out  with  these  seven  formulations. 
Of  these,  39  yielded  useful  data,  while  six  were  failures  due  to  breakup 
of  the  test  grain  (in  nozzleless  tests)  or  due  to  camera  failure.  The 
rationale  and  ballistics  analyses  used  in  selecting  the  test  conditions 
employed  were  discussed  in  detail  in  Reference  12.  Basically,  the  first 
three  tests  were  designed  to  yield  erosive  burning  data  for  Formulation  1 
over  a range  of  crossflow  velocities  of  180  to  350  m/sec  (600  to  1200 
ft/sec)  and  a range  of  pressures  of  1.4  to  8.2  MPa  (200  to  1200  Psla). 

The  next  three  tests  were  chosen  to  examine  the  same  formulation  over  a 
crossflow  velocity  range  of  approximately  600  to  850  m/sec  (2000  to  2800 
ft/sec)  and  a pressure  range  of  1 to  5 MPa  (150  to  750  psia).  Tests  7 and 
8 differed  from  Tests  1 and  3 only  in  having  no  test  grain  in  the  transi- 
tion section.  These  tests  were  aimed  at  determining  the  sensitivity  of 
erosive  burning  to  major  upstream  geometry  changes.  Tests  9 and  10  differed 
from  Tests  1 and  3 only  in  their  use  of  a hotter  (2400*K)  driver  formulation 
with  the  baseline  test  formulation  (1667*K  flame  temperature).  The  purpose 
of  these  tests  was  to  determine  whether  the  "core"  gas  temperature  affected 
the  erosive  burning  of  a given  formulation.  (Recall  that  Marklund  and 
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Lake^  found  that  It  did  not.)  Test  11  - 15,  16  - 20,  21  - 25,  26  - 30, 
31-35  and  36  - 40  were  designed  to  be  analogous  to  Test  1-5  with  re- 
placement of  Formulation  1 (4525)  by  Formulations  2,  3,  4,  5,  7 and  8 
(5051,  4685,  4869,  5542T,  5565T,  and  5555T).  Tests  41  - 45  were  added  to 
fill  in  data  gaps  revealed  by  earlier  tests. 

3.  Test  Results  and  Discussion 

Measured  burning  rate  augmentation  ratios  (e  * r/rQ)  f°r  two 
typical  nozzled  tests  are  plotted  against  time  in  Figure  II-9  and  11-10. 

For  these  tests,  the  augmentation  ratio  at  any  given  time  should  be  nearly 
the  same  from  window  to  window  (since  pressure  and  crossflow  velocity  do 
not  vary  significantly  with  location),  unless  boundary  layer  development 
effects  are  playing  an  important  role.  As  may  be  seen,  the  variation  from 
window  to  window  is  actually  fairly  small.  The  augmentation  ratio  versus 
time  predicted  with  the  first  generation  model  described  earlier,  using 
the  measured  pressure  and  crossflow  velocity  versus  time  data,  is  also 
presented  in  these  figures.  The  agreement  between  theory  and  experiment 
is  seen  to  be  reasonably  good. 

Results  of  the  tests  made  fox  study  of  the  effect  of  upstream  flow 
conditions  (two  tests  conducted  at  essentially  identical  conditions  to 
tests  in  the  main  test  series,  except  for  the  absence  of  test  grain  in 
the  transition  section)  are  presented  in  Figure  11-11  and  11-12.  As  may 
be  seen,  the  effects  of  the  upstream  flow  change  were  quite  small,  the 
differences  in  burning  rate  augmentation  ratio  between  corresponding  tests 
varying  essentially  only  to  the  degree  predicted  by  the  slight  difference 
in  pressure-crossf low  velocity-time  history  in  the  matched  tests.  Accordingly, 
it  is  concluded  that  the  erosive  burning  measured  at  the  viewing  ports  is 
not  particularly  sensitive  to  the  driver  grain-transition  section  contours 
in  the  test  apparatus.  This  result  is  consistent  with  the  observation  that 
the  augmentation  rates  do  not  vary  significantly  with  window  location 
for  the  nozzled  tests  (where  pressure  and  crossflow  velocity  are  nearly 
the  same  at  each  window  location  at  any  given  time). 

As  discussed  earlier,  erosive  burning  models  based  upon  increased 
heat  transfer  from  a "core"  gas  flow  accompanying  crossflow,  predict  that 
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Figure  11*9.  Burning  Rate  Augmentation  Ratio  Versus  Time  for  Test  Condition  No.  1, 
Formulation  No.  1 (4525),  Nozzled  Test,  Low  Pressure. 


Figure  11*10.  Burning  Rate  Augmentation  Ratio  Versus  Time  for  Test  Condition  No. 

Formulation  No.  1 (4525)  - Nozzled  Test,  High  Pressure.  No  Grain  in  T 


Figure  11-11.  Comparison  of  Erosive  Burning  with  and  Without  Test  Grain  Extending  Thru 
Transition  Section  to  Mata  with  Driver  Grain  - Low  Pressure. 


Figure  11-12.  Comparison  of  Erosive  Burning  with  and  Without  Test  Grain  Extending  Thru 
Transition  Section  to  Mate  with  Driver  Grain  - High  Pressure. 


with  a given  test  section  propellant,  variation  of  the  flame  temperature 
of  the  driver  propellant  should  lead  to  variation  in  the  erosive  burning 
augmentation  ratio  at  fixed  crossflow  velocity  and  pressure.  Two  pairs 
of  tests  (land  9,  3 and  10), in  which  the  driver  grain  flame  temperature 
was  varied  from  1667 °K  to  2425°K,  while  the  test  section  propellant  was 
held  constant  and  crossflow  velocity  and  pressure  versus  time  histories  were 
held  as  nearly  constant  as  possible,  were  carried  out  in  this  program. 
Results  are  presented  in  Figures  11-13  and  11-14.  In  each  figure,  measured 
burning  rate  augmentation  ratio  and  the  ratio  predicted  using  the  first 
generation  model  described  earlier  are  plotted  against  time  for  each  of 
the  paired  tests.  (The  predicted  values  are  presented  to  permit  a zeroing 
out  of  the  slight  differences  in  pressure  and  crossflow  velocity  versus  time 
histories  of  the  paired  tests.)  The  different  "core"  gas  temperatures  in 
the  paired  tests  are  seen  to  have  negligible  effect  on  the  erosive  burning 
characteristics  of  the  test  propellant.  This  result  is  consistent  with  the 
results  of  Marklund  and  Lake^  and  casts  further  doubt  on  the  validity  of 
models  in  Category  1 of  Table  II-II  (most  notably,  the  commonly  used  model 
of  Lenoir  and  Robillard). 

A rather  complete  set  of  data,  covering  a pressure  range  of 
1 to  5 MPa  (10  to  50  atmospheres)  and  a crossflow  velocity  range  of  180 
to  670  m/sec  (600  to  2200  ft/sec)  has  been  obtained  for  Formulation  4525. 
(This  is  the  baseline  HTPB  formulation,  containing  73  weight  percent  20 
micron  diameter  ammonium  perchlorate  and  27  weight  percent  hydroxyterminated 
polybutadiene  (HTPB)  binder,  with  a trace  of  carbon  black  added  to  opacify 
the  propellant.)  Experimental  results  and  theoretical  predictions  are 
presented  in  Figure  11-15  and  11-16.  As  may  be  seen,  agreement  between 
prediction  and  data,  while  not  as  good  as  with  the  Saderholm  propellant,  is 
reasonably  good.  The  predicted  curves  for  burning  rate  versus  pressure  at 
various  crossflow  velocities  (Figure  11-15)  do  seem  to  group  more  tightly 
than  the  data.  That  is,  as  shown  more  clearly  in  Figure  11-16,  the  model 
tends  to  slightly  overpredict  the  burning  rate  at  low  crossflow  velocities 
and  slightly  underpredict  it  at  high  velocities. 

Theoretical  predictions  and  experimental  measurements  of  erosive 
burning  rates  for  Formulations  5051,  4685,  4869,  5542T,  5565T,  and  5555T 
are  presented  in  Figures  11-17  through  11-22.  Formulation  5051,  which 
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Figure  11-14.  Comparison  of  Erosiva  Burning  of  Formulation  1 (4525)  with  1667°K 
Driver  Grain  and  2425°K  Driver  Grain  - High  Pressure. 
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Figure  11-16.  Burning  Rate  Versus  Crossflow  Velocity  Data  and  Predictions  for 
Formulation  4525  (73/27  AP/HTPB,  20  Micron  AP). 


PRESSURE  (MPa) 


(091/uio)  3iva  9NINU0S 


o 

o 

e 

in 

o 

■n 

o 

in 

CM 

CM 

P-’ 

o 

d 

o 

m 


O 

u 

t 

O 

Q 

< 

H 

< 

Q 


OUUU  IUUOO 

.8  .8  81  ^ a .8 

'e'e^I'e'e  e e e 

o«nnnoooo 

NNN<-'-IAIONtD 
II  II  H II  II  II  II  II  II 
>>>>>>>>> 

♦ X0300CIQO 


=b 

O 


o 

00 

d 


o 

to 


o 

CM 

o 


CM 


a* 

-Zo 

ew 

o«» 


■m 

(To 

• in 

Msr 

|l< 

a 3 C 

’ £ ° 

• r t. 

£ og 

00  • . 

• 53* 

Jof: 

C oh 

psr 

lii 

utsd 

_ S co 
■g  2^ 

• u c- 

« 3.2 

.3  o-g 

s5| 

8 >i 


3 

o» 


2 8 8 

d d d 


(3m/ui)  3ivu  ONiNbna 


64 


PRESSURE  (MPa) 


(oas/uiD)  3lVb  ONINdnS 


ooo  o o i/tom 


(sat/uj)  31VH  ONINdOS 


65 


I 

f 

i 

'4'*- 


differs  from  the  baseline  formulation  through  use  of  200  micron  AP  oxidizer 
in  place  of  20  micron  oxidizer,  is  predicted  to  be  somewhat  more  sensitive 
to  crossflow  than  the  baseline  formulations.  Except  at  low  pressure  and 
very  high  crossflow  velocities,  agreement  between  predicted  and  measured 
augmentation  ratio  is  fairly  good.  At  low  pressure  and  high  crossflow 
velocity,  however,  the  measured  burning  rates  considerably  exceed  the  pre- 
dicted values.  As  shown  in  Figure  11-18,  Formulation  4685,  which  differs 
from  the  baseline  formulation  by  replacement  of  20  micron  oxidizer  with 
5 micron  oxidizer,  exhibits  considerably  less  sensitivity  to  erosion  than 
that  baseline  formulation,  as  predicted.  Agreement  between  predicted  and 
observed  burning  rates  appears  to  be  good,  except,  again,  in  the  low  pressure, 
high  crossflow  velocity  region  (less  than  2 MPa  or  20  atmospheres,  greater 
than  300  to  600  m/sec  or  1000  to  2000  ft/sec  crossflow  velocity).  Break- 
down of  the  model  presented  herein  in*  this  pressure-crossflow  velocity  region 
is  not  unexpected  since,  in  this  region,  the  composite  propellant  begins 
to  behave  more  like  a homogeneous  propellant  than  a heterogeneous  propellant, 
and  the  model  only  considers  effects  of  crossflow  on  the  diffusional  mixing 
processes  of  oxidizer  and  fuel  streams.  In  order  for  the  model  to  be  useful 
in  low  pressure,  high  crossflow  velocity  regions,  it  appears  that  an  addi- 
tional mechanism  beyond  that  of  flame-bending  must  be  invoked.  With 
Formulation  4869  (Figure  11-19),  which  differs  from  the  baseline  formulation 
through  addition  of  two  percent  iron  oxide  catalyst,  data  and  theoretical 
predictions  agree  fairly  well  at  high  crossflow  velocities,  but  not  nearly 
as  well  at  low  crossflow  velocities  where  the  predictions  of  erosive  burning 
rate  augmentation  are  somewhat  higher  than  observed  in  the  experiments. 

An  explanation  of  this  discrepancy  has  not  yet  been  developed. 

With  Formulation  5542T  (analogous  to  the  baseline  formulation  but 
with  higher  oxidizer/fuel  ratio  and  consequently  higher  temperature  and  base 
burning  rate,  oxidizer  size  being  held  constant)  the  sensitivity  to  cross- 
flow  appears  to  be  somewhat  lower  than  predicted  (Figure  11-20)  though  the 
degree  of  disagreement  between  data  and  theory  is  not  gross.  The  data 
obtained  for  Formulation  5565T  (with  approximately  the  same  zero  crossflow 
burning  rate-pressure  behavior  as  the  baseline  formulation,  but  a considerably 
higher  oxidizer/fuel  ratio  and  flame  temperature)  presented  in  Figure  11-21 
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Figure  11-19.  Theoretical  and  Experimental  Burn  Rate  • Pressure  Relationship: 

for  Various  Crossflow  Velocities  for  Formulation  4869  (1667°K 
Formulation,  76/22/2  AP/HTPB/Fe203,  20  Micron  AP). 


BURNING  RATE  (irt/tec) 


PRESSURE  (MPa) 


PRESSURE  (atm) 


Figure  11-20.  Theoretical  and  Experimental  Burn  Rate  • Pressure  Relationships 
for  Various  Crossflow  Velocities  for  Formulation  5S42T  (2065  K 
Formulation,  77/23  AP/HTPB,  20  Micron  AP). 
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Figure  11*21.  Theoretical  and  Experimental  Burn  Rate  - Preuure  Relationships  for 

Various  Crossflow  Velocities  for  Formulation  5565T  (2575°K  Formulation, 
82/18  AP/HTPB,  Bimodal  With  Sizes  Chosen  to  Match  4525  Burning  Rate). 
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are  somewhat  limited,  but  indicate  reasonable  agreement  with  theory,  the 
formulation  being  quite  sensitive  to  crossflows.  Formulation  5555T 
(Figure  11-22),  a high  burning  rate  formulation  is  predicted  to  be  rather 
insensitive  to  crossflows:  the  data  corroborates  this  prediction. 

Next,  let  us  compare  results  for  the  various  formulations  to 
Identify  parameters  which  influence  the  sensitivity  of  composite  propellants 
to  crossflows.  Between  Formulations  4525,  5051,  and  4685,  the  only  indi- 
pendent  variable  changed  is  the  oxidizer  particle  size,  composition  being 
held  constant.  The  change  of  oxidizer  size,  of  course,  leads  to  a change 
in  base  (no  crossflow)  burning  rate  versus  pressure  characteristics. 
Formulation  5051,  containing  200  micron  diameter  AP  is  the  slowest  burning 
of  the  three  formulations,  with  Formulation  4685  (5  micron  AP)  being  the 
fastest  and  Formulation  4525  (20  micron  AP)  being  intermediate.  For 
instance,  at  5 MPa  (50  atmospheres)  the  base  burning  rate  of  5051  is  0.47 
cm/sec,  that  of  4525  is  0.68  cm/sec  and  that  of  4685  is  1.15  cm/sec. 
Examination  of  Figures  TI-15,  11-17,  and  11-18  indicates  that  the  sensitivity 
of  burning  rate  to  crossflow  increases  with  increasing  particle  size 
(decreasing  base  burning  rate).  For  example,  at  a crossflow  velocity  of 
200  m/sec  (650  ft/sec)  and  a pressure  of  5 MPa  (50  atmospheres),  the  aug- 
mentation ratio  for  4685  is  about  1.10,  that  for  4525  is  1.65,  and  that  for 
5051  is  2.0. 

Comparison  of  data  for  4525  and  4869,  two  formulations  of  essen- 
tially the  same  oxidizer/fuel  ratio,  flame  temperature,  and  oxidizer  particle 
size,  with  the  base  burning  rate  being  varied  through  use  of  catalyst  in  4869, 
again  shows  an  increase  in  sensitivity  of  burning  rate  to  crossflow  with  a 
decrease  in  burning  rate.  At  5 MPa  (50  atmospheres)  the  base  burning  rates 
for  4869  and  4525  are  1.40  cm/sec  and  0.68  cm/sec,  respectively.  At  this 
pressure,  with  a crossflow  velocity  of  200  m/ser.  (650  ft/sec)  their  r/rQ 
values  are  1.10  and  1.65  respectively,  while  at  600  m/sec  (1950  ft/sec), 

the  r/r  values  are  1,75  and  2.3.  Thus  base  burning  rate  is  seen  to  affect 
o 

the  erosion  sensitivity  of  composite  propellants  even  at  constant  oxidizer 
particle  size,  erosive  effects  increasing  with  decreasing  base  burning 
rate. 
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Figure  11-22.  Theoretical  and  Experimental  Burn  Rate  • Pressure  Relationships 
for  Various  Crossflow  Velocities  for  Formulation  5555T  (2575  * 
Formulation,  82/11  AP/HTPB,  High  Burn  Rate). 
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BURNING  RATE  (cm/sec) 


Formulations  4685  and  4869  have  approximately  the  same  base  burning 
rate  at  8 MPa  (80  atmospheres)  with  catalyst  and  oxidizer  particle  size 
effects,  on  base  burning  rate  roughly  cancelling.  Thus  comparison  of  erosion 
sensitivity  of  these  formulations  at  this  pressure  is  of  interest  in  that 
oxidizer  particle  size  is  varied  (5  micron  diameter  for  4685,  20  micron 
diameter  for  4869)  while  base  burning  rate  is  held  constant.  Comparison 
of  data  from  Figures  11-18  and  11-19  indicates  that  these  formulations  have 
roughly  the  same  sensitivity  to  the  lower  crossflow  velocities  tested  at 
8 MPa  (80  atmospheres) , with  the  catalyzed  propellant  being  somewhat  more 
sensitive  at  the  higher  crossflow  velocities  tested.  Thus  it  appears  that 
it  is  the  base  burning  rate  rather  than  the  oxidizer  particle  size  per  s 6 
which  dominates  the  sensitivity  of  composite  propellants  to  erosive  burning, 
though  oxidizer  size  does  have  some  further  residual  effect,  erosion  sensi- 
tivity decreasing  with  decreasing  particle  size  at  constant  base  burning  rate. 

Comparison  of  test  results  for  Formulations  4525,  5542T  and  5565T 
permits  study  of  the  effect  of  oxidizer/fuel  ratio  (and  thus  flame  temperature) 
on  erosion  sensitivity,  both  at  constant  oxidizer  particle  size  (5542T 
and  4525)  and  at  constant  base  burning  rate  (5565T  and  4525).  Formulation 
5542T  differed  from  4525  in  oxidizer/fuel  ratio  (77/23  versus  73/27)  and 
consequently  flame  temperature  (2065°K  vs.  1667°K).  Since  the  oxidizer 
particle  size  was  the  same  for  both  propellants,  the  higher  oxidizer/fuel 
ratio  for  5542T  led  to  high  base  burning  rate  (1.14  cm/sec  vs.  0.68  cm/sec 
at  5 MPa).  Study  of  Figures  11-15  and  11-20  reveals  that  the  erosion  sen- 
sitivity of  5542T  is  considerbly  less  thac  that  of  4525  over  the  entire 

range  of  crossflow  velocities  studied  (e.g.,  r/r  * 1.10  for  5542T  and  1.65 

o 

for  4525  at  200  m/sec,  5 MPa;  and  r/rQ  * 1.7  for  554 2T  and  2.9  for  4525  at 
800  m/sec,  5 MPa).  Thus  we  see  that  changing  oxidizer /fuel  ratio  from  very 
fuel-rich  to  less  fuel-rich,  with  accompanying  Increase  in  flame  temperature 
and  burning  rate  leads  to  decreased  sensitivity  to  erosive  burning.  Com- 
parison of  results  lor  5565T  and  4525,  which  differ  in  oxidizer/fuel  ratio 
but  not  in  base  burning  rate  (oxidizer  particle  size  having  been  adjusted 
to  compensate  for  the  burning  rate  change  with  changing  oxidizer/fuel)  permits 
separation  of  the  effects  of  varying  oxidizer/fuel  (and  thus  flame  temperature) 
from  the  effects  of  base  burning  rate.  As  may  be  seen  by  study  of  Figures 
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11-15  and  11-21,  the  sensitivity  of  Formulations  5565T  and  4525  to  crossflow 
are  nearly  the  same.  For  instance,  at  200  m/sec  (650  ft/sec)  crossflow 
velocity  and  5 MPa  (50  atmospheres),  the  augmentation  ratios  for  5565T 
and  4525  are  1.50  and  1.65,  respectively,  while  at  800  m/sec  (2600  ft/sec) 
and  3 MPa  (30  atmospheres),  they  are  2.65  and  2.50.  Accordingly,  we  may 
conclude  that  oxidizer/fuel  ratio  (and  consequently  flame  temperature) 
does  not  directly  afr,cr  ’•he  erosion  sensitivity  of  the  composition  studied 
to  date,  but  only  a' ; through  its  effect  on  base  burning  rate. 

Formulation.-'  5555T  and  5565T  had  the  same  composition,  differing 

only  in  oxidizer  particle  size,  which  was  adjusted  in  5555T  to  give  a very 

high  burning  rate.  Again,  the  effect  on  erosion  sensitivity  of  Increased 

hase  burning  rate  can  be  seen  in  comparison  of  Figures  11-21  and  11-22. 

At  5 MI’a  (50  atmospheres),  the  base  burning  rates  of  5555T  and  5565T  are 

2.94  and  0.70  cm/sec,  respectively.  At  200  m/sec  (650  ft/sec)  crossflow 

velocities,  the  respective  values  of  r/r  are  1.0  and  1.5,  while  at  700  m/sec 

o 

(2300  ft /sec),  they  are  1.2  and  2.4.  Thus,  once  again,  erosion  sensitivity 
is  seen  to  decrease  with  increasing  base  burning  rate. 

D.  Conclusions 

An  experimental  apparatus  for  measurement  of  erosive  burning  rates 
over  a wide  range  of  crossflow  velocities,  up  to  Mach  1 has  been  designed, 
constructed  and  checked-out.  Erosive  burning  characteristics  of  seven  formu- 
lations, with  systematically  varied  properties,  have  been  measured  in  this 
test  device  and  checked  against  predictions  of  a first  generation  composite 
propellant  erosive  burning  model  based  upon  the  bending  of  columnar  diffusion 
flames.  In  general,  the  model  appears  to  give  reasonably  goed  agreement 
with  measured  erosive  burning  data,  except  under  conditions  where  the 
heterogeneity  of  the  composite  propellant  is  unimportant  (low  pressure, 
high  crossflow  velocity) . Here,  it  appears  that  an  additional  mechanism(s) 
of  erosive  burning  will  have  to  be  considered.  The  data  Indicate  that  the 
base  (no  crossflow)  burning  rate  versus  pressure  characteristics  of  the 
propellant  have  a predominant  effect  on  its  sensitivity  to  erosive  burning, 
high  burning  rate  propellants  being  considerably  less  sensitive  to  crossflow 
than  low  burning  rate  formulations,  whether  the  burning  rate  alterations 


are  produced  by  oxidizer  particle  size  variation,  oxidlzer/fuel  ratio  varia- 
tion, or  use  of  catalysts.  Oxidizer  particle  size  appears  to  have  some 
effect  (but  not  a great  one)  beyond  its  effect  on  base  burning  rate,  aug- 
mentation ratio  increasing  with  increasing  particle  size.  Oxidizer/fuel 
ratio  (and  thus  flame  temperature)  appears  to  affect  erosion  sensitivity 
only  through  its  effect  on  base  burning  rate. 

The  first  generation  erosive  burning  model  described  in  this  report 
requires  no-crossflow  burning  rate  versus  pressure  data  as  input:  a second 
generation  model  for  composite  propellants  containing  unimodal  oxidizer 
(still  based  upon  bending  of  columnar  diffusion  flames)  which  does  not  require 
such  data,  but  only  composition  and  ingredient  particle  size  data,  has  been 
developed  and  appears  to  give  reasonable  burning  rate-pressure-crossflow 
velocity  predictions  for  the  one  formulation  against  which  it  has  been  tested. 
Optimization  of  this  model  with  respect  to  selection  of  several  kinetic 
constants  appearing  within  it  has  not  been  completed.  In  addition,  the 
model  has  not  yet  been  extended  to  the  more  useful  case  of  multimodal 
oxidizer  particle  sizes. 

A review  of  the  literature  indicates  that  the  boundary  layer 
profiles  in  rocket  motors  may  differ  significantly  from  those  in  typical 
erosive  burning  test  devices.  Comparison  of  erosive  burning  calculations 
using  the  first  generation  erosive  burning  model  described  in  this  report 
with  profiles  expected  to  prevail  in  the  test  apparatus  versus  those  esti- 
mated to  exist  in  cylindrically  perforated  motor  grains  indicate  that  erosive 
burning  may  be  considerably  less  for  a given  mainstream  crossflow  velocity 
in  such  a motor  than  in  a typical  erosive  burring  testing  apparatus,  in- 
dicating a strong  need  for  further  study  of  boundary  layer  profiles  in  the 
near-wall  region  in  rocket  motor  grain  ports. 
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TTI . MODELING  OF  ALUMINUM  PARTICLE  COMBUSTION  IN  CARBON  DIOXIDE-NITROGEN 
MIXTURE,  ALLOWING  FOR  FINITE-RATE  KINETICS  FOR  THE  ALUMINUM  VAPOR- 
CARBON  DIOXIDE  GAS-PHASE  REACTION 

A.  Introduction  and  Background 

Over  the  past  20  years,  numerous  models  have  been  developed  for 
the  combustion  of  metal  particles  or  droplets,  with  varying  degrees  of 
sophistication  of  the  treatment  of  oxide  product  condensation.^  In 

all  of  tiit. v models,  a flame  sheet  approximation  was  used;  that  is,  metal 
vapor  and  oxidizer  were  assumed  to  react  at  an  infinite  rate  in  an  infinitesi- 
mally thin  zone  removed  some  distance  from  the  particle  surface.  Thus, 
the  physical  processes  of  diffusion  of  the  metal  vapor  and  oxidizer  species 
were  assumed  to  limit  the  rate  of  combustion.  Such  an  assumption,  in  the 
absence  of  detailed  consideration  of  the  complex  temperature  dependency  of 
equilibrium  species  concentrations  in  the  A^O^  system,  results  in  a "d  - 
law"  for  combustion: 

d 2 - d2  - K0  (1) 

o 

d * initial  particle  diameter 
o 

d = particle  diameter  at  time  0 
0 = time  from  beginning  of  combustion 

K = proportionality  constant,  function  of  temperature  and 
gas  composition. 

Macek^  has  experimentally  obtained  burning  time  data  for  single  aluminum 

particles  as  a function  of  their  initial  diameter  which  seem  to  follow  such 

2 (7) 

a d -law.  Pokhil,'  on  the  other  hand,  has  observed  aluminum  combustion  to 
more  closely  follow  a d^'^-law.  (Since  kinetically  controlled  combustion 
of  metal  particles  may  be  shown  to  imply  a d^-law,  Pokhil1  s observations  may 
be  interpreted  as  implying  burning  behavior  which  is  partly  kinetics-controlled 
and  partly  diffusion-controlled:  that  is,  a finite  resistance  due  to  each 
process  must  be  considered.) 

The  question  of  what  particle  size  power  law  is  applicable  to 
aluminum  particle  combustion  is  important  since,  due  to  resolution  limita- 
tions, single  particle  combustion  experiments  have  been  conducted  mostly 
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with  particles  on  the  order  of  50  microns  diameter  or  larger,  while  aluminum 

particles  used  in  solid  propellants  are  typically  on  the  order  of  5 to  10 

microns  in  diameter  (in  the  absence  of  agglomeration).  Assuming  that  a 

50  micron  particle  burns  in  20  milliseconds  under  a given  set  of  conditions, 

1 1.5  2 

use  of  a d , a d or  a d -law  leads  to  predicted  combustion  times  for 

a 5 micron  particle  of  2.0,  0.63  and  0.20  msec,  respectively,  a rather  large 

range  of  times.  Accordingly,  predictions  of  aluminum  particle  combustion 

efficiency  in  typical  low  residence  time  tactical  rocket  motors  are  strongly 

dependent  on  the  assumed  particle-size  combust  ion- time  power  law. 

(8) 

Recently,  Fontijn  measured  kinetic  data  for  the  reactions  of 

aluminum  vapor  with  and  0^  in  a fast-flow  tubular  reactor  experiment. 

The  reaction,  Al(g)  + C02  -*•  A10  + CO,  was  found  to  be  first  order  with 

respect  to  aluminum  vapor  and  first  order  with  respect  to  carbon  dioxide, 

with  a complex  dependency  of  the  rate  constant  on  temperature,  as  shown  in 

Table  III-I.  Using  preliminary  Fontijn  data,  which  indicated  a rate  constant 

-11  3 

of  approximately  10  cm  /molecule  second,  independent  of  temperature,  this 
author  performed  calculations,  described  in  detail  in  Reference  9,  to  eval- 
uate the  validty  of  the  infinite-rate-kinetics  assumption  with  respect  to 
calculation  of  aluminum  particle  combustion  rates.  Two  levels  of  analysis 
were  carried  out,  the  first  being  a very  simple  estimate  of  the  flame  zone 
thickness  and  the  second  involving  development  of  a simplified  droplet 
combustion  model. 

In  the  first  analysis,  a rough  estimate  was  made  of  the  required 
volume  (and  thus  thickness)  of  flame  zone  required  for  quasi-steady-state 
consumption  of  C0^  by  aluminum  vapor  (using  Fontijn' s kinetic  data)  with 
the  C02  being  supplied  at  a diffusion-limited  rate  from  the  surroundings. 

In  this  calculation,  a mass  transfer  Nusselt  Number  of  2 (based  on  the  flame 
diameter)  was  assumed.  In  addition,  it  was  assumed  that  the  average  aluminum 
vapor  concentration  and  CC>2  concentration  in  the  reaction  zone  were  equal 
and  each  was  approximately  ten  to  thirty  percent  of  the  freestream  C02 
concentration.  With  these  assumptions,  the  following  approximate  expression, 
based  on  a mans  balance  on  C02>  was  derived  for  the  flame  zone  thickness: 
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82.05  D T 


(2) 


A * estimated  reaction  zone  thickness  (cm) 

2 

D = diffusivity  (cm  /sec) 

Tn  = ambient  temperature  (°K) 
y_„  ~ ambient  carbon  dioxide  mole  fraction 

*co ,» 

y 

C00 , rx  zone  * "average"  carbon  dioxide  mole  fraction  in  reaction  zone. 

P = pressure  (atm) 

3 

k * rate  constant  for  Al(g)  + C02  A10  + CO  (cm  /mole  sec) 

rg  = particle  radius  (cm) 

6 ■=  average  flame  offset  distance  (cm). 

Results  for  a set  of  calculations  with  pressure  equal  to  one  atmosphere, 
ambient  temperature  of  300°K,  ambient  CO2  mole  fraction  of  0.2,  average 

flame  radius  equal  to  1.10  times  the  particle  radius,  and  diffusivity  equal 

2 

to  2.1  cm  /second  are  presented  in  Table  III-II.  As  may  be  seen,  these 
calculations  indicate  that  at  least  for  particle  sizes  below  20  to  40  microns 
(depending  upon  what  one  accepts  as  a reasonable  value  of  the  ratio  of 
"average"  concentrations  in  the  flame  region  to  free-siream  COj  concentra- 
tion), the  assumption  that  the  reaction  zone  thickness  is  small  compared  to 
the  particle  size  is  poor.  However,  without  further  analysis,  it  is  impossible 
to  define  just  what  a meaningful  "average"  reaction  zone  concentration  is 
or  to  estimate  the  effect  of  a given  reaction  zone  thickness  on  burning  rate. 
Accordingly,  the  question  of  the  effect  of  use  of  finite-rate  kinetics  for 
the  metal  vapor  oxidation  process  on  predictions  of  aluminum  droplet  burning 
rates  was  examined  further  by  use  of  a hydrocarbon  droplet  burning  rate  model 
capable  of  treating  finite  rate  oxidation  kinetics. 

The  treatment  employed  was  quite  similar  to  that  used  by  Peskin 
and  coworkers  ^ for  the  analysis  of  the  combustion  of  hydrocarbon  fuel 
droplets.  The  major  shortcoming  involved  with  using  such  a model  for  treatment 
of  aluminum  combustion  was  that  it  does  not  properly  treat  product  condensa- 
tion effects  which  are  known  to  be  Important  in  the  combustion  of  aluminum. 
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TABLE  III-I.  Rate  Constant  Data  obtained  by  Fontijn  in  a Fast-Flow 

Tubular  Reactor  for  the  Reaction,  Al(g)  + CC^  -*■  A10  + CO. 


Rate  Constant.  k(cc/molecula  second)  Temperature  *K 


1.5 

m 

H 

1 

O 

H 

X-/ 

310 

6.9 

do"13) 

480 

1.6 

✓-x 

CM 

H 

1 

O 

H 

730 

9.0 

do-12) 

1470 

3.8 

(10-11) 

1830 

TABLE  III-II.  Estimated  Reaction  Zone  Thicknesses  for  Aluminum  Particles 
Burning  in  Carbon  Dioxide,  with  Use  of  Fontijn  Rate  Data  for 
Al(g)  + C02  -►  A10  + CO 

Pressure  ■ 1 atmosphere 
Ambient  Temperature  " 30G*K 
Ambient  CO^  Mole  Fraction  «■  0.2 

Reaction  Zone  Thickness  (p) 


Assuming  Assuming 

Y /v  y /Y 

COj,  Rx  tone  CO2,"  CO^,  Rx  zone  CO^,® 

Particle  Radius  (u)  ■ 0.1  ■ 0.3 

50  8 0. 89 

40  10  1.11 

30  13.3  1.48 

20  20  2.22 

10  40  4.45 

5 80  8.90 
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However,  it  was  felt  that  even  though  use  of  such  a model  would  not  give 
accurate  quantitative  burning  rate  predictions,  it  would  permit  examination 
of  the  effect  of  use  of  finite  rate  kinetics  for  the  vapor  phase  oxidation 
step  versus  use  of  infinite  kinetics  (flame  sheet  assumption),  by  comparison 
of  predicted  burning  rates  with  use  of  the  Fontijn  rate  expression  in  the 
derived  equations  versus  burning  rates  calculated  with  use  of  a very  high 
rate  constant  (approaching  infinity)  in  the  seme  equations. 

In  this  preliminary  model  developed  for  study  of  the  effects  of 
finite  kinetics  on  predicted  aluminum  particle  burning  rate,  the  surrounding 
gases  were  assumed  to  be  a mixture  of  CO2  and  N2*  The  overall  oxidation 
reaction,  as  far  as  heat  release  was  concerned,  was  assumed  to  be  3/2  CO2  + 

A1  (g)  - 3/2  CO  r 1/2  A1203  (l),  with  the  reaction  Al(g)  + C02  -*•  A10  + CO 
being  the  limiting  reaction,  all  succeeding  steps  being  fast  in  comparison. 
(This  assumption  of  fast  succeeding  steps,  particularly  condensation,  was 
ohviously  a poor  one  and  must  be  removed  in  any  detailed  analysis.  This 
question  was  studied  in-depth  as  a major  par.t  of  the  aluminum  particle  com- 
bustion modeling  part  of  the  current  program.)  In  addition,  the  product, 
A12°3’  althouRh  bein8  treated  as  a liquid  as  far  as  heat  release  was  con- 
cerned, was  treated  as  a gas  in  terms  of  its  density  and  transport  properties 
in  order  to  permit  use  of  the  Shvab-Zeldovich  equation  formulation.  The 
Lewis  numbers  of  all  species  were  assumed  to  be  unity  (also  necessary  for 
use  of  the  Shvab-Zeldovich  formulation.  The  purpose  of  these  assumptions 
was  to  permit  reduction  of  n + 1 second  order  differential  equations  (where 
n is  the  number  of  chemical  spenies  involved)  to  one  second  order  differential 
equation  plus  n algebraic  equations  (the  Shvab-Zeldovich  formulation). 

Without  the  assumption  that  all  species  behave  in  terms  of  transport  as  gases 
with  a Lewis  Number  of  unity,  the  second  order  differential  equations  were  not 
similar,  and  thus  could  not  be  reduced  to  one  second  order  differential 
equation.  Since  solution  of  these  equations  involves  solution  of  a two- 
boundary eigen-value  problem,  the  reduction  from  n + 1 differential  equations 
to  one  differential  equation  resulted  in  a tremendous  simplification  of  the 
solution  procedure.  In  addition,  for  further  mathematical  simplification, 
properties  such  as  specific  heat,  thermal  conductivity,  and  density-diffusi- 
vlty  product  were  assumed  to  be  Independent  of  position  in  the  analysis. 
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Details  of  the  equation  development  and  solution  procedure  for 
this  simplified  model  are  presented  in  Reference  9.  From  this  analysis, 
particle  burning  rates  were  calculated  as  a function  of  particle  size  for 
a given  ambient  temperature,  pressure,  and  ambient  gas  composition.  These 
calculations  resulted  in  plots  of  burning  rate  versus  particle  size  for 
finite  and  infinite  kinetics.  Numerical  integration  was  then  performed 
to  calculate  the  burning  time  of  a particle  of  given  initial  size  as: 

/rs,o  dr 

s 

(linear  burning  rate) 

for  both  finite  and  infinte  kinetics  cases.  Results  of  a set  of  calculations 
for  P = 1 atmosphere,  Tam^ent  = 2000°K,  and  y^Q^  t ambient  = 0.5  are  presented 
in  Figure  TTI-1  in  the  form  of  burn  time  versus  initial  particle  diameter 
on  a log-log  plot. 

As  may  be  seen,  for  infinite  rate  kinetics  the  particle  combustion 
2 

is  predicted  to  follow  a d -law.  However,  with  the  finite  rate  kinetics 

treatment  using  Fonfijn's  preliminary  data,  the  burn  time  versus  initial 

particle  diameter  predicted  behavior  is  somewhat  more  complex,  with  a shift 

2 

from  approximately  d -law  behavior  for  particles  in  the  lOOp  initial  diameter 
range  to  approximately  d^-law  behavior  in  the  10  to  20p  initial  diameter 
range.  Thus,  the  retarding  effect  of  finite  rate  kinetics  is  quite  important 
for  the  smaller  particle  sizes.  In  fact,  use  of  an  infinite  rate  oxidation 
kinetics'  (flame  sheet)  model  to  extrapolate  lOOp  diameter  data  to  lOp  initial 
particle  diameter  results  in  underprediction  of  the  burn  time  compared  to  that 
calculated  using  Fontijn's  preliminary  data  by  a factor  of  A to  5 for  the 
cases  studied. 

Thus  it  appears  that  the  flame-sheet  approximation  for  the  reaction 
of  metal  vapor  and  oxidizer  used  in  existing  aluminum  combustion  models  is 
a poor  one,  resulting  in  large  errors  in  predicted  burning  times  of  small 
(less  than  20  to  30  micron  diameter)  aluminum  particles,  the  degree  of  error 
increasing  rapidly  as  the  particle  size  is  further  decreased.  Accordingly, 
as  described  below,  an  advanced  model  including  treatment  of  the  kinetics 
of  the  metal  vapor  oxidiation  process  (along  with  the  condensation  processes 
which  was  not  really  treated  in  the  above  preliminary  analysis)  was  developed 
on  this  program. 
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BURN  TIME  (msec) 


Figure  111*1.  Predicted  Burn  Time  u a Function  of  Particle  Size  For  Aluminum 
Particle*,  TAMBIENT  - 2000  K,  Vco2,  AMBI  ENT  “ 0 6 
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B. 


Development  of  Physical  Picture 


As  mentioned  above,  In  the  preliminary  simplified  analysis,  all 
reactions  subsequent  to  the  gas-phase  oxidation  of  aluminum  vapor  via 
Al(g)  + C02  -*■  A10  + CO,  including  condensation  of  were  assumed  to 

be  infinitely  fast.  Among  other  problems  resulting  from  such  an  assumption 
is  the  prediction  of  very  high  temperature  regions  (with  temperatures  well 
in  excess  of  5000°K)  due  to  release  of  heat  in  the  condensation  of  A^O^l). 
At  such  temperatures,  however,  AljO^U)  will  be  endothermically  dissociated 
into  such  species  as  A10,  A102>  &1,0,  etc.,  with  a resultant  lowering  of 
temperature.  That  is,  the  temperature  is  limited  by  the  dissociative  be- 
havoir  of  A^O^Ci),  and  this  feature  must  somehow  be  incorporated  into  any 
physically  realistic  model. 

As  part,  of  the  development  of  a physical  picture  of  processes 
occurring  during  aluminum  particle  combustion  in  carbon  dioxide  atmospheres, 
a procedure  described  by  Henderson^)  was  used  to  identify  likely  reactions 
in  the  gas-phase.  The  key  point  to  this  procedure  is  the  imposition  of  two 
criteria  to  identify  likely  important  reactions.  First,  only  two-body  colli- 
sion reactions  are  considered  because  of  their  much  higher  probability  than 
three-body  collisions.  Secondly,  it  is  required  that  the  absolute  value  of 
the  heat  of  reaction  be  small  in  order  for  a reaction  to  be  likely.  (If 
the  reaction  is  highly  endothermic,  the  fraction  of  sufficiently  energetic 
reactant  molecules  will  be  extremely  low,  while  if  the  reaction  is  highly 
exothermic,  the  energy  release  will  result  in  dissociation  of  the  product 
molecule[s] . ) Under  these  guidelines,  four  possibly  significant  gas-phase 
reactions  are  identified  for  the  aluminum-carbon  dioxide-nitrogen  system: 

^5  kcal/mole  endothermic 


Al(g)  + C02  i A10  + CO 
A10  + A10  *■  AljO  + 0 
A10  + C02  t A102  + CO 

ai2o  + co2  i ai2o2  + co 


^5  kcal/mole  exothermic 
'v 5 kcal/mole  endothermic 
^ kcal/mole  exothermic 


As  seen,  none  of  theae  reactions  yield  Al20^(i):  the  path  to  this  final 
product  is  not  clear.  Both  homogeneous  and  heterogeneous  processes  may 
be  postulated.  For  homogeneous  condensation,  Henderson  has  proposed  a 
series  of  polymerization  reactions  as  follows: 


82 


(A10)2  + (A10)2  -►  (A10)3  + A10 
(Alp)  3 + (A10)2  (A10)a  + AlO 

(A10)4  + (A10)4  - (A10)s  + (AlO) 3 
(A10)5  + (A10)4  -#■  (A10)6  + (AlO)  3 


AH 


3000 


AH 


3000 


AH 


3000 


AH 


3000 


-37  kcal 
-6kcal 
+2  kcal 
-4  kcal 


in  combination  with  formation  of  liquid  Al203  by  evolution  of  aluminum 
vapors  from  the  AlO  polymers: 

(A10)n  -►  (Al2°3)n/3()l)  + n/3  A1 

where  n must  be  large  to  keep  the  heat  of  reaction  small.  This  author  is 
not  aware  of  any  experimental  evidence  supporting  this  scheme.  Moreover, 
it  is  a rather  complex  multistep  scheme  which  as  a consequence  is  probably 
rather  slow.  (Quantitative  estimates  cannot  be  made  since  there  are  no 
rate  data  available  for  any  of  the  reactions  postulated.)  Consequently, 
it  appears  more  likely  that  most  or  all  of  the  condensation  occurs  by  hetero- 
geneous processes  as  discussed  below,  with  reaction  sites  being  supplied  by 
initial  aluminum  oxide  shed  from  the  parent  aluminum  particle,  further 
condensed  phase  alumina  product  trapped  in  a ring  around  the  particle  during 
its  burning,  and  possibly  carbon  particles  resulting  from  partial  breakdown 
of  CO  in  the  fuel-rich  regions  near  the  particle. 

Returning  to  the  four  likely  gas-phase  reactions  listed  above, 
rate  data  obtained  by  Fontijn  for  the  first  reaction  are  available,  as  men- 
tioned earlier.  (See  Table  II1-I.)  For  this  model  development,  Fontljn's 
data  have  been  fit  to  a modified  Arrenhius  expression. 


k - 7.8(10)'13  T1/2  exp  (-1500/T) 


(3) 


3 

k - rate  constant  (cm  /molecule  second) 
T ■ temperature  (#K) 


with  the  rate  being  first  order  with  respect  to  C02  concentration  and  first 
order  with  respect  to  aluminum  vapor  concentration.  Unfortunately,  no  kinetic 
data  are  available  for  the  other  three  reactions.  Accordingly,  it  was  ori- 
ginally Intended  to  examine  two  limiting  cases  with  respect  to  these  reactions: 
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(1)  reactions  negligibly  slow;  and  (2)  reactions  sufficiently  fast  to  be 
treated  as  equilibriun  reactions.  Unfortunately,  time  and  financial  limita- 
tions permitted  examination  of  only  one  of  these  cases.  Comparative  esti- 
mates of  rates  of  all  four  of  the  reactions  based  on  the  procedures  of 
(13) 

Benson  indicate  that  the  second,  third,  and  fourth  reactions  should  be 
two  to  three  orders  of  magnitude  slower  than  that  of  the  Al(g)  + CC>2  -*• 

A10  + CO  reaction.  Accordingly,  it  was  decided  that  the  limiting  case  of 
A10  + A10  **■  A120  + 0,  A10  + C02  -►  A102  + CO,  and  A120  + C02  -*•  A1202  + CO 
all  being  negligibly  slow  be  modeled. 

Next,  we  must  examine  the  question  of  heterogeneous  condensation 
of  AL20^(1).  One  possible  condensation  site  Is  obviously  the  surface  of  Lhe 
parent  aluminum  particle  Itself.  In  addition,  one  may  logically  postulate 
the  existence  of  a fairly  thin  shell  region  some  distance  from  the  particle 
surface  at  which  virtually  all  of  the  A1X0  gas  species  condense  (with  C0? 
supplying  makeup  oxygen  as  needed  to  satisfy  stoichiometry  constraints) 
to  A120^(£).  Provided  that  no  more  than  about  one-third  of  the  aluminui.. 
burning  condenses  at  the  surface  (see  later  discussion) , there  is  a net  ga^ 
flow  away  from  the  surface  of  the  particle  as  it  burns.  However,  at  large 
distances  from  the  particle,  outside  the  region (s)  of  oxide  condensation, 
the  net  gas  flow  is  toward  the  particles.  As  a result,  there  is  some  radius 
at  which  the  net  gas  mass  flux  is  zero.  Any  particles  formed  inside  or 
outside  this  region  will  tend  to  be  convected  to  this  radius.  Thus,  one 
can  easily  visualize  a sink  radius  where  condensed-phase  species  will  tend 
to  gather.  These  species  will  then  serve  as  sites  for  heterogeneous  conden- 
sation at  that  radius.  Now  let  us  consider  the  location  of  this  "condensa- 
tion radius."  If  we  choose  a radius  too  close  to  the  particle,  condensation 
of  aluminum  oxide  there  in  a relatively  small  area  (the  surface  area  associated 
with  the  shell  at  the  condensation  radius  is  proportional  to  the  square  of 
that  radius)  will  lead  to  very  high  temperatures,  above  the  dissociation 
temperature  of  the  aluminum  oxide  and  thus  condensation  will  not  be  able 
to  occur  there.  As  the  chosen  radius  is  increased,  the  temperature  associated 
with  condensation  of  oxide  at  that  radius  decreases.  It  seems  likely  that 
some  supercooling  of  the  vapor  will  be  required  for  condensation  (certainly 
so  if  one  wishes  to  simplify  the  picture  by  considering  a thin  rather  than 
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broad  condensation  zone  for  mathematical  simplification).  Accordingly, 

In  the  model  described  below,  the  "condensation  radius"  was  chosen  as  one 
which  would  result  in  a temperature  at  that  radius,  with  condensation  of  all 
AlxOy  species  there,  of  a specified  value.  (3000°K  and  3500°K  cases  were 
examined.)  Treatment  of  a broad  condensation  zone  rather  than  an  infinites- 
imally thin  one,  with  the  temperature  being  held  constant  across  the  zone 
through  variation  of  the  fraction  of  A1  0 condensed  across  the  zone  would 

X y 

probably  be  more  realistic  than  the  thin  condensation  treatment  employed 
herein;  however,  it  does  not  appear  likely  that  the  calculated  burning  races 
would  be  strongly  changed,  and  the  mathematical  analysis  is  greatly  s'jnplified 
through  use  of  the  thin  condensation  zone  treatment . 

Next,  let  us  examine  possible  processes  occurring  at  the  surface. 
Within  the  constraints  of  limiting  the  gas  phase  chemistry  to  Al(g)  + C02  -+• 
A10  + CO,  there  are  three  possible  major  surface  reaction  processes: 


(1)  Vaporization  of  aluminum  liquid. 

(2)  Reaction  of  C02  and/or  A10  with  liquid  aluminum  at  the  surface 
to  yield  Al20(g). 

(3)  Reactive  condensation  of  A120  and/or  AlO  at  the  surface  to 
yield  AljO^l)  which  then  remains  on  the  surface. 

The  first  process,  vaporization  of  liquid  aluminum,  is  obviously  of  importance 
and  must  be  considered  in  any  model  of  aluminum  droplet  combustion.  Heat 
.nust  be  supplied  either  from  heat  feedback  from  the  outer  zones  (most  impor- 
tantly the  outer  condensation  zone)  or  from  the  exothermic  reactions  in 
Categories  2 and  3. 

The  reaction  of  AlO  with  Al(t)  at  the  surface  via  the  reaction: 


AlO  + Al(4)  -►  Al20(g) 

not  only  aids  in  the  consumption  of  the  aluminum  particle  through  direct 
removal  of  liquid  aluminum  from  the  surface  but,  being  exothermic  by  approxi- 
mately 55  kcal/mole,  also  releases  heat  for  vaporization  of  additional 
liquid  aluminum.  The  reaction  of  C02  with  Al(t)  at  the  surface  via: 

C02  + 2A1(£)  ■*>  CO  + Al20(g) 

directly  removes  two  moles  of  aluminum  from  the  particle  per  mole  of  C02 
reacting  there,  but  is  slightly  endothermic,  thus  requiring  heat  from  the 


85 


A10-A1  reaction  or  from  the  outer  regions.  Two  cases  were  examined  in  the 
course  of  this  modeling  effort:  in  the  first,  both  the  A10-A1  and  CO^-Al 
surface  reactions  were  assumed  to  occur  with  infinite  Kinetics;  while  in  the 
second,  the  CC^-Al  surface  reaction  was  prohibited.  It  is  noL  clear  which 
of  these  scenarios  is  most  .ealistic. 

The  question  of  condensation  of  aluminum  oxide  at  the  surface 
appears  to  be  quite  a controversial  one.  Preliminary  thermodynamic  analysis 
with  estimated  values  of  surface  temperature  indicated  that  thermodynamically 
AljO^  should  not  be  formed  at  the  surface,  the  reaction  3Al20(g)  t A^O^fJl)  + 
4A1(£)  being  shifted  all  the  way  to  the  left  for  reasonable  values  of  A^O 
partial  pressure  adjacent  to  the  surface.  However,  there  are  considerable 
references  in  the  literature  to  the  existence  of  A^O^  on  burning  aluminum 
particle  surfaces.  While  in  most  of  these  references  it  is  not  clear  whether 
the  A^O^  on  the  surface  comes  from  initial  A^O^  formed  in  preignition  and 
ignition  phases  or  whether  it  is  actually  built  up  during  the  combustion 
process,  several  (e.g..  References  14  - 16)  clearly  state  that  oxide  accumu- 
lates on  the  surface  during  the  self-sustained  combustion.  If.  in  interesting, 
however,  to  note  that  no  such  oxide  formation  on  the  surface  is  observed  in 

In 


(17) 


combustion  of  particles  in  0- /Argon  and  oxygen-rich  CL /CO  flames. 

(18)  ^ 
line  with  this,  Prentice'  suggest'-  '.hat  nitrogen  and  carbon  appear  to 

be  important  factors  in  causing  accumulation  of  oxide  (and  nitride  ?)  on 

the  surface.  It  should  be  pointed  out  that  these  observations  and  conclusions 

regarding  aluminum  oxide  accumulation  on  burning  aluminum  particles  are, 

in  general,  based  on  experiments  with  rather  large  (greater  than  100  micron 

diameter)  aluminum  particles.  A particularly  interesting  conclusion  regarding 

the  dependency  of  oxide  accumulation  on  initial  particle  size  is  presented 

in  Reference  19.  In  this  report,  a plot  of  weight  fraction  of  initial  aluminum 

appearing  in  the  final  residue  globule  after  burning  versus  original  particle 

diameter  is  presented  for  the  case  of  aluminum  particles  burning  in  1 atmosphere 

50/50  O2/N2  environment.  This  plot  indicates  that  this  weight  fraction 

decreases  with  decreasing  Initial  particle  diameter.  Extrapolation  of  the 

data  indicates  that  for  initial  particle  diameters  below  about  100  microns, 

no  accumulation  will  occur.  Since  the  model  development  on  this  program  was 

aimed  at  small  particles  (60  micron  diameter  or  less)  where  kinetic  effects 
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are  more  likely  to  be  important  relative  to  diffusion  effects,  it  was  decided 
that  oxide  condensation  on  the  burning  particle  surface  be  neglected  in  this 
modeling.  Outputs  of  the  model(s)  discussed  below  include  such  parameters 
as  surface  temperature  and  gas  phase  composition  adjacent  to  the  particle 
surface.  With  these,  the  validity  (in  terms  of  thermodynamics)  of  the  assump- 
tion of  no  condensed  oxide  formation  at  the  surface  can  be  checked.  Interest- 
ingly enough,  the  assumption  was  found  to  break  down  only  for  the  larger 
particles,  in  line  with  the  trend  observed  by  Prentice  of  fractional  surface 
oxide  accumulation  increasing  with  increasing  initial  particle  size. 

The  final  physical  picture  modeled  is  sketched  in  Figure  III-2. 
Actually,  two  models  were  developed,  differing  only  in  the  assumption  of 
infinitely  fast  or  infinitely  slow  reaction  kinetics  for  the  reaction  of  CO^ 
with  liquid  aluminum  at  the  particle  surface  to  form  A^O  and  CO.  Each  model 
is  divided  into  two  regions,  an  inner  region  and  an  outer  region,  separated 
by  a condensation  zone.  There  are  no  reactions  in  the  outer  region,  where  the 
only  species  considered  are  CO,  CO2 , and  N2.  CO  and  heat  are  transferred 
outward  from  the  condensation  zone  in  this  region  while  CO2  is  transferred 
inward.  At  the  surface,  aluminum  vapor  is  boiled  off,  using  heat  supplied 
by  feedback  from  the  outer  regions  and  by  the  infinitely  fast  exothermic 


(A10)  + Al(£)  AI2O  reaction.  In  addition,  aluminum  is  removed  from  the  surface 
as  AI2O  by  this  reaction  and,  in  one  of  the  model  versions,  by  the  reaction 
CO2  + 2A1(£.)  •+•  A^O  + CO.  In  the  inner  region,  between  the  particle  surface 
and  the  condensation  radius,  aluminum  vapor  reacts  with  CO2  at  a rate 
determined  by  the  previously  described  rate  expression  derived  from  the 
data  of  Fontijn.  Convection  and  diffusion  of  all  species  and  conduction 


and  convection  of  sensible  enthalpy  are  treated.  All  Alx0^  species  trans- 
ported to  the  condensation  radius  react  infinitely  rapidly  there  with  CO2 
supplied  from  the  outer  region  to  form  which  is  then  trapped  there. 

Heat  released  from  these  condensation  reactions  is  transferred  in  either 


direction  away  from  this  zone.  The  condensation  radius  is  chosen  by  a trial 
and  error  procedure  to  yield  a specified  temperature  at  that  radius.  (The 
further  out  the  condensation  radius,  the  lower  the  temperature.)  Two  values 
(3000°K  and  3500°K)  of  condensation  temperature,  somewhat  below  the  atmos- 
pheric decomposition  temperature  of  AljO^  (approximately  3700  to  3900°K, 
depending  rn  composition)  were  chosen  for  parametric  study. 
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Figure  111-2.  Schematic  of  Single  Particle  Aluminum  Combustion  Model(s)  Developed  in  This  Study. 


c. 


Model  Development 


As  mentioned  above,  two  models  of  single  particle  aluminum  combustion 
were  developed  during  the  course  of  this  program,  the  only  difference  being 
allowance  or  prohibition  of  an  infinitely  fast  reaction  of  C02  with  aluminum 
liquid  at  the  burning  particle  surface.  As  a result,  the  conservation  equa- 
tions in  the  inner  and  outer  regions  are  the  same  in  each  model  as  are  the 
condensation  radius  boundary  conditions  used  to  connect  the  inner  and 
outer  regions.  The  only  difference  occurs  in  the  boundary  conditions 
at  the  particle  surface. 

In  the  outer  rtgion  (r  > r , the  condensation  radius)  there  are 

c 

only  three  chemical  species,  C02>  CO,  and  N2<  Since  the  problem  is  formulated 
as  a quasi-steady-state  one,  it  may  easily  be  shown  that  the  overall  molar 
flux  in  this  region  is  zero  since  there  must  be  no  net  flux  of  carbon  or 
nitrogen  atoms  to  r = r^  and  since  CO  and  C02  each  contain  one  carbon  atom. 

In  addition,  there  are  no  reactions  allowed  in  the  outer  regions.  Accordingly, 
the  species  conservation  equations  for  CO,  C02  and  N2  in  the  outer  region 
reduce  simply  to: 


(4) 


where : 


r 

c 


D 


radius 

total  molar  concentration 
mole  fraction  of  species  j 
dif fusivity. 


Substituting  y 


p/MW  (where  is  the  partial  pressure  of 


Pj/P  and  C 

species  j , p is  the  overall  gas  density,  and  MW  is  the  overall  gas  molecular 
weight)  and  assuming  that  the  density-diffusivity  product  is  essentially 
independent  of  position,  we  arrive  at,  following  integration  of  Equation  4: 


2 9P 

— E 1 - k (5) 

P(MW)  3r  j K , 

where  is  a constant  to  be  evaluated  from  application  of  boundary  conditions 
and  is  different  for  each  species.  (Since  N2  is  the  only  nitrogen-containing 
species  and  the  net  molar  flux  of  N2  must  be  zero  for  quasi-steady-state 
conditions  to  apply,  it  may  easily  be  shown  that  k^  ■ 0 for  N2«)  The  average 
molecular  weight  varies  with  composition  and  may  be  expressed  as: 
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MW 


(6) 


44  + 28(P  - P^) 


Substitution  of  this  expression  into  Equation  5,  integration  of  that  equation 
and  application  of  the  boundary  conditions  that  the  various  partial  pressures 
at  r + • are  equal  to  ambient  values  (0.0  in  the  case  of  CO)  leads  to: 


PC02-  <PC02,  A„b  + 1-75P)<!'16  Vr-1.75P 


(7) 


?nn  - 2.75?  - P„  . , - (P„  . . + 1.75  P)e"16  k2/r  (8) 

CO  N2.  Amb  C02,  Amb  v 


P “ P 
N2  N2,  Amb 


(9) 


for  r > rc,  where  k2  must  be  evaluated  through  boundary  condition  matching 

at  r ■ r . 

c 

The  energy  conservation  equation  in  the  outer  region  (r  > r ) 

c 

may  be  written  as: 


1_  _9_ 
2 3r 


(10) 


T ■ temperature 
A • thermal  conductivity 
C^  ■ specific  heat 

u ■ linear  velocity  relative  to  the  particle  surface 


where  the  first  term  is  a conduction  term  and  the  second  term  is  a convection 

term.  For  integration  of  this  equation  it  is  necessary  that  an  expression 

for  the  mass  flux,  pu,  as  a function  of  radius  be  substituted,  with  this 

expression  preferably  relating  pu  to  the  total  mass  flux  of  aluminum  leaving 

the  particle  surface  (m  „„) . Since  all  of  the  aluminum  leaving  the  surface 

p.TOT 

is  assumed  to  disappear  as  A120^(£)  into  a single  sink  at  r - r^,  it  may 
easily  be  shown  from  a mass  balance  and  stoichiometry  considerations  that 
the  net  mass  flux,  pu,  at  r ■ r£+  is  given  by: 


2 

rc  pu 


• A 

" ““p.TOT  p \ 


Mass  of  oxygen  per  mole  of  A120 
Mass  of  aluminum  per  mole  of  Al2 
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-0.899  m Tnxr 
p,TOT  p 
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Then,  since  there  are  no  other  mass  sinks,  the  mass  flux  at  any  other 

radius  greater  than  r - r is  related  to  that  at  r ■ r + by  the  inverse 

c c J 

square  of  the  radius: 


pu  - 


0.889  m Tn_r2 
p .TOT  p 

2 

r 


(11a) 


• 

m “ mass  flux  away  from  particle  surface 

p , TOT 

r * particle  radius. 

P 

Thus,  Equation  10  becomes: 
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3r 
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(12) 


With  the  assumption  that  the  thermal  conductivity.  A,  is  independent  of  posi- 
tion (analogous  to  the  earlier  assumption  that  the  density-dif fusivity  product 
is  independent  of  position),  this  equation  may  be  integrated  to  yield: 
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0.889m  r C /A 
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TOT  p p 


(13) 


As  mentioned  earlier,  it  is  assumed  that  the  condensation  radius, 

r * rc>  is  a location  at  which  all  aluminum  containing  species  arriving 

there  react  infinitely  rapidly  to  form  condensed  A^O^,  any  needed  oxygen 

being  supplied  by  diffusing  and/or  convecting  in  from  the  outer  region. 

Thus,  one  may  immediately  write  as  part  of  the  boundary  conditions  to  be 

satisfied  at  r ■ r , 
c 

*Al|  +“°  (14) 

1 r ■ r 

c 


ai2o 


(15) 

(16) 


In  addition,  applications  of  conservation  equations  for  enthalpy,  oxygen 
atoms,  carbon  atoms,  and  nitrogen  atoms  at  r - r£  yield,  when  combined  with 
equations  7-9  and  13: 
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where  Q^q  is  the  heat  release  associated  with  Al(g)  + 3/2  C02  -*-1/2  A^O^fc) 
+ 3/2  CO,  Q is  the  heat  release  associated  with  AlO  + 1/2  CO^  1/2 
A1203(£)  + 1/2  CO,  and  Q^2  ifl  the  heat  release  associated  with  Al20(g)  + 

2C02  -*•  A1203(£)  + 2CO. 

In  the  inner  region  between  the  particle  radius  (r  ■ r ) and  the 

P 

condensation  radius  (r  ■ r ),  the  species  to  be  considered  are  (Figure  III-2) 

N- , CO,  CO,,  Al,  AlO,  and  A1,0.  Applying  mass,  enthalpy  and  species  conser- 

1 L (20) 

vation  equations,  using  the  procedures  of  Williams,  we  find: 


1 3 l 2,  3T»  _ 3T  , „ . 

7 a?  {r  x37)'  p'uCP  3F  + Q1R1  “ ° 


(21) 


1312  3<P</MW> 


3(P/MW) 


- pu- 
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- PRj  Vj 


(22a-f ) 


j « N , CO,  C02,  Al,  AlO,  A120 


2 2 
» TflTr  “pur 
p.tot  p 


(23) 


where  the  v^'s  are  the  stoichiometric  coefficients  in  the  reaction  Al(g)  + 
C02  -*•  AlO  + CO  (vN2  - 0.0,  vA1  - 1.0,  vCo2  " i-0,  vA10  “ 1.0,  vCQ  * -1.0, 
Uai2o  “ 0.0),  is  the  rate  of  the  reaction  (moles/volume/time),  and 
is  the  heat  release  associated  with  the  reaction.  Substitution  of  Equa- 
tion 23  into  Equations  21  and  22  yields  seven  second-order  differential 
equations  to  be  solved  over  the  inner  region  subject  to  boundary  conditions 
at  the  condensation  radius  given  above  (Equations  14  - 20)  and  the  boundary 
conditions  at  the  particle  surface  which  are  discussed  below. 
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As  mentioned  earlier,  two  slightly  differing  models  were  examined 
during  the  course  of  this  program,  with  infinitely  fast  reaction  of  CC^ 
with  liquid  aluminum  at  the  particle  surface  being  allowed  in  the  first  and 
prohibited  in  the  second.  Let  us  first  examine  the  inner  boundary  conditions 
for  the  case  where  C02  is  allowed  to  react  at  the  surface.  In  this  case, 
we  can  immediately  write  as  two  boundary  conditions: 


AlO.r  - r “ ° 

P 

(26) 

C0_,  r - r = 0 

2 P 

(27) 

In  addition,  conservation  equations  for  the  other  four  species  permit  deri- 
vation of  four  additional  boundary  conditions: 
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A seventh  boundary  condition  is  obtained  from  an  enthalpy  balance  at  the 
particle  surface  as: 


mevapLevap  ,3T  , D(MW) 
27  " 3r  RT 


, ^PA!0/MW) 
5 3r 


3(pco2/mw) 

+ V ‘ ' 3r  ' 


(32) 


r - r 


where  m is  the  evaporation  mass  flux  of  aluminum  from  the  particle  surface 
evap 

L is  the  heat  of  vaporization  of  aluminum  and  Qr  and  Q_  are  heat  releases 
evap  ^5  v 


(31) 
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associated  with  Al(£)  + A10  -►  Al^O  and  2A1(£)  + CC^  AljO  + CO,  respectively. 
Since  aluminum  leaves  the  surface  by  both  evaporative  and  reactive  processes, 
one  more  equation  relating  the  total  mass  flux  to  the  evaporative  flux  is 
given  by: 
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p ,TOT  evap 


27D(MW) 

RT 


9<pa,o/mw) 

9r 


• 0*co  /MW> 


+ 2- 


9r 


(33) 


rp 


Finally,  t>  close  the  problem,  a rate  expression  using  Fontijn's  kinetic 
data  (Equation  3)  is  substituted  for  in  Equations  24  and  25,  and  an 
expression  relating  the  partial  pressure  of  aluminum  vapor  at  the  particle 
surface  to  the  surface  temperature  is  employed.  This  latter  expression 
was  obtained  by  curve-fitting  JANNAF  thermodynamic  data,  yielding: 


T - 34860/(12.537  - In  PA,  ) (34) 

r = r Al,r  **  rp 

P 


In  the  second  model,  where  CO2  Is  not  allowed  to  react  with  aluminum 
liquid  at  the  surface,  the  only  changes  in  the  particle  surface  (inner) 
boundary  conditions  are  the  substitution  of: 


D(MW) 


mp,T0TRTs 

P(MW) 


(35) 


for  equation  27,  and  the  elimination  of  the  9 (P^q^/MW) /3r  term  in  Equations 
29,  31,  32,  and  33. 

The  problem  at  this  point  is  to  find  solutions  to  seven  coupled 
non-linear  second-order  differential  equations  (24,  25a-f)  which  will  satisfy 
the  seven  outer  boundary  condition  relationships  given  by  Equations  14  - 20 
and  the  seven  inner  boundary  conditions  relationships  and  ancillary  equations 
given  by  Equations  26  - 34  (or,  for  the  second  model.  Equations  26,  28,  30, 

34,  35,  and  modified  versions  of  29,  31,  32,  and  33).  In  order  to  simplify 
this  problem,  two  major  approximations  were  made.  First,  the  Lewis  Numbers 
(ratio  of  thermal  diffusivity  to  species  diffusivity)  were  assumed  to  be  unity 
for  all  species,  permitting  use  of  the  Shvab-Zeldovich  procedure  to  reduce 
the  seven  second-order  differential  equations  to  one  second-order  differen- 
tial equation  plus  six  algebraic  equations  relating  the  various  species 
concentrations  to  one  another  and  to  temperature  at  any  given  position. 
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In  addition,  the  thermal  conductivity  (and  thus  the  density-dif fusivity 
product)  were  assumed  to  be  independent  of  position. 

In  the  employment  of  the  Shvab-Zeldovich  procedure,  the  following 
new  variables  were  defined: 
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■ <PA10  + PA!  + 2PA120>'«M 

(36) 
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and  substituted  into  Equations  24  and  24a-f  to  yield: 
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where : 


4.70(10U)T1/2e"150°/TC  C 

A1  CO. 


6.98(107)  P^fco^'150071/!1'5 


(44) 


With  the  approximation  that  the  thermal  conductivity,  A,  is  independent  of 
position,  Equations  43a-f  can  then  be  integrated  to  yield  the  closed-form 
expressions : 
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(48) 
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+ d6 
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* 2 

where  K * m r C and  c and  d.  are  constants  of  integration  to  be 

1 p,TOT  P P j j 

evaluated  using  the  boundary  condition  equations. 

Equations  36  - 41  are  also  substituted  into  the  appropriate  boundary 

condition  equations  to  convert  them  into  the  new  variables.  Values  for  the 

condensation  radius  (r£)  and  the  total  mass  flux  at  the  particle  surface 

(m  ) are  then  estimated  for  given  ambient  conditions  and  particle  radius 

p ,TOT  6 

and  the  inner  and  outer  boundary  conditions,  ancillary  equations,  and  Equations 

45  - 50  are  solved  simultaneously  (a  miserable  algebra  problem)  to  yield 

first-cut  values  for  c thru  c , d thru  d , temperature  at  the  condensation 

1 6 1 6 

radius,  and  temperrture,  composition  and  temperature  gradient  at  the  particle 

surface.  Successive  estimates  are  made  on  the  condensation  radius  (at  a 

fixed  estimate  of  m _._)  until  the  desired  condensation  radius  temperature 

p , TOT 

is  calculated.  At  this  point,  the  calculated  surface  temperature,  composition 
and  surface  temperature  gradient  are  used  to  start  numerical  integration  of 

Equation  42  from  the  particle  surface.  Equations  40,  41,  49  and  50  are  used 

to  relate  aluminum  vapor  and  carbon  dioxide  partial  pressures  to  temperature 

at  each  point  to  permit  evaluation  of  from  Equation  44  at  each  step  in 

the  integration.  Equation  42  is  integrated  from  r ■ r^  to  r « where  the 

calculated  parameter  values  are  tested  to  see  whether  they  satisfy  the  outer 

boundary  conditions.  If  not,  a new  value  of  m is  estimated, 

p.TOT 

and  the  entire  procedure  (including  searching  for  the  condensation  radius 
which  yields  the  desired  condensation  temperature)  is  repeated. 
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Application  of  this  procedure  to  the  second  model  (in  which  the 

reaction  of  CO^  with  aluminum  liquid  at  the  particle  surface  was  prohibited) 

was  relatively  straightforward  (discounting  the  horrendous  algebra).  However, 

when  the  CO^  surface  reaction  was  permitted  (with  infinite  kinetics)  an 

interesting  (and  upon  reflection,  to  be  expected)  result  occurred;  namely, 

the  temperature  gradient  at  the  surface  was  calculated  to  be  infinity  (from 

a combination  of  Equations  32  and  33  with  m eliminated)  at  all  estimated 

evap 

values  of  ra^  except  one,  at  which  it  was  calculated  to  be  indeterminate 
(zero  over  zero).  Examination  of  the  problem  revealed  that  the  value  of 

m which  gave  an  indeterminate  value  of  this  gradient  was  the  correct 

answer  and  was  independent  of  the  kinetics  of  the  gas-phase  Al(g)  + CC^ 

A10  + CO  reaction.  Reflection  indicated  that  this  was  logical,  given  the 
assumptions  made  in  this  model.  Since  it  was  assumed  that  all  CO^  reactant 
which  reaches  the  surface  without  reacting  with  aluminum  vapor  reacts  infin- 
itely fast  there  to  produce  the  same  product  as  any  A10  from  the  gas-phase 
reaction  which  reaches  the  surface  (namely  Al^O)  and  since  it  was  furthermore 
assumed  that  all  Lewis  Numbers  were  unity,  the  total  mac u rin  from  the  surface 
is  naturally  independent  of  the  rate  of  the  gas-phase  reaction.  While  the 
details  of  the  temperature  and  concentration  profiles  in  the  gas-  phase  depend 
on  the  gas-phase  kinetics,  the  change  in  sensible  enthalpy  feedback;  from  the 
gas-phase  due  to  the  gas-phase  reaction  is  lust  offset  by  reactive  oxidizer 
species  gradient  changes  resulting  from  the  gas-phase  reaction.  Thus,  we 
arrive  at  the  very  interesting  result  that  if  we  allow  all  products  of  any 
gas-phase  reactions  (including  not  only  Reaction  1,  but  Reactions  2,  3,  and 
4 as  listed  earlier)  and  all  unreacted  oxidizer  (CC^)  to  react  infinitely 
fast  at  the  surface  to  the  same  product  (and  if  the  Lewis  Number  ■ 1 assumption 
is  reasonably  good)  the  burning  rate  of  an  aluminum  particle  is  independent 
of  gas-phase  kinetics  as  implicitly  assumed  in  past  models.  If,  however, 
oxidizer  which  does  not  react  in  the  gas-phase  cannot  react  at  the  surface 
while  products  (e.g.,  A10)  of  gas-phase  reactions  can,  the  burning  rate  will 
depend  on  the  gas-phase  kinetics  as  shown  by  the  results  presented  in  the 
next  section.  As  would  be  expected,  if  the  gas-phase  kinetic  rates  are 
assumed  to  be  infinite,  results  of  the  second  model  collapse  to  those  of 
the  first. 
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Further  discussion  will  be  mainly  limited  to  the  second  model, 
in  which  CC^  is  prohibited  from  reacting  with  liquid  aluminum  at  the  surface. 

A computer  program  was  written  to  solve  the  equations  presented  earlier 
using  the  procedure  outlined  above,  giving  predicted  surface  mass  fluxes 
as  a function  of  ambient  conditions,  particle  size,  and  assumed  condensation 
temperature.  These  surface  mass  fluxes  were  then  straightforwardly  converted 
to  linear  burning  rates.  Burning  times  for  particles  of  various  initial 
sizes  were  then  calculated  by  integrating  burning  rate  versus  particle  size 
predictions  (for  given  ambient  conditions  and  assumed  condensation  temperatures): 


T 


burn 


initial 

dr 

(Linear  Burning  Rate) 


o 

Results  are  presented  and  discussed  in  the  next  section. 


(51) 


D.  Results 

As  indicated  above,  only  results  from  the  second  model  (COj  surface 
reactions  prohibited)  will  be  discussed  here,  since  the  first  model  (infin- 
itely fast  reaction  of  CC^  with  liquid  aluminum  at  particle  surface  allowed) 
collapses  to  a diffusion-limited  result,  particle  burning  rate  being  inde- 
pendent of  gas-phase  kinetics.  That  is,  no  matter  what  rate  expression  is 
used  in  the  first  model  for  the  Al(g)  + CO^  •+  A10  + CO  reaction,  the  predicted 
burning  rate  is  the  same  as  that  which  would  be  predicted  with  the  second 
model  using  an  infinite  rate  for  that  reaction. 

Predicted  temperature  and  species  concentration  profiles  around 
a burning  aluminum  particle  are  presented  for  three  typical  cases  in  Figures 
III-3  through  III-5.  In  all  of  these  cases,  the  temperature  at  the  conden- 
sation radius  was  chosen  to  be  3500°K.  (As  mentioned  earlier,  we  specify 
the  temperature  at  the  condensation  shell  and  then  solve  for  its  location.) 
Since  the  condensation  shell  is  postulated  to  be  an  infinitesimally  thin 
region  in  which  considerable  reaction  occurs,  there  are,  as  would  be  expected, 
discontinuities  in  the  temperature  derivative  and  the  CO  and  CO2  partial  pres- 
sure derivatives  at  this  radius.  For  the  case  depicted  in  Figure  III-3 
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Figure  111-4.  Predicted  Temperature  and  Concentration  Profiles  Around  Burning  Aluminum  Particle. 

Case  17,  rp=  TAMB  = 1000°K,  P=  1 atm,  Yco  AMB  * 0.5,  Trc  = 3500,  m/m00=  0.53. 


(Particle  radius  ■ 15  microns,  ambient  temperature  - 1000#K,  pressure  ■ 

1 atmosphere,  ambient  CO2  mole  fraction  ■ 0.5),  the  diffusion  processes 

are  slow  enough  relative  to  the  kinetics  of  the  AI-CO2  gas  reaction  that 

this  reaction  occurs  to  an  appreciable  extent  as  may  be  seen  from  the  A10 

and  Al^O  partial  pressure  curves.  (Under  the  assumptions  of  this  model, 

A^O  can  be  generated  at  the  surface  only  if  A10  is  produced  by  the  initial 

Al-CO^  gas  reaction,  since  A10  is  the  only  species  permitted  to  react  at 

the  surface  to  form  A^O.)  For  this  particular  case,  the  burning  rate  of 

the  aluminum  particle  is  predicted  to  be  74%  of  what  it  would  be  if  the 

Fontiln  kinetics  were  replaced  with  an  infinite  rate  constant  (m_  . , / 

Fontijn 

m.  ..  . . « 0.74).  The  fact  that  the  A1-C0.  gas  reaction  occurs  to  an 

infinite  kin  2 

appreciable  extent  in  this  case  is  also  shown  in  the  CO  partial  pressure  profile, 
which  shows  a peak  value  between  the  particle  radius  and  the  condensation 
shell  due  to  its  production  by  the  AI-CO2  gas  reaction  in  this  inner  zone. 

In  contrast.  Figures  III-4  and  III-5  depict  predicted  profiles  for 
cases  in  which  the  amount  of  AI-CO2  gas-phase  reaction  is  very  small.  In 
both  these  cases,  the  particle  radius  is  considerably  smaller  (1  micron) 
than  that  of  the  Figure  III-3  case,  leading  to  a much  higher  ratio  of  charac- 
teristic reaction  time  to  diffusion  time  (lower  value  of  the  Damkohler  Number). 

As  may  be  seen,  considerably  less  AlO  and  A^O  are  produced  in  these  cases, 
with  the  evaporation  of  aluminum  thus  basically  being  driven  by  heat  feedback 
from  the  condensation  zone,  most  of  this  heat  being  produced  by  the  hetero- 
geneous reaction  of  Al(g)  and  CC^  at  particle  sites  in  this  thin  zone  to 
A^O^*.).  In  these  cases,  the  CO  concentration  peaks  at  the  condensation 
radius  rather  than  in  the  region  between  the  particle  surface  and  the  conden- 
sation radius.  In  addition,  the  ratio  of  predicted  burning  rate  to  that 
which  would  be  predicted  with  infinite  kinetics  for  the  Al-CO^  gas-phase 
reaction  drops  to  approximately  0.5. 

Burning  rate  predictions  have  been  carried  out  using  the  model 
described  above  with  Fontijn's  reaction  rate  data  for  the  Al(g)  + CO2  ■+• 

AlO  + CO  reaction  for  a number  of  cases  covering  a wide  range  of  particle 
sizes  and  ambient  conditions.  Results  of  these  calculations  are  summarized 
in  Table  III-III  and  Figures  III-6  and  III-7.  In  Table  III-III,  the  inde- 
pendent variables  tabulated  in  Columns  2 to  5 are  particle  radius,  ambient 
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PARTICLE  RADIUS  (micron*) 


Figure  1 1 1-7.  Predicted  Ratio  of  Burning  Rate  with  Fontijn  Kinetics  to 
That  with  Infinite  Kinetics. 


temperature,  pressure,  and  CO^  ambient  partial  pressure  (with  the  remaining 
ambient  gas  being  nitrogen).  Two  values  of  condensation  radius  temperature, 
3000°K  and  3500°K,  were  studied.  The  predicted  mass  burning  rate  (mass  per 
unit  area  per  time)  was  found  to  be  very  nearly  independent  of  this  tempera- 
ture over  this  range  for  the  case  of  infinite  kinetics  for  the  Al(g)  + 

A10  + CO  reaction  (or  for  the  case  where  is  allowed  to  react  infinitely 
rapidly  with  liquid  aluminum  at  the  particle  surface,  as  discussed  earlier). 
This  predicted  mass  burning  rate  for  infinite  kinetics  is  presented  in  Column 
6.  With  finite  kinetics  (and  CO^  reaction  with  liquid  aluminum  prohibited), 
however,  the  predicted  burning  rate  was  found  to  depend  on  condensation  radius 
temperature,  decreasing  with  decreased  temperature  (logically  enough).  In 
columns  7 through  10,  the  mass  burning  rate,  ratio  of  that  burning  rate  to 
the  infinite-kinetics  rate,  linear  burning  rate,  and  condensation  radius 
location,  are  tabulated  for  the  3500°K  condensation  temperature  case,  while 
the  same  outputs  for  the  3000°K  case  are  presented  in  columns  11  through  14. 

As  may  be  seen,  the  condensation  shell  location  is  not  particularly  dependent 
on  the  condensation  temperature.  (If  the  burning  rates  were  the  same,  the 
condensation  shell  location  would  be  further  out  for  the  lower  temperature 
case,  but  as  it  turns  out  the  decrease  in  burning  rate  associated  with  the 
decreased  condensation  temperatures  results  in  a shift  of  the  condensation 
shall  back  inwards.)  In  some  cases,  the  condensation  shell  radius  is  slightly 
less  for  the  3000°K  condensation  temperature;  in  other  cases  it  is  slightly 
greater. 

In  Figures  III-6  and  III-7,  the  ratio  of  burning  rate  predicted  using 
Fontijn  kinetics  to  that  predicted  using  infinite  kinetics  for  the  aluminum 
vapor  reaction  with  CC^  is  plotted  against  particle  radius  for  various  sets 
of  ambient  conditions  (temperature,  pressure,  CC^  tuole  fraction).  The 
ratios  in  Figure  III-6  a*.e  calculated  for  a condensation  shell  temperature 
of  3500°K,  while  those  in  Figure  III— 7 are  for  3000°K.  As  may  be  seen,  the 
ratios  (and  thus  the  absolute  rates,  since  the  infinite-kinetics  burning  rate 
is  essentially  independent  of  condensation  shell  temperature  over  this  range) 
are  lower  for  the  3000°K  case.  As  particle  radius  becomes  larger,  diffusion 
becomes  slower  relative  to  kinetics,  and  the  burning  rate  ratio  accordingly 
approaches  unity  at  very  large  particle  radii.  Logically  enough,  since  the 
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aluminum  vapor  reaction  with  CC^  is  second  order  with  respect  to  pressure, 
this  approach  of  the  burning  rate  ratio  to  unity  occurs  at  lower  values  of 
particle  radius  for  higher  pressures.  (Compare  the  1,  5 and  10  atmosphere 
curves  for  Tamjj^ent  * 1000°K,  Y^Q^ambient  “ 0.5.)  At  very  small  values  of 
particle  radius,  on  the  other  hand,  where  the  characteristic  diffusion  times 
are  small  compared  to  the  characteristic  reaction  times,  one  might  expect 
the  ratio  of  burning  rate  calculated  using  Fontijn  kinetics  to  that  calculated 
using  infinite  kinetics  to  drop  to  very  low  (approaching  zero  as  particle 
radius  approaches  zero)  values.  This  does  not  occur,  however,  since  as  the 
rate  of  the  gas-phase  reaction  drops  well  below  the  diffusion  rate,  the  loc- 
ation of  the  condensation  zone  begins  to  dominate  the  burning  rate.  (That 
is,  the  heat  feedback  from  the  condensation  zone  begins  to  dominate  the  rate 
of  removal  of  aluminum  from  the  particle  surface,  with  removal  by  the  A10(g)  + 
Al(£)  ■+■  A^^Hg)  reaction  becoming  comparatively  small.)  Examination  of  Table 
III-III  reveals  that  as  particle  radius  becomes  very  small  (less  than  about 
5 microns  for  the  1 atmosphere  cases)  the  ratio  of  condensation  radius  to 
particle  radius,  which,  in  general,  decreases  with  decreasing  particle  radius, 
approaches  a constant  value  (as  is  the  case  for  all  particle  radii  for  infin- 
ite kinetics):  as  a result,  the  ratio  of  burning  rate  calculated  using  Fontijn 
kinetics  to  that  calculated  using  infinite  kinetics  asymptotically  approaches 
a constant  value  at  small  particle  sizes,  the  magnitude  of  this  value  depending 
on  ambient  conditions.  For  the  1 atmosphere  cases  studied,  this  asymptotic 
value  is  approached  at  particle  radii  of  5 to  10  microns  (depending  on  ambient 
conditions,  while  for  the  5 atmosphere  case  it  appears  that  it  is  being 
approached  at  about  a 1 to  2 micron  radius,  and  for  the  10  atmosphere  case, 
it  may  not  be  approached  until  particle  radius  is  less  than  1 micron.  The 
burning  ratio  (Fontijn  kinetics/infinite  kinetics)  is  seen  to  decrease  with 
increasing  ambient  temperature  or  increasing  ambient  CC^  mole  fraction. 
Examination  of  Table  III-III  reveals  that  most  of  the  variation  of  the  burning 
rate  ratio  with  ambient  temperature  is  due  to  Increase  in  the  infinite-kinetics 
burning  rate  with  ambient  temperature,  with  the  Fontijn-kinetics  burning  rate 
decreasing  with  increasing  temperature  for  small  particles  (5  micron  radius) 
and  slightly  increasing  with  increasing  temperature  for  large  particles 
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(20  micron  radius) . With  respect  to  the  calculated  dependence  of  burning 
rate  ratio  on  ambient  C02  mole  fraction,  this  is  due  to  the  dependence  on 
this  parameter  of  the  burning  rate  calculated  from  Fontijn  kinetics  being 
somewhat  less  than  that  calculated  using  infinite  kinetics,  both  calculated 
rates  increasing  with  increasing  ambient  COj  mole  fraction.  The  dependence 
of  the  burning  rate  calculated  using  Fontijn  kinetics  on  this  parameter 
decreases  with  decreasing  particle  radius. 

One  of  the  more  critical  assumptions  made  in  development  of  this 
aluminum  particle  burning  model  was  the  assumption  (based  partly  on  thermo- 
dynamic equilibrium  calculations  at  estimated  particle  surface  temperatures) 
that  aluminum  oxide  condensation  at  the  particle  surface  could  be  neglected. 
With  the  outputs  of  the  final  model  (notably  aluminum  vapor  concentration 
and  AI2O  concentration  adjacent  to  the  surface  and  surface  temperature) 
this  assumption  can  be  checked,  at  least  in  terms  of  thermodynamic  equili- 
brium considerations,  more  rigorously.  Such  calculations  have  been  per- 
formed, with  the  predicted  AljO  partial  pressure  at  the  surface  being  compared 
to  the  saturation  Al^O  partial  pressure  as  calculated  from  the  equilibrium 
constant  for  3Al20(g)  Z Al^U)  + 4A1(£.)  at  the  surface  temperature. 

Results  of  calculations  performed  for  the  cases  where  the  condensation 
shell  temperature  was  set  at  3500°K  are  presented  in  Table  ITI-IV.  Similar 
calculations  were  also  performed  for  the  cases  where  the  condensation 
shell  temperature  was  set  at  3000°K.  With  the  3500°K  cases,  thermodynamically 
predicted  surface  condensation  of  A120^  was  only  encountered  for  the  high 
pressure-high  C02  concentration  (P  - 10  atm,  Yqq2  m 0*8)  Case  Number  12, 
while  it  was  marginal  for  large  particles  for  the  high  ambient  temperature 
and  high  pressure  intermediate  C02  concentration  cases.  With  3000°K  as  the 
condensation  shell  temperature,  potential  surface  condensation  problems 
extended  to  more  of  the  high  pressure  cases,  being  predicted  for  the  larger 
particle  sizes  but  not  for  the  smaller,  while  the  largest  particle  sizes 
(30  micron  radius)  gave  problems  for  high  C02  mole  fractions  even  at  1 
atmosphere.  The  trend  that  larger  particles  are  indicated  to  be  most  likely 

to  have  oxide  condensation  on  the  surface  is  quite  interesting  in  that,  as 

(19) 

discussed  earlier,  Prentice  ' observed  that  the  amount  of  initial  aluminum 
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TABLE  III-IV.  Thermodynamic  Evaluation  of  the  Assumption  of  No  Al-O^U) 
Condensation  on  the  Burning  Aluminum  Particle  Surface 


Case  No. 

T , (CK) 

surface 

P*  , (atm) 

AI2O  surface 

P* 

A^O  saturation 

1 

2389 

0.0469 

0.40 

2 

2408 

0.0096 

0.46 

3 

2357 

0.1020 

0.32 

4 

2300 

0.0129 

0.44 

5 

2365 

0.1097 

0.34 

6 

2345 

0.1647 

0.30 

7 

2559 

0.697 

1.20 

8 

2641 

1.571 

1.90 

10 

2368 

0.101 

0.35 

12 

2551 

2.330 

1.18 

13 

2372 

0.078 

0.36 

14 

2348 

0.117 

0.30 

15 

2342 

0.126 

0.29 

16 

2409 

0.0074 

0.46 

17 

2412 

0.0011 

0.48 

18 

2383 

0.044 

0.39 

19 

2370 

0.067 

0.36 

20 

2404 

0.003 

0.46 

21 

2404 

0.002 

0.46 

22 

2405 

0.001 

0.46 

23 

2400 

0.048 

0.44 

24 

2330 

0.162 

0.27 

25 

2381 

0.106 

0.39 

26 

2303 

0.216 

0.22 

27 

2626 

1.695 

1.72 

28 

2668 

1.330 

2.20 

29 

2699 

1.035 

2.65 

30 

2769 

0.202 

0.38 

31 

2302 

0.194 

0.22 

32 

2325 

0.165 

0.26 

33 

2408 

0.029 

0.47 

34 

2416 

0.010 

0.50 

35 

2420 

0.0013 

0.51 

26 

2317 

0.169 

0.24 

37 

2339 

0.139 

0.29 

38 

2378 

0.076 

0.38 

39 

2407 

0.021 

C.46 

40 

2413 

0.008 

0.48 

41 

2415 

0.003 

0.49 

42 

2381 

0.05C 

0.39 

43 

2389 

0.035 

0.40 

44 

2401 

0.013 

0.44 

45 

2406 

0.003 

0.46 

46 

2406 

0.0004 

0.46 

47 

2406 

0.007 

0.46 

*Based  on  equilibrium  calculations  for  SA^CKg)  £ ^2^3^  + 4A1(£) 

at  T - T , 

surface. 
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appearing  In  Che  final  residue  globule  after  burning  decreased  with  decreasing 
particle  size  and,  in  fact,  extrapolation  of  his  data  indicates  that  for 
sufficiently  small  particle  sizes,  no  such  oxide  residue  globule  will  occur. 

The  predicted  burning  rates  as  functions  of  particle  radius 

presented  in  Table  III-III  and  Figures  III-6  and  III-7  for  various  sets  of 

ambient  conditions  were  plugged  Into  Equation  51  for  graphical  integration 

to  yield  predicted  burn  times  as  functions  of  initial  particle  radius. 

Results  of  these  calculations  are  presented  in  Figures  III-8  through  III-13. 

In  each  figure  (representing  one  set  of  ambient  conditions,  as  labeled) 

predicted  burn-time  is  plotted  on  log-log  coordinates  against  initial  particle 

radius.  Lines  representing  three  sets  of  assumptions  are  plotted  on  each 

figuia:  (1)  infinite  kinetics  for  the  Al-CC^  gas-phase  reaction;  (2)  Fontijn 

kinetics,  condensation  shell  temperature  - 3500°K;  and  (3)  Fontijn  kinetics, 

condensation  shell  temperature  * 30008K.  The  average  slope  of  each  line, 

representing  the  exponent  in  the  d n - d°  ■ kt  burning  time  power  law  is 

0 

included  on  these  figures.  While  the  predicted  exponent  for  the  infinite- 

2 

kinetics  case  is  2.0  (classic  diffusion-limited  d law)  the  exponents  for 
the  Fontijn  kinetics  cases  are  somewhat  lower,  ranging  from  approximately 
1.35  to  1.85  for  the  3000°K  condensation  temperature  cases  and  from  1.4  to 
1.9  for  the  3500°K  condensation  cases.  The  predicted  burn  times  for  the 
3500°K  cases  are  anywhere  from  10  to  120  percent  greater  than  for  the  infinite- 
kinetics  cases  (the  largest  differences  occur ing  for  the  smallest  particles) 
while  those  for  the  3000°K  cases  are  from  30  to  200  percent  greater  than  for 
the  infinite-kinetics  cases.  As  would  be  expected,  the  high  pressure  (10  atm) 
cases  exhibit  less  difference  between  burn-time  predictions  with  infinite  and 
finite  rate  klnetcs  than  the  low  pressure  (1  atm)  cases. 

In  Figures  III-14  and  III-15,  predicted  burn-times  are  cross-plotted 
against  ambient  CO^  mole  fraction  for  small  (5p  radius)  and  large  (30y  radius) 
particles  for  each  of  the  three  sets  of  assumptions  listed  above.  As  may 
be  seen,  for  the  small  particles  use  of  the  finite-rate  kinetics  of  Fontijn 
considerably  reduces  the  predicted  dependency  of  burn-time  on  COj  mole  frac- 
tion for  small  particles,  but  does  not  appreciably  affect  this  dependence 
for  large  particles  where  the  burning  is  more  nearly  diffusion-limited. 
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Figure  111-8.  Burn-Time  Versus  Particle  Radius,  T^mb  " tOOO°K,  P ■ 1 atm, 


■ 0.8. 


Figura  111*9.  Burn-Timt  Varsus  Partlcla  Radius,  ■ 1000°K,  P ■ 1 atm,  Y qq  ■ 0.( 
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Figura  IIMO.  Burn-TIma  Varsus  Particla  Radius,  T^mb  ■ 1000°K,  P ■ 1 atm,  Yqq^  " 
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Figure  111-11.  Burn-Time  Versus  Particle  Radius, 


Tamb-  1000%  P -1  atm,  VCo2 
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Figure  111*12.  Burn-Tima  Versus  Particle  Radius,  T^mb  ■ 1000°K,  P»  10  atm,  YqQj  " 
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In  Figures  III-16  and  III-17  predicted  burn-times  are  plotted  against  pressure 
for  small  (5p  radius)  and  large  (20p  radius)  particles.  While  with  infinite 
kinetics  (diffusion-limited  burning)  the  burn-time  is  predicted  to  be 
independent  of  pressure,  with  Fontijn  kinetics  burn-time  is  predicted  to  de- 
crease with  increasing  pressure  (an  expected  result).  Finally,  in  Figures 
III-18  and  111-19,  predicted  burn-time  is  plotted  against  ambient  temperature 
for  5p  and  20p  radius  particles.  In  both  cases,  for  infinite  kinetics,  the 
predicted  burn-time  decreases  with  increasing  ambient  temperature.  However, 
with  Fontijn  kinetics,  for  the  5p  particles,  burn-time  is  predicted  to  increase 
with  increasing  ambient  temperature,  probably  due  to  the  burn  rate  in  these 
cases  being  mainly  determined  by  the  location  of  the  condensation  shell, 
which  moves  out  with  increasing  ambient  temperature.  This  same  trend  is  seen, 
though  to  a much  lesser  extent  for  the  20y  particles  for  the  3000°K  conden- 
sation temperature  case  though  not  for  the  3500°K  case  where  the  effect  is 
sufficient  to  essentially  eliminate  the  decrease  in  burn  time  with  increased 
ambient  temperature  predicted  for  infinite  kinetics,  but  not  to  reverse  it. 

E.  Conclusions 

A model  of  aluminum  particle  combustion  in  CC>2  which  incorporates 
measured  kinetic  data  for  the  gas-phase  reaction  of  aluminum  with  CC^  rather 
than  assuming  a flame-sheet  (infinite  gas-phase  kinetics)  as  in  past  models 
has  been  developed,  programmed  and  used  to  parametrically  study  the  effect 
of  various  parameters  on  particle  burning  rate  and  burn-time.  This  model 
treats  aluminum  oxide  condensation  as  occuring  in  an  infinitesimally  thin 
shell  whose  location  is  determined  by  specification  of  a condensation 
temperature.  Condensation  of  aluminum  oxide  at  the  particle  surface  is  not 
permitted,  but  A10  produced  by  the  gas-phase  reaction  of  aluminum  vapor  with 
CC>2  is  permitted  to  react  with  infinite  kinetics  with  aluminum  liquid  at 
the  surface  to  form  A^O  gas,  thus  removing  aluminum  from  the  surface  as 
well  as  supplying  heat  to  vaporize  additional  aluminum.  Heat  feedback  from 
the  condensation  shell  also  provides  additional  energy  for  aluminum  vapori- 
zation. Two  variations  of  this  model  were  developed:  in  one,  any  CC^ 
reaching  the  surface  was  permitted  to  react  infinitely  rapidly  with  aluminum 
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Figure  111-17.  Burn-Time  Versus  Pressure,  Y^q  * 0.5,  * 1000° K, 
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Figure  111-18.  Bum  Time  Versus  ^AMBIENT’  P * ' atm. 


Figure  111-19.  Burn  Time  Versus  ^AMBIENT’  ^ ' *tmf 


liquid  to  produce  A^O  gas,  while  in  the  other  this  reaction  was  prohibited. 
With  the  similarity  approximations  made  in  the  model (s),  the  former  variant 
leads  to  a degenerate  solution  in  which  the  predicted  burning  rate  is  inde- 
pendent of  the  gas-phase  kinetics,  while  the  latter  variant  leads  to  burning 
rates  predicted  with  the  measured  gas-phase  kinetic  data  which  are  less  than 
those  predicted  with  infinite  kinetics  assumed.  The  ratio  of  burning  rates 
predicted  with  the  measured  kinetic  data  to  those  predicted  with  infinite 
kinetics  for  the  same  conditions  varies  from  about  0.3  to  0.95  over  the  range 
of  conditions  studied,  increasing  with  increased  particle  size,  increased 
pressure,  decreased  ambient  CO^  mole  fraction,  and  decreased  ambient  tempera- 
ture. This  ratio  asymptotically  approaches  unity  for  large  particle  sizes  and 
asymptotically  approaches  a value  between  zero  and  unity,  which  depends  on 
ambient  conditionsjat  very  small  particle  sizes  (where  the  location  of  the 
condensation  shell  dominates  relative  to  gas-phase  kinetics)  . Burn-time  versus 

initial  particle  diameter  calculations  indicate  a reduction  in  the  exponent  in 

2 2 

the  d -burning  law  (dQ  - d * kt)  from  2 for  the  case  of  infinite  kinetics 
(diffusion-limited  burning)  to  1.35  - 1.90,  depending  on  ambient  conditions. 
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